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Abstract 

Lithium  tetraborate  (Li2B407  or  LTB)  is  a  promising  material  for  both  radiation 
dosimetry  and  neutron  detection  applications.  LTB  crystals  can  be  grown  pure  or  doped 
with  different  impurities  including  transition-metal  and  rare-earth  ions.  Research  in  this 
dissertation  focuses  on  undoped  LTB  crystals  and  LTB  crystals  doped  with  copper  and 
silver.  Electron  paramagnetic  resonance  (EPR)  and  electron-nuclear  double  resonance 
(ENDOR)  are  used  to  characterize  point  defects  in  the  lithium  tetraborate  crystals. 
Thermoluminescence  (TL),  photoluminescence  (PL),  photoluminescence  excitation 
(PLE),  and  optical  absorption  (OA)  are  also  used.  An  intrinsic  hole  trap  associated  with 
lithium  vacancies  is  characterized  with  EPR  and  ENDOR  and  its  thermal  stability  is 
determined  using  thermoluminescence.  A  “perturbed”  hole  trap  due  to  Ag2+  ions  is 
characterized  in  doped  crystals  using  EPR  data  alone.  This  method  is  tested  on  a 
previously  studied  hole  center  where  both  EPR  and  ENDOR  were  used.  New  x-ray 
induced  centers  are  identified  in  copper-doped  crystals.  These  include  two  Cu2+  hole 
centers  and  two  Cu°  electron  centers.  These  centers  are  characterized  with  EPR  and  their 
thermal  stability  explains  TL  peaks  in  glow  curves.  Finally,  a  comprehensive  study 
utilizing  EPR,  OA,  PL,  and  PLE  data  provide  convincing  explanations  for  the  absorption 
and  emission  features  of  silver-doped  crystals. 
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EPR  AND  ENDOR  STUDIES  OF  POINT  DEFECTS  IN  LITHIUM 

TETRABORATE 


I.  Introduction 


1.1.  Motivation 

Since  the  early  1990s,  United  States  national  defense  strategy  has  shifted  from  the 
cold-war  era  threat  of  the  Soviet  Union  to  an  asymmetric  threat  environment.  One  of 
these  threats  has  been  the  rise  in  terrorism.  Terrorist  groups,  namely  al-Qaeda,  have 
shown  an  ability  to  be  clever  and  furtive  in  their  plots  to  attack  the  US  and  its  allies. 
Terrorist  groups  have  shown  interest  in  seeking  and  deploying  weapons  of  mass 
destruction  and  mass  disruption— weapons  that  might  kill  no  one  but  would  create 
widespread  psychological  trauma.  In  testimony  before  the  Senate  Foreign  Relations 
Committee  on  March  6,  201 1,  some  of  the  nation’s  top  nuclear  officials  said  “dirty 
bombs”  could  become  the  terrorist’s  weapon-of-choice.  A  dirty  bomb,  also  called  a 
radiological  weapon,  would  use  conventional  explosives  to  spread  radioactive  materials. 
Likely  radioactive  elements  for  these  devices  would  be  cesium,  cobalt,  and  iridium 
isotopes  since  they  are  widely  used  for  industrial  purposes  and  are  easily  acquired.  These 
materials  could  also  be  used  in  a  dirty  bomb  or,  even  worse,  a  nuclear  bomb.  In  June 
2011,  conventional  bombs  disguised  in  computer  ink  cartridges  were  timed  to  blow  up 
aboard  commercial  aircraft  while  flying  over  the  US.  The  material  in  these  devices  could 
easily  have  been  “laced”  with  radioactive  material  to  magnify  the  effect.  It  is  a  national 
priority  to  detect  clandestine  radioactive  materials  before  they  can  be  transported  and 
exploited  by  terrorists  within  the  US.  Border  officials  and  other  security  forces  must 
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possess  the  proper  equipment  to  be  able  to  detect  these  unauthorized  radiological 
materials. 

Radioactive  decay  is  a  process  where  an  unstable  nucleus  loses  energy  by  emitting 
particles  and  high  energy  photons.  Typical  decay  processes  involve  the  emission  of  alpha 
and  beta  particles  and  gamma  rays.  There  are  commercial  instruments  that  are  capable  of 
detecting  these  decay  products,  but  shielding  reduces  the  effectiveness  of  these  devices  in 
detecting  radioactive  decays.  On  the  other  hand,  neutron  emission,  which  is  also  a 
product  of  radioactive  decay  and  spontaneous  fission  can  be  used  as  an  indicator  of  the 
presence  of  radioactive  material.  Thus,  their  detection  can  be  an  indication  of  the 
presence  of  nearby  radioactive  material.  Because  neutrons  have  no  charge,  detecting 
them  needs  to  be  done  through  a  reactionary  process. 

Currently,  neutron  detectors  are  commercially  available,  but  they  are  neither  compact 
nor  cost  effective.  He-3  detectors  and  BF3  tubes  are  two  examples.  Some  detectors 
operate  by  monitoring  visible  luminescence  when  neutron  capture  occurs  [1].  They  are 
most  efficient  in  detecting  neutrons  with  energy  of  30  meV  or  less  [2].  A  typical  fissile 
decay  from  radiological  material  emits  neutrons  of  1.5  MeV,  so  detectors  use  thick 
moderators  to  slow  the  neutrons.  These  moderators,  however,  reduce  the  flux  of 
neutrons.  Some  current  detectors  also  need  large  chambers  and  high-voltage  power 
sources  which  make  them  bulky  and  less  portable. 

Solutions  to  the  limitations  of  current  neutron  detectors  are  emerging.  One  possibility 
is  to  use  a  luminescent  material  that  has  a  high  neutron  capture  cross  section.  Several 
isotopes  have  been  experimentally  found  to  have  high  neutron  capture  cross  sections: 

-  Gd,  B,  ’He,  and  Li.  Their  capture  cross  sections,  specified  in  bams,  are  plotted  as  a 
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function  of  neutron  energy  in  Fig.  1.1.  The  capture  cross  sections  of  these  nuclei  are 
orders  of  magnitude  larger  than  most  other  nuclei  at  the  meV  (thermal  neutron)  range. 
The  present  dissertation  focuses  on  single  crystals  of  lithium  tetraborate,  a  material  that 
has  concentrations  of  lithium  and  boron  nuclei  and  thus  is  potentially  useful  for  neutron 
detection. 


Figure  1.1.  Neutron  capture  cross-sections  of  various  isotopes  [3]. 


Lithium  tetraborate  (LLB4O7),  commonly  referred  to  as  LTB,  is  a  versatile  insulating 
crystal  with  potential  use  in  neutron  detection.  LTB  can  be  grown  using  enriched  6Li  and 
10B  to  increase  its  neutron  opacity  [4].  When  doped  with  an  appropriate  optically  active 
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element,  LTB  could  be  used  as  an  effective  scintillation  material.  Candidate  elements  are 


cerium  or  erbium  because  of  their  excellent  luminescent  properties  in  other  crystals  [5]. 
Lithium  tetraborate,  has  been  grown  pure  and  doped  with  many  different  elements 
including  transition  metals,  actinides,  and  rare  earth  elements.  The  research  in  this 
dissertation  addresses  undoped  lithium  tetraborate  as  well  as  LTB  doped  with  copper  and 
silver. 

The  following  nuclear  reactions  for  lithium  and  boron  involve  neutron  capture. 


+  n  -»  7Li(0.84  MeV)  +  4He  (1.47  MeV)  +  y  (0.48  MeV) 

(94%) 

(1.1) 

+  n  -»  7Li(1.01  MeV)  +  4He  (1.78  MeV) 

(6%) 

(1.2) 

6Li  +  n  -»  3H( 2.73  MeV)  +  4He  (2.05  MeV) 

(1.3) 

In  these  reactions,  both  lithium  and  boron  isotopes  absorb  neutrons  and  then  decay  by 
releasing  alpha  particles  and  lithium  or  tritium  nuclei,  along  with  energy.  Since  10B  only 
has  a  20%  natural  abundance  and  6Li  has  an  even  smaller  7%  abundance,  materials  need 
to  be  enriched  with  10B  and  6Li  to  give  them  significantly  higher  neutron  capture  cross 
sections. 

Consider  an  enriched  LTB  crystal.  If  a  thermal  neutron  is  captured  by  a  lithium  or 
boron  atom  in  the  crystal  lattice,  an  energetic  alpha  particle  would  be  released  according 
to  Eqs  1.1  -  1.3.  Alpha  particles  will  not  travel  far,  on  the  order  of  microns  in  pure 
silicon  for  example  [2],  During  this  time,  the  positive  charge  of  the  high  speed  alpha 
particle  will  attract  electrons  and  cause  them  to  leave  their  lattice  bonds  and  move  about 
the  crystal  lattice  as  conduction  electrons.  Many  of  these  electrons  will  absorb  energy 
from  the  alpha  particles.  The  alpha  particles  will  thus  lose  kinetic  energy  and  eventually 
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combine  with  two  electrons  to  form  helium  atoms.  The  conduction  electrons  formed  may 
be  detected  indicating  neutron  capture. 

As  alpha  particles  travel  in  the  crystal  lattice,  there  are  two  features  that  could  aid  in 
neutron  detection.  First,  electrons  that  absorb  energy  from  the  alpha  particles  become 
conducting  electrons  in  the  crystal.  If  the  material  were  biased  with  a  voltage,  a  change 
in  the  resulting  current  could  be  measured  indicating  neutron  absorption.  This  principle 
is  the  basis  of  solid  state  detectors.  Second,  many  scintillators  emit  photons  when  the 
crystal  is  excited  by  energy  depositions.  In  scintillating  materials,  “excited”  electrons 
recombine  with  “holes”  in  the  crystal  lattice  and  release  electromagnetic  radiation  that 
can  be  detected  using  a  photomultiplier  tube.  If  a  neutron  absorbing  material  such  as 
lithium  tetraborate  is  attached  to  a  photomultiplier  tube  and  placed  in  the  vicinity  of 
radioactive  materials,  alpha  particles  would  be  created  when  neutrons  react  with  boron 
and  lithium  atoms  in  the  crystal.  These  alpha  particles  will,  in  turn,  create  electrons  in 
the  conduction  band.  These  electrons  could  participate  in  radiative  recombination 
processes  in  the  lattice,  i.e.,  scintillate,  and  produce  a  measureable  signal  in  a 
photomultiplier  tube. 

To  explain  scintillation  (shown  in  Fig.  1.2),  first  consider  that  semiconducting  and 
insulating  materials  have  a  band  gap,  or  forbidden  energy  region.  The  energy  levels 
below  this  gap  are  called  the  valence  band  and  represent  those  electrons  that  are  bound  to 
lattice  sites  or  form  bonds  between  ions.  When  an  electron  gains  energy,  perhaps  from 
ionizing  radiation  (See  step  A  in  Fig.  1.2),  it  can  be  excited  across  this  gap  into  the 
conduction  band  (See  step  B  in  Fig.  1.2).  Electrons  in  the  conduction  band  are  no  longer 
localized  and  have  enough  energy  to  move  about  the  crystal  lattice.  The  missing  electron 
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Electron  in  conduction  band 


Figure  1.2.  Energy  band  structure  illustrating  scintillation. 


in  the  valence  band  is  called  a  hole,  where  holes  can  conduct  about  the  crystal  lattice  as 
well.  Many  materials  contain  impurities  that  are  present  either  intentionally  or 
unintentionally.  Impurities  may  allow  energy  levels  to  exist  in  the  band  gap.  These 
impurities  are  sometimes  called  activators.  A  conducting  electron  may  fall  into  an 
activator  site  after  first  giving  up  some  of  its  energy,  i.e.,  through  a  phonon  interaction. 
From  the  activator  location,  the  electron  then  emits  energy  via  a  photon  of  energy  that  is 
less  than  the  band  gap  of  the  material  and  returns  to  the  valence  band  (See  step  C  in  Fig. 
1.2).  In  other  words,  the  electron  recombines  with  a  hole  in  the  valence  band  and  thus 
returns  the  crystal  to  its  pre-irradiated  equilibrium  state. 

Once  an  electron  has  returned  to  the  valence  band  via  an  acceptor  site,  the  emitted 
photon  could  be  measured  by  converting  it  to  an  electron  using  the  photoelectric  effect. 
Photomultiplier  tubes  are  designed  to  first  convert  incident  photons  into  a  current  of 
electrons.  These  tubes  then  multiply  the  electrons  produced  millions  of  times  by  the 
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process  of  secondary  emission.  The  efficiency  of  current  production  in  a  photomultiplier 
tube  depends  on  the  wavelength  of  incident  photons.  Photomultiplier  tubes  are  most 
efficient  in  the  visible  spectrum.  With  this  very  high  sensitivity,  a  photomultiplier  tube 
could  be  combined  with  a  properly  doped  lithium  tetraborate  crystal  known  to  scintillate 
in  the  visible  spectrum  and  become  an  inexpensive  and  compact  solution  in  detecting 
neutron  radiation. 

Electrons  can  release  energy  through  processes  other  than  scintillation  in  order  to 
return  to  the  valence  band.  Various  unintentional  impurities  and  dislocations  create 
energy  levels  within  the  band  gap  where  nonradiative  recombinations  can  occur.  For  this 
reason,  unintentional  contamination  or  poor  crystal  growth,  i.e.,  crystals  that  do  not  grow 
uniformly,  must  be  avoided  or  corrected  with  post  growth  techniques  such  as  annealing. 
Materials  are  commonly  “doped”  with  specific  impurities  that  allow  energy  levels  to  exist 
inside  the  band  gap  and  to  increase  the  probability  of  scintillation. 

In  order  to  elaborate  further  on  lithium  tetraborate’s  potential  as  a  scintillating 
material,  consider  the  alpha  particles  produced  in  the  reactions  described  in  Eqs.  1.1  -  1.3. 
These  particles  range  in  energy  from  1.47  to  2.05  MeV.  For  a  wide  category  of  materials, 
energy  on  the  order  of  three  times  the  band  gap  energy  is  required  to  create  one  electron- 
hole  pair  [2].  If  an  alpha  particle  of  1.47  MeV  gave  up  two  thirds  of  its  energy  in  LTB 
with  an  approximate  10  eV  band  gap,  32,700  electron-hole  pairs  would  be  produced. 
Various  experimental  determinations  have  shown  that  the  absolute  scintillation  efficiency 
of  thallium-activated  sodium  iodide  is  about  12%.  If  LTB  were  only  8%  efficient,  an 
alpha  particle  would  provide  78,400  eV  in  total  light  energy,  or  28,500  photons  that  have 
a  450  nm  wavelength.  At  8%  efficiency,  LTB  would  still  have  one  emitted  photon  for 
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each  electron  hole  pair  created  making  it  a  good  candidate  as  a  scintillating  material  to  be 
used  in  neutron  detection.  Also,  a  thin  sample  of  LTB  with  its  low  Z  would  make  it 
nearly  gamma  blind,  i.e.,  gamma  radiation  would  not  likely  produce  additional  electrons 
in  the  conduction  band  because  the  gammas  would  pass  through  the  thin  sample  without 
interacting  with  the  crystal  nuclei. 

LTB  has  other  valued  commercial  uses.  It  has  been  applied  to  dosimetry  and 
scintillation  detectors  [6-7]  and  also  has  potential  uses  in  the  fields  of  acousto-electronics, 
and  optics  [6,  8-10].  LTB  exhibits  piezoelectric  coupling  values  that  fall  between  those 
of  lithium  niobate  and  quartz,  and  possesses  orientations  for  which  the  temperature 
coefficient  of  frequency  and  delay  time  is  zero  for  bulk  and  surface  acoustic  waves  [11]. 
Despite  LTB’s  potential  as  either  a  primary  or  companion  material  in  future  electronic 
devices,  there  is  still  a  lack  of  fundamental  characterization,  particularly  regarding  the 
defect-related  properties. 

The  ideal  lithium  tetraborate  crystal  belongs  to  the  space  group  I4icd  (tetragonal  I4cd, 
the  point  group  is  4mm)  and  possesses  104  atoms  per  unit  cell.  LTB  has  lattice  constants 
a  =  9.475  A  and  c  =  10.283  A  at  room  temperature  [12-16].  On  the  left  side  of  Fig.  1.3,  a 
basic  (B4O9)6  building  block — one  of  eight  as  part  of  the  conventional  unit  cell — is 
displayed  along  with  a  neighboring  lithium  ion.  The  c  axis  is  oriented  upward  and 
toward  the  left  (see  the  right  side  of  Fig.  1.3).  A  tetragonal  structure  like  LTB  has  two  a 
axes  (a  =  b  in  Fig.  1.3)  due  to  its  symmetry.  The  full  unit  cell,  on  the  right  hand  portion 
of  Fig.  1.3,  is  meant  to  illustrate  the  complexity  of  the  crystal.  The  left  portion  of  Fig.  1.3 
shows  the  basic  building  block  of  an  LTB  crystal  structure  with  the  appropriate 
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Figure  1.3:  Basic  structural  component  of  the  LEB4O7  crystal  lattice  (left)  and  a  unit  cell 

of  Li2B407  (right). 

equivalent  atoms  labeled.  The  extra  02,  03,  and  04  on  the  lithium  are  additional  to  the 
basic  structure.  They  are  shown  to  aid  visualizing  the  other  lithium  atoms  in  the  lattice. 

All  crystals  have  point  defects  that  occur  during  growth.  These  defects  often  affect 
the  bulk  properties  of  solids  such  as  sub-band-gap  optical  absorption  or  emission  and 
affect  the  mechanical  and  electrical  properties  as  well.  For  example,  the  F  center  in 
potassium  chloride  consists  of  an  electron  trapped  in  a  chloride  ion  vacancy.  These  F 
centers  cause  the  normally  colorless  potassium  chloride  crystals  to  appear  dark  blue. 

Electron  paramagnetic  resonance  (EPR)  and  electron-nuclear  double  resonance 
(ENDOR)  are  powerful  techniques  in  determining  defect  structure  [17].  Their  restriction 
is,  of  course,  that  the  defects  must  be  paramagnetic.  Fortunately,  EPR  has  demonstrated 
that  an  x-ray  irradiated  sample  of  undoped  LTB  produces  spectra  due  to  oxygen  and 
lithium  vacancies.  The  lithium  vacancies  trap  holes  and  are  the  subject  of  Chapter  Three. 
Further,  there  are  x-ray  irradiated  spectra  produced  in  both  silver  and  copper  doped  LTB 
crystals.  EPR  spectra  are  also  present  in  copper  doped  LTB  crystals  without  the  need  for 
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irradiation.  Chapters  Four  and  Five  describe  the  EPR  behavior  of  silver  and  copper  in 
LTB.  Chapter  Six  presents  the  optical  properties  of  the  silver  defects  characterized  with 
EPR  in  Chapter  Four. 

The  characterization  of  point  defects  in  lithium  tetraborate  that  are  described  in  this 
dissertation  will  aid  in  developing  devices  utilizing  this  material.  These  studies  may  also 
help  clarify  crystal-growth  properties  and  suggest  improvements.  For  example,  the  large 
concentrations  of  lithium  vacancies  in  LTB  (discussed  in  Chapter  Three)  could  possibly 
be  reduced  by  heating  crystals  close  to  their  melting  point  while  immersed  in  a  lithium- 
containing  powder.  This  process,  known  as  vapor  phase  equilibration,  has  been 
successfully  used  to  eliminate  lithium  vacancies  in  lithium  tantalate  and  lithium  niobate 
crystals  [18]. 
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II.  Experimental  Techniques 


2.1.  Electron  Spin 

O.  Stem  and  W.  Gerlach,  in  1922,  used  a  beam  of  silver  atoms  to  demonstrate  that  the 
magnetic  moment  of  an  electron  could  only  have  discrete  orientations  in  a  magnetic  field. 
In  their  experiment,  silver  atoms  passed  through  a  region  of  magnetic  field  and  struck  a 
target  placed  a  short  distance  beyond  the  magnet.  A  continuous  line  along  the  horizontal 
axis  of  their  target,  due  to  the  slit  shape  of  the  collimator,  was  predicted  for  the  silver 
atoms  since  they  do  not  possess  orbital  angular  momentum,  i.e.,  L  =  0,  and  thus  were 
assumed  to  have  zero  magnetic  moment.  To  Stern  and  Gerlach’ s  surprise,  the  atoms 
struck  the  target  in  two  specific  regions  above  or  below  their  predicted  horizontal  line. 
Later  in  1925,  G.  Uhlenbeck  and  S.  Goudsmit  explained  Stem  and  Gerlach’s  experiment 
by  assuming  that  the  electron  possessed  an  intrinsic  magnetic  moment  that  was  due  to  an 
intrinsic  angular  momentum,  thus  the  creation  of  a  quanta  of  spin  [19].  Silver,  with  one 
unpaired  electron,  has  a  spin  of  S  =  1/2.  This  intrinsic  angular  momentum  is  called  the 
spin  of  an  electron.  In  a  magnetic  field,  the  electron  spin  can  only  take  on  discrete 
values,  whereas  in  the  absence  of  a  magnetic  field,  the  energies  of  the  different  spin  states 
are  degenerate.  The  application  of  a  magnetic  field  lifts  this  degeneracy  and  splits  the 
energy  levels.  The  split  in  the  energy  levels  of  an  electron  in  the  presence  of  a  magnetic 
field  is  called  the  Zeeman  splitting. 

Single  unpaired  electrons  have  a  spin  value  of  S  =  1/2,  which  can  take  on  discrete 
values  of  ms  =  ±  1/2.  In  general,  a  spin  system  with  a  net  spin  S  can  take  on  the 
quantized  values  ms  =  -S,  -S  +1,  -S+2  to  +S.  The  Zeeman  energy  associated  with  the 
magnetic  field  splitting  is  given  by, 
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(2.1) 


where  g  is  the  magnetic  moment.  If  the  spin  states  were  placed  in  a  magnetic  field 
chosen  to  be  the  z  direction,  the  corresponding  Zeeman  energy  levels  would  be  given  by 

E  =  -gz  Bz  (2.2) 

where  gz  =  -ge  /?  ms  (2.3) 

The  z  direction  is  chosen  to  be  in  the  direction  of  the  magnetic  field  in  this  chapter  to 
avoid  unnecessary  complications  arising  from  vector  notation  and  dot  products.  The 

6  ft  94  2 

quantity  /?  is  the  Bohr  magneton,  - or  9.274015  x  10""  Am",  and  ge  is  the  electron 
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spin  g-factor  (more  often  simply  called  the  electron  g-factor).  The  electron  g  factor  has  a 
value  of  2.002319  for  a  free  electron  and  is  known  to  extremely  high  precision  [20].  The 
energy  given  in  Eq.  2.2  is  directly  proportional  to  the  magnetic  field.  The  Zeeman 
splitting  of  a  single  electron,  with  S  =  1/2,  is  demonstrated  in  Fig.  2.1.  In  this  figure,  the 
separation  between  the  two  ms  energy  levels  increases  linearly  as  the  magnitude  of  the 
magnetic  field  increases.  The  transition  energy,  AE,  is  the  energy  difference  between  the 
two  spin  states  of  the  electron.  A  transition  between  these  energy  levels  can  be  induced 
when  an  amount  of  energy,  AE  =  hv ,  is  provided  by  an  external  electromagnetic  wave 
oriented  with  its  magnetic  field  perpendicular  to  the  static  applied  magnetic  field.  This 
energy  is  absorbed  by  the  electron  and  causes  it  to  switch  spin  states,  going  from  -1/2  to 
+1/2.  The  electrons  can  then  “relax”  back  by  giving  energy  up,  usually  by  emitting 
phonons.  The  basis  of  an  EPR  experiment  involves  inducing  a  change  in  the  spin  state  of 
an  unpaired  electron  and  detecting  the  absorption  of  energy  that  takes  place. 
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Figure  2.1.  Illustration  of  Zeeman  splitting  for  an  electron  in  a  magnetic  field. 


Unpaired  electrons  in  crystals,  free  radicals,  and  many  other  paramagnetic  entities 
have  electron  g  factors  that  differ  from  the  free-electron  value.  The  difference  in  the  g 
factor  is  due  to  the  environment  in  the  vicinity  of  the  electron.  This  can  cause  shifts  and 
splittings  in  the  Zeeman  energy  levels.  The  reasons  for  these  shifts  and  splittings  in 
crystals  will  be  discussed  later  in  this  chapter.  A  spin  Hamiltonian  is  used  to  represent 
the  spin-related  energy  of  a  paramagnetic  system.  The  discrete  energy  levels  associated 
with  an  electron  spin  can  be  measured  using  electron  paramagnetic  resonance  (EPR) 
spectroscopy.  These  energies  that  an  EPR  experiment  measures  represents  the  diagonal 
elements  of  the  modeled  spin  Hamiltonian.  From  these  energies,  the  parameters 
describing  the  spin  Hamiltonian  can  be  determined  using  a  least-squares  fitting  process. 
The  spin  Hamiltonian  is  often  too  complex  to  fit  without  the  aid  of  a  computer.  These 
extracted  spin  Hamiltonian  parameters  provide  information  about  the  local  spin 
environment  and  are  useful  in  establishing  defect  models. 
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2.2.  Electron  Paramagnetic  Resonance  (EPR) 

Electron  paramagnetic  resonance  spectroscopy  can  be  used  to  determine  the  local 
environment  surrounding  an  unpaired  electron.  Electron  paramagnetic  resonance  (EPR), 
or  as  it  is  sometimes  called  electron  spin  resonance  (ESR),  is  a  technique  for  studying 
point  defects  in  inorganic  crystals  and  amorphous  solids,  acceptors  and  donors  in 
narrow-band-gap  and  wide-band-gap  semiconductors,  and  chemical  species  that  have 
one  or  more  unpaired  electrons  such  as  organic  and  inorganic  free  radicals  and  spin 
labels.  EPR  was  first  observed  at  Kazan  State  University  by  Soviet  physicist  Yevgeny 
Zavoisky  in  1944,  and  was  developed  independently  following  World  War  II  by  Brebis 
Bleaney  at  the  University  of  Oxford  [21]. 

In  a  simplified  description,  EPR  works  by:  (1)  placing  a  paramagnetic  material  in  a 
static  magnetic  field,  (2)  applying  a  microwave  frequency  (n)  with  the  microwave 
magnetic  field  perpendicular  to  the  static  magnetic  field  such  that  the  spin  states  will 
undergo  a  transition  due  to  the  absorption  of  energy,  A E  —  hv,  and  finally  (3)  detecting 
this  absorption.  The  g  factor  can  then  be  determined  by  solving  the  following  resonance 
equation. 

A E  =  hv  =  gzp  Bz  (2.4) 

Conceptually,  the  simplest  way  to  conduct  an  EPR  experiment  would  be  to  hold  the 
static  magnetic  field  fixed  at  some  value  and  then  sweep  through  a  range  of  frequencies 
while  monitoring  for  the  absorption  of  energy.  This  is  how  many  optical  spectrometers 
operate.  However,  sweeping  the  frequency  is  not  feasible  in  an  EPR  spectrometer 
because  that  would  require  changing  the  physical  dimensions  of  the  microwave  cavity  as 
the  sweep  progresses.  (A  microwave  cavity  only  resonates  at  one  frequency  for  a  given 
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set  of  dimensions.)  In  practice,  the  static  magnetic  field  is  slowly  swept  in  EPR 
spectrometers  while  the  microwave  frequency  is  held  constant. 

EPR  spectrometers  typically  operate  with  magnetic  fields  up  to  1.3  Tesla  (13,000 
gauss).  For  a  static  magnetic  field  value  of  3,400  gauss  and  a  g  value  of  2.0023,  the 
necessary  frequency  would  be  9.5284  GHz,  which  is  in  the  X-band  for  microwaves  and  is 
the  frequency  where  many  EPR  spectrometers  operate.  As  the  magnetic  field  is  swept 
through  a  region  of  interest,  the  spin  system  will  go  through  one  or  more  resonances 
where  microwave  energy  is  absorbed.  The  output  of  an  EPR  spectrometer  is  absorption 
versus  magnetic  field.  In  actual  experiments,  the  output  of  the  spectrometer  does  not 
show  absorption  curves,  but  rather  their  first  derivatives.  This  is  because  of  an  applied 
modulation,  usually  with  a  frequency  of  100  kHz  that  produces  an  AC  contribution  to  the 
static  magnetic  field.  This  modulation,  when  coupled  with  a  phase-sensitive  detector, 
dramatically  increases  the  sensitivity  of  the  EPR  spectrometer.  Descriptions  of  how  this 
modulation  process  produces  the  first  derivative  can  be  found  in  literature  [22].  The 
output  of  an  EPR  spectrometer  is  thus  the  first  derivative  of  the  absorption  curve,  as 
illustrated  in  Fig.  2.2. 

EPR  spectra  are  usually  more  complicated  than  just  a  single  derivative  spectrum. 
Electrons  interact  with  other  magnetic  moments  or  electric  quadrupole  moments  that  are 
part  of  the  same  atom  or  are  nearby  and  this  causes  the  spin  related  energy  levels  to  split 
and  possibly  shift.  When  an  electron  interacts  with  a  nearby  nucleus  that  has  a  non-zero 
spin,  the  electron  Zeeman  levels  will  be  split.  This  is  called  the  hyperfine  interaction. 
Combining  Eqs.  2.2  and  2.3,  electrons  have  Zeeman  energy  levels  described  by 
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Figure  2.2.  An  EPR  absorption  curve  (top)  and  its  corresponding  first  derivative 

(bottom)  [23]. 


£(ms)  =  g/3msBz  (2.5) 

where 

£(+;)  =~2sPBz 
)  =  -\sPB2 

When  an  electron  spin  S  =  1/2  interacts  with  a  nuclear  spin  I  =  1/2,  the  energy  is 
described  by 

^(ms,  irq)  =  g  /?  ms  Bz  +  /iims  m;  (2.6) 

where  the  second  term  in  Eq.  2.6  is  called  the  hyperfine  term  or  hyperfine  interaction. 

The  four  possible  energy  levels  of  this  S  =  1/2, 1  =  1/2  spin  system  are  now  described 
using  E(ms,  mi). 
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f(+2'“2j  =  -29PB2--hA 

/  1  1\  1  1 

Errr-2^8!+i'1J1 

/  1  1\  1  1 

ff(“2'  +  2)= 

EPR  transitions  occur  when  Ams  =  ±1  and  have  the  restriction  that  Ami  =  0.  Figure  2.3 
describes  these  energy  levels  due  to  the  hyperfine  splitting  with  the  possible  EPR 
transitions  shown  using  red  arrows.  The  resulting  spectrum  corresponds  to  two  EPR 
transitions  centered  about  the  electron’s  g-factor  and  separated  by  A  (the  hyperfine 
parameter).  The  two  peaks  shown  in  Fig.  2.3  occur  when 

/I  1\  /  1  1\  1  1 

A  El  =  hv  =  E[-,-)-E[--.-)=-gpBlz+-hA 

and 

/I  1\  /  1  1\  1  1 

AE2  =  2iv  =  E  —  i;  17  /?  B2z  —  —  h  A 

If  the  hyperfine  interaction  is  large  and  the  EPR  linewidths  are  narrow,  a 
hyperfine  parameter  can  be  determined  directly  from  an  EPR  spectrum.  If,  however,  the 
hyperfine  interaction  is  small  and  the  EPR  linewidths  are  large,  the  hyperfine  splitting  is 
not  resolved  in  an  EPR  spectrum.  In  this  latter  case,  electron-nuclear  double  resonance 
(ENDOR)  may  be  used  to  determine  the  hyperfine  parameter.  ENDOR  will  be  discussed 
later  in  this  chapter.  Figure  2.3  is  the  simplest  case  of  a  hyperfine  interaction  (an  electron 
of  spin  1/2  interacting  with  a  nucleus  of  spin  1/2).  An  EPR  spectrum  becomes  more 
complicated  for  systems  that  have  larger  electron  and  nuclear  spin  values.  The  number 
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Resulting  EPR  Spectrum 


Figure  2.3.  Energy  level  splitting  due  to  an  I  =  1/2  hyperfine  interaction  (upper).  The 

resulting  EPR  spectrum  (lower). 


of  spin-related  energy  levels  is  (2S  +  1)- (21  +  1)  and  the  number  of  EPR  transitions  is 
(2S)-(2I+1),  where  the  selection  rules  Ams  =  ±1  and  Ami=  0  have  been  invoked.  For 
example,  an  Fe  ion  has  5  unpaired  electrons  (S  =  5/2)  and  has  an  I  =  1/2  nuclear  spin 
value.  This  would  result  in  twelve  level  energy  diagram  with  ten  possible  EPR 
transitions.  So,  rather  than  observing  a  large  single  EPR  transition  in  the  spectrum,  there 
would  be  ten  transitions  (or  lines)  with  the  same  total  area  under  them  as  a  large  single 
transition.  For  this  reason  low  concentration  spin  systems  with  large  S  and  I  values  may 
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have  EPR  signals  with  low  intensities  that  do  not  significantly  rise  above  the  background 
noise.  There  are  other  interactions  that  further  complicate  an  EPR  spectrum.  This  will  be 
discussed  in  the  section  titled  “Spin  Hamiltonian.” 

2.3.  EPR  Spectrometer 

An  EPR  spectrometer  collects  data  by  using  a  microwave  “bridge,”  schematically 
shown  in  Fig.  2.4.  The  role  of  each  major  component  in  the  diagram  will  be  discussed 
and  related  to  the  operation  of  the  EPR  spectrometer.  To  begin,  the  microwave  source  is 
located  in  the  lower  left  portion  of  Fig.  2.4.  Historically,  klystrons  produced  the 
microwaves  but  present  day  instruments  use  solid  state  Gunn  diodes  as  the  microwave 
source.  After  leaving  the  Gunn  diode,  the  waveguide  branches  into  two  arms.  The  signal 
arm  contains  the  microwaves  that  continue  down  to  the  resonant  cavity  containing  the 
sample  to  be  studied.  There  is  an  attenuator  in  the  signal  arm,  shown  by  a  variable 
resistor  symbol,  that  controls  the  amount  of  microwave  energy  reaching  the  sample 
cavity.  The  actual  attenuator  is  a  lossy  dielectric  material  that  is  moved  in  and  out  of  a 
section  of  the  waveguide.  In  many  EPR  experiments  it  is  important  to  be  able  to  vary  the 
microwave  power  level  with  the  attenuator  because  the  EPR  signal  may  “saturate”  at 
higher  power  levels.  Saturation  occurs  when  the  microwaves  are  too  intense  upon  the 
sample  cause  the  population  of  the  upper  and  lower  energy  states  to  become  nearly  equal; 
this  results  from  a  long  spin-lattice  relaxation  time  that  prevents  spins  in  the  upper  state 
from  quickly  relaxing  to  the  lower  state.  The  power  level  of  the  microwaves  is  a  crucial 
element  in  collecting  good  EPR  spectra  and  may  need  to  be  adjusted  when  looking  at 
different  defects  in  the  same  crystal  or  different  orientations  of  the  same  defect. 
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Microwave  Bridge 


Figure  2.4.  EPR  Spectrometer  diagram  [23]. 


Next,  the  microwaves  in  the  signal  arm  continue  from  the  attenuator  to  the  three -port 
circulator,  a  ferromagnetic  device.  Entering  from  port  1,  the  microwaves  continue  out 
port  2  and  on  to  the  resonant  cavity  containing  the  sample.  The  cavity  is  “tuned”  such 
that  the  microwaves  are  all  absorbed  (in  other  words,  the  cavity  perfectly  terminates  the 
waveguide).  Effectively,  the  cavity  is  impedance  matched  to  the  waveguide.  The 
microwave  cavity  is  said  to  be  “critically  coupled”  when  there  is  no  reflected  microwave 
power.  Absorption  of  microwaves,  which  occurs  when  the  magnetic  field  satisfies  the 
resonance  condition,  causes  the  microwave  cavity  to  become  slightly  uncoupled  and  thus 
reflects  microwave  power  back  up  the  waveguide  and  into  port  2  of  the  circulator.  These 
reflected  microwaves  from  the  cavity  then  pass  out  of  port  3  of  the  circulator  and 
continue  on  to  the  detector  diode.  The  detector  diode  is  a  semiconductor  junction  that 
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creates  current  proportional  to  the  incident  microwave  power  and  produces  the  measured 
EPR  signal. 

Separately,  there  is  the  reference  arm  of  the  microwave  bridge.  This  portion  of  the 
waveguide  performs  two  functions.  First,  it  provides  microwave  energy  that  serves  to 
“bias”  the  detector  diode.  This  “bias”  energy  incident  on  the  detector  diode  does  not 
come  from  the  cavity  containing  the  sample,  but  instead  comes  directly  from  the  Gunn 
diode  through  the  reference  arm  to  the  detector  diode.  The  purpose  of  the  detector  bias  is 
to  increase  the  spectrometer’s  detection  sensitivity.  In  an  EPR  spectrometer,  the  detector 
is  a  Schottky  barrier  diode  that  can  operate  either  in  the  square  law  region  (where  the 
electrical  power  is  proportional  to  the  square  of  voltage  or  current)  or  in  the  linear  region. 
The  square  law  region  occurs  at  powers  of  less  than  a  microwatt.  At  higher  incident 
powers,  the  diode  operates  in  the  linear  region.  This  is  the  region  of  greatest  sensitivity. 
Thus,  biasing  the  diode  (via  the  Reference  Arm)  keeps  the  detector  operating  in  the  linear 
region.  Second,  there  is  a  phase  shifter  in  the  reference  arm  that  adjusts  the  phase  of  the 
propagating  microwaves.  The  phase  of  the  microwaves  in  the  reference  arm  are  adjusted 
to  match  the  phase  of  the  microwaves  reflected  from  the  cavity  containing  the  sample. 
With  proper  bias,  sample  power,  and  phase,  the  EPR  microwave  bridge  and  sample 
cavity  will  yield  optimal  EPR  spectra. 

2.4.  EPR  Sensitivity 

An  EPR  spectrometer  is  able  to  detect  very  small  concentrations  of  unpaired  spins. 
Under  optimal  conditions  it  is  easy  to  observe  an  EPR  signal  representing  1011  spins  [17]. 
Due  to  the  size  of  the  microwave  cavity,  samples  are  typically  no  more  than  0.03  cm'  in 
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volume.  Thus  the  10  spins  corresponds  a  concentration  of  approximately  3x10  cm  . 
Temperature  is  the  most  important  experimental  parameter  in  an  EPR  experiment 
(usually,  lower  temperature  means  larger  EPR  signals).  Consider  a  collection  of 
independent  electrons  with  an  upper  energy  state  associated  with  the  spin  +1/2  electrons 
and  a  lower  energy  state  associated  with  the  spin  -1/2  electrons.  The  populations  of  the 

-A  E 

two  states  follows  the  classical  Boltzmann  distribution  of  Nupper  —  Niowere  ksT,  where 
A E  is  the  energy  difference  between  the  two  spin  levels,  ks  is  the  Boltzmann  constant, 
and  T  is  the  temperature  in  Kelvin  units.  Greater  population  differences  for  the  spin 
states  (i.e.,  Niower  -  NUpper)  will  occur  at  lower  temperatures.  For  this  reason,  EPR  is 
conducted  at  lower  temperatures  because  the  intensity  of  the  EPR  signal  is  proportional 
to  the  population  difference  between  the  spin  levels.  Many  experiments  are  conducted  at 
temperatures  between  20  and  100  K.  At  extremely  low  temperatures,  the  spin-lattice 
relaxation  time  becomes  longer  and  this  may  result  in  microwave  power  saturation  of  the 
EPR  signal.  Avoiding  saturation  is  important  in  an  EPR  experiment,  as  it  will  limit  the 
detection  sensitivity.  Each  paramagnetic  defect  may  have  a  different  temperature 
dependence  of  its  spin-lattice  relaxation  time.  Therefore,  each  distinct  EPR  signal  will 
have  an  optimal  temperature  that  maximizes  its  intensity.  Each  optimal  temperature  must 
be  determined  experimentally. 

2.5.  Electron-Nuclear  Double  Resonance  (ENDOR) 

Another  experiment  that  can  be  performed  with  an  EPR  spectrometer  is  electron- 
nuclear  double  resonance,  referred  to  as  ENDOR.  ENDOR  uses  the  sensitivity  of  the 
EPR  spectrometer  to  perform  nuclear  magnetic  resonance  (NMR)  spectroscopy.  Many 
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nuclei  have  a  non-zero  magnetic  moment,  i.e.,  they  have  a  nuclear  spin.  In  general,  their 
magnetic  moments  are  three  orders  of  magnitude  smaller  than  those  of  electrons.  While  a 
frequency  in  the  gigahertz  range  (microwaves)  induces  EPR  transitions,  a  frequency  in 
the  megahertz  range  (radio  waves)  induces  nuclear  spin  transitions.  Nuclear  transitions 
are  described  by  the  same  principles  previously  discussed  for  electrons,  except  the 
selection  rule  is  now  Ams  =  0  and  Ami=  ±1.  The  energy  splitting  of  the  nucleus  due  to  its 
intrinsic  spin  is  called  the  nuclear  Zeeman  effect.  In  performing  ENDOR  with  an  EPR 
spectrometer,  one  “sits”  at  a  specific  magnetic  field  location  where  an  EPR  transition 
occurs.  Microwaves  are  provided  to  the  sample  such  that  the  EPR  transition  becomes 
saturated,  i.e,  the  populations  of  the  upper  and  lower  spin  states  become  nearly  equal. 
While  saturating  the  EPR  transition,  the  magnetic  field  of  the  radio  waves  and  the 
microwave  magnetic  field  must  both  be  perpendicular  to  the  static  magnetic  field.  The 
frequency  of  the  radio  waves  is  swept  across  a  region  of  interest.  Nuclear  transitions  will 
occur  at  specific  frequencies  of  the  radio  waves  when  the  difference  between  two  nuclear 
Zeeman  levels  equals  the  energy  of  the  radio  wave  photons.  During  these  nuclear 
transitions,  the  saturated  EPR  signal  will  change  in  magnitude  due  to  the  slight  variations 
in  populations  caused  by  the  nuclear  magnetic  resonances.  The  changes  in  the  intensity 
of  the  EPR  signal  as  a  function  of  the  radio  wave  frequency  is  the  ENDOR  spectrum. 

Like  EPR,  ENDOR  spectra  also  appear  as  the  first  derivative  of  an  absorption  spectrum. 
The  first  derivative  shape  results  from  modulating  the  frequency  of  the  radio  waves, 
usually  with  a  modulation  frequency  of  20  kHz  or  less. 

To  illustrate  ENDOR,  consider  the  example  of  an  electron  with  spin  S  =  1/2 
interacting  with  a  nucleus  of  spin  I  =  1/2.  The  energy  of  this  spin  system  was  described 
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by  Eq.  2.6.  Because  of  the  different  selection  rules  for  ENDOR,  a  nuclear  Zeeman 
interaction  term  must  be  included  in  the  spin  Hamiltonian.  The  resulting  energy  equation 
becomes 

£(ms,  rrq)  =  g  /?  ms  Bz  +  h  A  ms  nq  -  gNgNmlBz  (2.7) 

which  can  be  rewritten  as 

E(ms,  nq)  —  g  /?  ms  Bz  +  hA  irq  ms  —  hvN  nq  (where  vN  =  gN  PN  B/h)  (2.8) 
The  four  possible  energy  levels  are  now,  using  E(ms,  mi) 


/  1  1\  1  1  1 
E{+2'+2 )  "  28PBz  +  4hA~2l"’K 


Using  the  ENDOR  selection  rule  that  Ams  =  0  and  Anq=  ±1,  the  following  transitions 
occur. 
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The  energies  have  been  divided  by  h.  The  hyperfine  parameter  A  and  vN  have  units  of 
frequency,  typically  MHz.  Depending  on  the  values  of  A/2  and  vN,  there  are  two 
possible  forms  that  the  ENDOR  spectra  may  take. 
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Consider  first  the  case  where  A/2  >  vN.  ENDOR  transitions  occur  at  the  following 
frequencies  given  by  Eq  2.9. 

A  A 

Vr=--vN  andv2=-  +  vN  (2.9) 

This  gives  a  pair  of  ENDOR  lines  centered  on  A/2  and  separated  by  2  vN.  By  measuring 
the  separation  of  the  two  ENDOR  lines,  vN  can  be  directly  determined.  The  two  ENDOR 
lines  will  be  centered  on  A/2. 

The  second  case,  for  vN  >  A/2,  gives  ENDOR  transitions  at 

A  A 

Vi  =  VN  -  -  and  v2  =  vN  +  -  (2.10) 

This  results  in  a  pair  of  lines  centered  on  vN  and  separated  by  A.  These  two  possible 
cases  are  shown  in  Fig.  2.5.  ENDOR  can  measure  hyperfine  parameters  that  are  not 
resolved  in  EPR  experiments.  Each  element  in  the  periodic  table,  with  a  non-zero  nuclear 
spin,  has  its  own  unique  value  of  vN .  Therefore,  the  ENDOR  technique  can  be  used  to 
determine  the  identity  of  the  nucleus  responsible  for  a  hyperfine  interaction. 

A  major  drawback  to  ENDOR,  however,  is  the  fact  that  not  all  EPR  signals  can  be  easily 
microwave  power  saturated.  When  attempting  to  obtain  an  ENDOR  spectrum,  the 
temperature  and  microwave  power  must  be  optimized  if  success  is  to  be  achieved.  In 
some  cases,  ENDOR  spectra  are  easily  observed,  while,  in  other  cases  no  amount  of 
searching  will  reveal  an  ENDOR  spectrum.  The  intensity  of  an  ENDOR  spectrum  is 
usually  one  to  two  percent  of  the  intensity  of  the  corresponding  EPR  spectrum. 

2.6.  Spin  Hamiltonian 

An  isolated  electron  would  produce  a  single  peak  in  an  EPR  spectrum.  In  crystals, 
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Figure  2.5.  The  two  possible  cases  for  an  ENDOR  spectrum  where  S  =  1/2  and  I  =  1/2. 


however,  unpaired  electrons  are  influenced  by  their  environment  which  leads  to  their 
having  a  host  of  possible  energies,  and  thus  more  complex  EPR  spectra.  Some  of  these 
influences  on  spin  have  already  been  discussed.  It  is  difficult  to  establish  a  Hamiltonian 
that  describes  all  of  the  influences  on  an  electron,  then  solve  Schrodinger’s  equation 
using  this  Hamiltonian,  and  finally  relate  the  resulting  wave  function  to  the  experimental 
EPR  data.  Therefore,  simplifications  are  introduced  which  lead  to  the  concept  of  a  spin 
Hamiltonian.  The  spin  Hamiltonian  includes  only  the  spin-dependent  terms  in  the 
Hamiltonian.  It  is  used  to  derive  the  discrete  spin  energy  levels  and  determine  the 
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possible  transitions  between  these  levels.  The  energy  levels  are  the  eigenvalues  of  the 
spin  Hamiltonian  matrix  representing  the  total  energy  of  the  spin  system. 

Point  defects,  which  are  typical  of  unpaired  spin  systems  in  crystals,  often  exhibit 
anisotropic  behavior.  Defect  anisotropy  requires  that  the  spin  Hamiltonian  be  written  in 
matrix  form  using  the  appropriate  set  of  spin  vectors  as  a  basis  set.  The  spin  Hamiltonian 
must  be  expressed  in  a  specific  coordinate  system.  Usually,  a  coordinate  system  is 
chosen  with  z  being  the  direction  of  the  static  magnetic  field  (referred  to  as  the  magnetic- 
field  coordinate  system).  The  g  and  A  matrices,  containing  the  spin-Hamiltonian 
parameters,  are  diagonal  in  their  principal-axes  coordinate  systems.  When  constructing 
the  spin-Hamiltonian  matrix,  the  g  and  A  matrices  must  be  rotated  into  the  magnetic-field 
coordinate  system.  For  most  spin  systems,  a  computer  must  be  used  to  diagonalize  the 
spin-Hamiltonian  matrix. 

Each  term  appearing  in  a  general  spin  Hamiltonian  will  now  be  discussed.  There  are 
several  excellent  books  available  that  provide  in  depth  discussions  of  a  spin  Hamiltonian 
[17,  24,25].  Thus  far  in  this  dissertation,  the  electron  Zeeman  term  containing  the  g 
matrix,  the  hyperfine  term  containing  the  A  matrix,  and  the  isotropic  nuclear  Zeeman 
term  have  been  introduced.  The  equations  describing  the  energies  associated  with  these 
terms  (Eqs.  2. 5-2. 8)  could  be  written  as  spin  Hamiltonians,  by  replacing  the  E  with  an  H 
and  including  vector  notation. 

The  goal  of  an  EPR  experiment  performed  on  a  crystal  is  to  characterize  a  defect 
(with  one  or  more  unpaired  spins)  by  proposing  a  spin  Hamiltonian  and  determining  its 
parameters.  One  begins  by  noting  the  number  of  lines  in  an  EPR  spectrum  and  the 
relative  spacing  of  these  lines.  This  suggests  the  appropriate  S  value  and  the  I  values 
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representing  resolved  hyperfine  interactions,  and  thus  determines  the  general  form  of  the 
spin  Hamiltonian  for  that  defect.  Next,  an  angular  dependence  study  of  the  spectrum  is 
conducted  and  the  results  are  used  as  input  data  for  a  fitting  program  that  determines  the 
final  values  of  the  spin  Hamiltonian  parameters.  After  the  spin  Hamiltonian  parameters 
are  determined,  a  model  for  the  defect  is  established  that  is  consistent  with  the  crystal 
structure  and  these  parameters.  This  will  lead  to  a  better  understanding  of  the  defect 
chemistry  of  the  material. 

The  complete  spin  Hamiltonian  for  an  EPR  experiment  can  be  written  as 

H  =  Hez  +  Hhf  +  HFS  +  Hnz  +  Hq  (2.11) 

where 


— '  — * 

Hez  =  HES  ■  g  ■  B  electron  Zeeman  interaction  (2.12) 

— >  < - >  — > 

Hhf  —  I-A-S  hyperfine  (nuclear-electron  spin)  interaction  (2.13) 

HFS  —  S  ■  D  ■  S  fine  structure  (electron-electron  spin)  interaction  (2.14) 
Hnz  =  gNfiNl  ■  B  nuclear  Zeeman  interaction  (2.15) 

Hq  —  l  ■  Q  ■  1  nuclear  electric  quadrupole  interaction  (2.16) 


In  the  above  expressions,  g".  A ,  and  D  are  three-dimensional  symmetric  matrices.  Each 
matrix  is  described  by  six  parameters,  i.e.,  three  principal  values  and  three  Euler  angles 
defining  the  principal  axes  directions.  Q  is  a  traceless  three-dimensional  symmetric 
matrix  and  is  described  by  five  parameters. 

HEz ,  the  electron  Zeeman  term,  describes  the  interaction  between  the  electron  spin 
and  the  applied  magnetic  field.  This  is  typically  the  dominant  (largest  value)  term  in  a 
spin  Hamiltonian.  The  electron  Zeeman  Hamiltonian  term  can  be  rewritten  as: 


Hez=  gBS-g-B=pB-(L  +  geS)  +  AL-S  (2.17) 
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There  are  two  terms  on  the  right  side  of  Eq.  2.17.  The  first  term  represents  the  energy 
associated  with  the  separate  coupling  of  the  orbital  angular  momentum  L  and  the  spin 
angular  momentum  S  with  the  magnetic  field,  while  the  second  term  represents  the  spin- 
orbit  energy  associated  with  the  coupling  of  L  and  S.  The  constant  A  is  the  spin  orbit 
coupling  constant.  With  some  mathematical  manipulation  of  Eq.  2.17,  the  g  matrix  can 
be  rewritten  as 

g=geT+2AA  (2.18) 

where  /  is  the  identity  matrix.  If  the  angular  momentum  of  a  system  is  due  solely  to  spin 
angular  momentum,  g’  should  be  isotropic  with  the  value  ge.  Any  anisotropy  in  g’  is  due 

to  the  presence  of  orbital  angular  momentum  which  couples  higher  excited  states  to  the 

< — » 

ground  state.  This  anisotropy  is  represented  by  the  A  matrix  [25]. 

Hhf  is  the  hyperfine  term.  It  represents  the  interaction  between  the  electron  spin  and 
the  nuclear  spin.  This  A  matrix  can  be  decomposed  into  isotropic  and  anisotropic  parts, 

< — »  /  < — »  < — k 

described  by  A  —  ( al  +  B)  where  a  is  the  scalar  Fermi  contact  term  and  represents  the 

isotropic  part  of  A.  The  Fermi  contact  term  is  the  magnetic  interaction  between  the 
electron  and  the  atomic  nucleus  when  the  electron  is  inside  that  nucleus.  It  is 
proportional  to  \ip2\  evaluated  at  the  nucleus  (only  s  orbitals  contribute  to  the  Fermi 
contact  term).  The  anisotropic  portion  of  the  A  matrix  corresponds  to  a  dipole-dipole 
interaction  between  the  electron  and  nuclear  magnetic  moments.  It  is  represented  by  the 
traceless  matrix  B  and  has  two  hyperfine  interaction  constants  (b  and  b').  In  its  principal 
axis  system,  B  is  written  as 
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f-b  +  b'  0  0  \ 

B  =  [  0  -b-b'  0 

\0  02 b) 

It  follows  that 


b  ^  Bzz  and  b  (Bxx  Byy ) 

where  b'  describes  the  deviation  from  axial  symmetry. 

The  fine  structure  term,  Hfs,  in  the  spin  Hamiltonian  accounts  for  interaction 
between  electrons  and  is  present  when  S  >  1/2.  Examples  where  this  term  is  important 
include  crystal  field  effects  on  transition-metal  ions,  spin  triplets,  and  biradicals.  The  D 
matrix  is  traceless  and,  in  its  principal  axis  system,  can  be  represented  by  an  axially 
symmetric  parameter  D  and  an  asymmetry  parameter  E.  It  can  be  written  as  HF s  = 
DxxSx2  +  DyySy2  +  DZZSZ2  and  expressed  in  matrix  form  as 
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This  dissertation  does  not  describe  an  experiment  where  a  fine- structure  matrix  was 
needed  in  the  spin  Hamiltonian. 

The  nuclear  Zeeman  term,  HNZ,  occurs  for  nuclei  with  1^0.  This  term  is  the  same  as 
the  electron  Zeeman  term,  except  that  now  it  is  the  nuclear  spin  interacting  with  a 
magnetic  field.  This  energy,  even  though  it  is  small,  needs  to  be  accounted  for  in  the  spin 
Hamiltonian.  It  is  especially  important  in  ENDOR  experiments. 

The  nuclear  electric  quadrupole  interaction,  Hq  is  the  interaction  of  a  nucleus  with  the 
electric  field  gradient  at  the  nucleus.  When  a  nucleus  with  I  >  1/2  interacts  with  an 
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electric  field  gradient,  the  nuclear  spin  levels  are  no  longer  equally  spaced  in  energy. 

This  unequal  spacing  of  the  energy  levels  is  clearly  seen  in  ENDOR  spectra.  The  nuclear 
electric  quadrupole  term  is  represented  by  the  coupling  matrix  Q  (which  is  traceless)  and 
contains  the  product  of  the  nuclear  electric  quadrupole  moment  and  the  gradient  of  the 
electric  field  at  the  nucleus  due  to  the  surrounding  electrons.  For  a  symmetric  charge 
distribution  such  as  in  a  cubic  crystal,  the  electric  field  gradient  is  zero  at  the  “regular” 
lattice  sites.  However,  if  either  the  overall  symmetry  of  the  lattice  or  the  local  symmetry 
of  the  point  defect  is  lower,  then  the  electric  field  gradient  may  be  nonzero.  The 
quadrupole  coupling  parameters  are  then  a  measure  of  the  relative  distribution  of  the 
surrounding  electric  charge  [24] .  The  nuclear  electric  quadrupole  term  can  have  a 
significant  impact  in  ENDOR  spectra.  The  effect  of  the  quadrupole  term  will  be 
demonstrated  in  Chapter  3  where  boron  nuclei  interact  with  an  S  =  1/2  trapped  hole 
center. 

2.7.  EPR  and  ENDOR  Instrumentation 

The  EPR  and  ENDOR  spectra  presented  in  this  dissertation  were  taken  with  an  X- 
band  EMX  Bruker  spectrometer  and  a  Bruker  ESP  300  spectrometer  (used  primarily  for 
ENDOR).  Bruker,  a  German  company,  is  the  primary  manufacturer  of  versatile  and 
sensitive  EPR  and  ENDOR  spectrometers  in  the  world.  Both  systems  are  EPR  units  with 
a  single  computer  control  interface.  For  low  temperatures,  an  Oxford  Instruments 
helium-gas  flow  system  was  used  to  maintain  the  samples  at  selected  temperatures  in  the 
4  -  200  K  range. 
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Three  types  of  resonant  cavities  were  used  with  the  EPR  spectrometers.  The  EMX 
spectrometer  uses  either  a  standard  rectangular  TE102  cavity  or  a  high  sensitivity  TE011 
cylindrical  cavity  while  the  Bruker  ESP  300  uses  a  cylindrical  TEon  ENDOR  cavity. 
The  sharpness  of  response  of  any  resonant  system  is  commonly  described  by  a  factor  of 
merit,  universally  represented  by  the  symbol  Q.  Q  is  defined  as 

_  27t  (maximum  microwave  energy  stored  in  the  resonator) 

^  energy  dissipated  per  microwave  cycle  (2-19) 


A  more  practical  “working”  definition  of  the  Q  factor  is  given  by 


Q  =  — 
^  Av 


(2.20) 


where  v  is  the  resonant  frequency  of  the  cavity  and  Av  is  the  full  width  at  half  maximum 
of  the  resonant  peak.  The  standard  rectangular  cavity  for  the  EPR  system  and  the 
ENDOR  cavity  both  have  Q  values  between  2000  and  5000.  The  high  sensitivity  cavity 
has  a  Q  value  greater  than  12,500.  A  higher  Q  value  corresponds  to  a  better  signal-to- 
noise  ratio.  Metal  and  water,  including  water  vapor,  are  some  examples  of  materials  that 
will  drastically  reduce  the  Q  of  the  cavity.  For  this  reason,  nitrogen  gas  continually  flows 
through  the  cavity  during  low  temperature  operation  to  prevent  moisture  buildup  on  the 
surface  of  the  glassware  within  the  cavity.  The  cylindrical  ENDOR  cavity  is  less 
sensitive  because  it  contains  the  metal  coil  that  produces  the  radio  waves  needed  to 
induce  the  ENDOR  transitions.  This  metal  coil  is  attached  to  the  cold  temperature 
glassware  within  the  cavity.  The  cavities  are  also  designed  such  that  the  microwave 
magnetic  field  is  a  maximum  and  the  microwave  electric  field  is  a  minimum  at  the 
sample  position. 
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2.8.  Optical  Absorption  Spectroscopy 

Optical  absorption  is  another  spectroscopic  technique  used  to  characterize  point 

he 

defects  in  crystals.  Photons  with  energy  E  =  hv  =  —  are  absorbed  by  the  point  defects, 

A 

i.e.,  transitions  occur  from  the  ground  state  to  one  or  more  excited  states.  The  electronic 
structure  of  the  defect  (and  the  crystal)  determines  the  possible  absorption  transitions. 
Because  the  electronic  energy  levels  in  matter  are  quantized  and  there  is  a  forbidden 
energy  band  in  semiconductors  and  insulators,  only  specific  photon  energies  are 
absorbed.  When  energy  is  absorbed  by  an  atomic  scale  defect,  it  is  excited  to  a  higher 
electronic  state.  The  defect  emits  its  absorbed  energy  and  returns  to  the  ground  state. 
There  are  two  types  of  emission  processes,  radiative  and  non-radiative.  Radiative 
emissions  occur  when  the  defect  emits  a  photon  and  nonradiative  occurs  when  the  defect 
emits  one  or  more  phonons  to  the  surrounding  lattice.  These  photons  and  phonons  have 
discrete  values  of  energy  that  are  controlled  by  the  electronic  structure  of  the  defect  and 
the  surrounding  lattice.  The  principles  of  absorption  and  emission  are  best  described  by 
the  Franck-Condon  principle  with  the  help  of  a  configuration-coordinate  energy  diagram, 
illustrated  in  Fig.  2.6. 

For  a  given  electronic  state  of  a  point  defect,  the  surrounding  ions  have  an 
equilibrium  configuration,  i.e.,  this  is  the  positions  of  the  ions  that  corresponds  to  a 
minimum  energy  for  the  defect.  There  will  be  one  equilibrium  configuration  when  the 
defect  is  in  its  electronic  ground  state  and  a  different  equilibrium  configuration  when  it  is 
in  an  excited  electronic  state.  For  this  reason,  the  horizontal  axis  of  Fig.  2.6  is  labeled 
“nuclear  configuration”  and  is  often  referred  to  as  the  configuration  coordinate  Q.  The 
parabolic  curves  represent  the  potential  wells  of  the  electronic  ground  state  (So)  and  the 
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Figure  2.6.  Configuration-coordinate  energy  diagram. 


first  excited  state  (Si)  while  the  horizontal  lines  in  the  wells  represent  the  discrete 
vibrational  modes  of  the  structure.  The  Franck-Condon  principle  states  that  the  change  in 
the  electronic  structure  of  the  defect  occurs  very  quickly,  before  the  surrounding  ions 
have  a  chance  to  adjust  to  their  new  equilibrium  positions.  This  is  shown  by  the  vertical 
arrow  in  Fig.  2.6  labeled  Absorption  (0,1).  The  (0,1)  label  on  the  vertical  arrow  indicates 
that  the  electron  goes  from  the  lowest  vibrational  mode  (n  =  0)  in  the  ground  state  to  the 
second  vibrational  mode  (n  =  1)  in  the  excited  state.  After  this  upward  electronic 
transition  (absorption),  the  defect  will  release  its  added  energy  in  three  steps.  First,  the 
surrounding  ions  shift  to  their  new  equilibrium  configuration  and  lower  the  energy  of  the 
defect  by  emitting  a  phonon  (shown  with  a  red  arrow  in  Fig.  2.6).  Then,  following  the 
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Franck-Condon  principle,  the  defect  emits  a  photon  (fluorescence)  as  it  returns  to  the 
ground  state.  As  it  reaches  the  ground  state,  the  surrounding  ionic  configuration  still 
corresponds  to  the  excited  state  equilibrium  configuration.  Finally,  the  ionic 
configuration  returns  to  the  original  ground  state  equilibrium  configuration  by  emitting  a 
phonon  (shown  with  a  blue  arrow  in  Fig.  2.6). 

Absorption  measurements  provide  useful  information  about  the  electronic 
structure  of  the  defect  and  often  can  be  correlated  with  EPR  measurements,  as  will  be 
seen  later  in  this  dissertation  for  LTB  samples  doped  with  silver.  A  Varian-Cary 
spectrophotometer  was  used  to  measure  room  temperature  optical  absorption  in  both 
undoped  and  doped  LTB  crystals.  Absorption  measurements  were  also  conducted  on 
pre-  and  post-x-ray  irradiated  LTB  crystals.  The  optical  absorption  data  presented  in  this 
dissertation  was  obtained  by  passing  electromagnetic  radiation,  with  wavelengths  from 
the  near  ultraviolet  to  the  near  infrared  region,  through  the  crystal  and  thus  exciting 
transitions  between  electronic  states. 

A  diagram  that  is  similar  to  the  spectrophotometer  used  for  the  present  study  is 
shown  in  Lig.  2.7.  The  spectrometer  contains  two  light  sources,  one  for  wavelengths 
between  190  and  350  nm  and  one  for  wavelengths  between  350  nm  and  3.2  pm.  Light 
passes  through  a  monochromator  that  selects  the  wavelength  range  for  which  absorption 
is  to  be  measured.  The  monochromator  contains  a  diffraction  grating  that  disperses  the 
light  into  individual  wavelengths.  The  grating  allows  only  a  narrow  range  of 
wavelengths  to  pass  through  a  slit  and  into  the  sample  chamber.  A  chopper  causes  the 
beam  to  alternately  pass  through  the  sample  and  reference  areas.  Both  beams  are  then 
incident  on  the  detector.  Two  detectors  are  used,  a  photomultiplier  tube  for  wavelengths 
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from  190  to  850  nm  and  a  semiconducting  PbS  cell  for  wavelengths  from  850  nm  to  3.2 
|um.  The  detector  is  used  to  compare  the  intensities  of  the  alternating  beams  and 
produces  an  output  signal  proportional  to  their  difference.  This  signal,  obtained  as  a 
function  of  wavelength,  represents  the  absorbed  light.  The  output  can  be  displayed  as 
either  a  transmission  curve  or  an  absorption  curve. 


Light  Source  Narrow  slit 


Figure  2.7.  Diagram  of  a  spectrophotometer. 


Optical  absorption  spectra  are  often  plotted  as  optical  density  (or  absorbance) 
versus  wavelength.  Optical  density  (OD)  is  defined  as 

OD  =  log10  0  (2.21) 
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where  Io  is  the  incident  light  intensity  and  I  is  the  transmitted  light  intensity.  The 
Lambert-Beer  law  (Eq.  2.22)  relates  I  and  Io  to  the  absorption  coefficient  a  and  the 
optical  path  length  x  in  the  sample. 


/  =  IQe~ax  (2.22) 

Combining  Eqs.  2.21  and  2.22  gives  the  following  expression  for  the  absorption 
coefficient,  shown  in  Eq.  2.23. 


a  = 


OD 

xlog(e) 


(2.23) 


Absorption,  fluorescence  spectroscopy,  and  EPR  are  combined  to  show  that  Ag2+  ions  are 
responsible  for  an  optical  absorption  peak  found  in  x-ray  irradiated  LTB  crystals. 


2.9.  Fluorescence  Spectroscopy 

Fluorescence  spectroscopy  is  a  similar  technique  to  that  of  absorption  described  in 
Section  2.8  with  the  exception  that  a  fluorescent  spectrometer,  or  fluorospectrometer,  can 
measure  and  monitor  specific  absorptions  and  emissions.  This  is  accomplished  using  two 
different  optical  monochromators,  one  specifically  for  each  function  of  emission  or 
excitation.  The  instrument,  illustrated  in  Fig.  2.8,  produces  output  that  is  given  as 
intensity  (photomultiplier  tube  counts)  versus  wavelength.  The  spectrometer  used  in  this 
dissertation  is  a  Fluorolog  -3  fluorospectrometer  manufactured  by  the  Horiba  Scientific 
Corporation. 

The  Fluorolog-3  uses  a  xenon  high  voltage  lamp  as  the  light  source.  The  light 
then  passes  through  a  narrow  slit  into  the  excitation  monochromator.  In  the  excitation 
monochromator,  the  light  is  separated  into  individual  wavelengths  using  two  sets  of 
mirrors  and  two  diffraction  gratings  to  provide  a  narrow  band  of  wavelengths  to  the 
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sample.  There  are  two  slits  in  the  excitation  monochromator  that  can  be  used  to  adjust 
the  resolution  of  the  exciting  light,  from  1  to  14  nm.  After  passing  through  the  final  slit 
of  the  excitation  monochromator,  light  is  incident  on  the  sample.  The  sample  holder  used 
in  this  dissertation 


Sample 


Figure  2.8.  Diagram  of  a  fluorescence  spectrometer  (also  referred  to  as  a  fluorolog). 


was  attached  to  an  movable  stage  such  that  the  beam  could  be  focused  to  a  small  spot  size 
on  the  sample.  Once  light  reaches  the  sample,  a  portion  will  be  reflected  and  a  portion 
will  be  absorbed  and  then  emitted  as  fluorescent  light.  The  emitted  light,  and  any 
scattered  light,  passes  through  the  emission  monochromator  and  is  collected  at  the 
detector,  i.e.,  a  photomultiplier  tube  in  the  case  of  the  Fluorolog-3  instrument.  The 
emission  monochromator  also  has  two  diffraction  gratings  and  two  controlled  slits. 
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A  computer  is  used  to  control  the  wavelength  settings  for  both  monochromators 
as  well  as  the  slit  resolutions.  The  Fluorolog  can  be  set  up  to  run  either  an  excitation  or 
emission  experiment.  In  excitation,  the  emission  monochromator  can  be  set  to  monitor  a 
specific  wavelength.  Conversely,  in  emission,  the  excitation  monochromator  is  set  to 
provide  a  specific  wavelength.  Normally,  a  series  of  intensity  vs  emission  spectra  are 
taken  at  different  excitation  wavelengths  in  order  to  find  the  optimum  excitation 
wavelength  for  a  particular  emission  peak.  Excitation  spectra  tend  to  be  slightly  narrower 
bands  while  emission  curves  tend  to  be  broader  due  to  phonon  interaction.  The 
broadening  of  the  emission  and  excitation  bands  can  be  described  by  the  Huang-Rhys 
factor  [26],  but  is  not  discussed  in  this  dissertation. 

While  an  absorption  spectrometer  provides  data  that  does  not  need  correction, 
there  are  several  corrections  that  must  be  made  in  data  taken  with  a  fluorospectrometer. 
First,  the  intensity  of  the  light  provided  by  the  xenon  lamp  varies  with  wavelength.  The 
data  collected  at  the  photomultiplier  tube  must  be  adjusted  for  a  constant  photon  flux. 
Second,  the  photomultiplier  tube  has  a  nonconstant  response  to  wavelength  as  well.  Both 
of  these  adjustments  may  cause  shifts  in  peak  positions  and  intensities  of  the  acquired 
spectra. 
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III.  Characterization  of  the  Intrinsic  Hole  Center  in  Lithium 


Tetraborate 

3.1.  Introduction 

Point  defects  have  been  widely  studied  in  various  oxide  materials  [27-30]. 
Depending  on  their  effective  charge  relative  to  the  lattice,  these  defects  may  either  trap 
electrons  or  holes  when  a  crystal  is  exposed  to  ionizing  radiation.  Examples  of  common 
point  defects  in  oxides  are  anion  vacancies,  cation  vacancies,  and  substitutional 
transition-metal  impurities.  Optical  spectroscopy,  electron  paramagnetic  resonance 
(EPR),  and  electron-nuclear  double  resonance  (ENDOR)  are  often  used  to  study  these 
point  defects. 

Previous  research  at  the  Air  Force  Institute  of  Technology  investigated  intrinsic 
point  defects  in  nominally  undoped  single  crystals  of  lithium  tetraborate  (LTB)  [31]. 

EPR  was  used  to  identify  oxygen  and  lithium  vacancies  as  the  primary  point  defects  in 
the  as-grown  LTB  crystals.  To  be  observed  in  EPR  experiments,  the  defects  need  to  be  in 
a  paramagnetic  charge  state.  In  the  case  of  LTB,  this  is  accomplished  by  placing  the 
crystal  in  liquid  nitrogen,  thus  cooling  the  crystal  to  77  K,  and  then  while  in  liquid 
nitrogen,  irradiating  the  crystal  with  x  rays  (60  kV  and  30  mA)  produced  by  a  high-flux 
x-ray  tube.  Energy  provided  by  the  x  rays  moves  electrons  from  the  valence  band  to  the 
conduction  band  and  results  in  the  creation  of  electron-hole  pairs.  Many  of  these 
electrons  and  holes  recombine  immediately,  while  others  move  separately  through  the 
lattice  and  become  “trapped”  at  defects.  Defects  that  trap  electrons  are  referred  to  as 
electron  centers  and  defects  that  trap  holes  are  referred  to  as  hole  centers.  Trapped 
electron  centers  have  gained  electrons  from  the  lattice  and  trapped  hole  centers  have  thus 
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lost  electrons.  The  concentration  of  these  two  types  of  defects  must  balance  to  maintain 
charge  neutrality.  These  electron  and  hole  centers  form  at  different  locations  in  the 
crystal.  In  LTB,  the  electrons  are  trapped  at  oxygen  vacancies  and  the  holes  are  trapped 
near  lithium  vacancies.  Following  an  irradiation  at  77  K,  these  trapped  electron  and  hole 
centers  are  stable  as  long  as  the  sample  remains  below  ~90  K.  Above  ~90  K,  the 
electrons  and  holes  recombine  and  the  crystal  returns  to  its  non-paramagnetic  charge 
state.  The  present  chapter  describes  an  in-depth  EPR  and  ENDOR  study  of  the  intrinsic 
trapped  hole  center  in  undoped  LTB  crystals. 

3.2.  Thermoluminescence  Results 

A  thermoluminescence  (TL)  experiment  was  performed  to  establish  the  thermal 
stability  of  the  intrinsic  trapped  hole  center  in  LTB .  Though  there  is  considerable 
literature  describing  thermoluminescence  from  lithium  tetraborate  crystals  containing 
various  dopants  [32-34],  thermoluminescence  from  undoped  lithium  tetraborate  at  low 
temperatures  has  not  been  reported.  Figure  3.1  is  a  TL  glow  curve  obtained  from  an 
undoped  LTB  crystal.  These  data  were  acquired  from  a  77  K  irradiated  crystal  that  was 
quickly  removed  from  liquid  nitrogen  and  placed  at  the  end  of  an  optical  fiber.  As  the 
crystal  rapidly  warmed  to  room  temperature,  light  was  collected  through  the  optical  fiber 
and  passed  into  an  optical  spectrometer.  Temperatures  were  determined  by  assuming  a 
linear  heating  rate.  The  sample  was  at  77  K  immediately  upon  placement  at  the  end  of 
the  fiber,  then  a  timer  was  started,  and  the  elapsed  time  was  noted  when  the  crystal 
reached  273  K  (the  moment  when  frost  converted  to  moisture  on  the  crystal  surface). 

The  horizontal  axis  of  Fig.  3.1  reflects  this  time-to-temperature  conversion.  The  strong 
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peak  just  above  100  K  agrees  with  the  90  K  EPR  results  reported  by  Swinney  et  al.  [31]. 
The  heating  rate  used  in  the  TL  experiment  was  much  faster  (13  K/s)  than  the  heating  rate 
(<  1  K/s)  used  in  the  EPR  experiment,  thus  explaining  why  the  TL  peak  occurs  at  a 
slightly  higher  temperature.  The  glow  from  the  crystal  was  intense  enough  to  be 
observed  with  the  naked  eye.  The  emission  at  100  K  in  the  TL  curve  is  assigned  to 
electron-hole  recombination,  where  the  participating  hole  is  the  subject  of  the  present 
chapter. 


Fig.  3.1.  TL  glow  curve  of  an  LTB  crystal  irradiated  with  x  rays  at  77  K. 
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3.3.  Electron  Paramagnetic  Resonance  Results 

In  the  earlier  study  by  Swinney  et  al.  [31],  the  intrinsic  paramagnetic  electron  and 
hole  centers  were  shown  to  have  g  values  near  2.0,  hence  the  EPR  spectra  of  the  x-ray 
induced  defects  overlap.  The  g  shifts  associated  with  these  two  defects  are  relatively 
small  and  all  directions  of  magnetic  field  have  overlapping  electron  and  hole  center 
spectra.  With  the  magnetic  field  parallel  to  the  c  axis  of  the  crystal,  or  [001]  direction, 
the  electron  center  has  a  spectrum  that  extends  over  400  Gauss  and  is  centered  at  g  = 
2.0049  while  the  hole  center  extends  over  45  Gauss  and  is  centered  at  g  =  2.0 1 1 6.  By 
acquiring  the  hole  center  spectrum  at  40  K,  instead  of  the  typical  25  K  used  for  the 
electron  center,  microwave  power  saturation  effects  are  reduced  and  the  well-resolved 
structures  that  are  associated  with  the  hole  center  can  be  clearly  seen.  In  practice,  the 
hole  center  spectrum  obscures  a  portion  of  the  electron  center  spectrum.  Figure  3.2 
shows  the  EPR  spectrum  of  the  intrinsic  hole  center  in  LTB.  The  three  traces  in  Fig.  3.2 
were  taken  with  the  magnetic  field  parallel  to  the  [001],  [100],  and  [110]  directions, 
respectively,  in  the  crystal. 

The  seven  line  spectrum  observed  in  Fig.  3.2(a)  is  representative  of  an  S  =  1/2  spin 
interacting  with  two  1  'B  nuclei  where  each  has  a  nuclear  spin  of  I  =  3/2.  The  energy  level 
diagram  in  Fig.  3.3  explains  the  EPR  spectrum  (shown  in  Fig.  3.2(a))  taken  when  the  c 
axis  of  the  crystal  is  aligned  parallel  to  the  magnetic  field.  The  hole  has  a  spin  of  S  =  1/2. 
In  a  magnetic  field,  this  spin  degeneracy  is  lifted  and  creates  ms  =  +1/2  and  ms  =  -1/2 
energy  levels.  One  boron  atom  with  a  nuclear  spin  of  I  =  3/2  causes  each  of  these  levels 
to  split  into  four.  These  levels  are  labeled  mi(l)  in  Fig.  3.2.  Next,  each  mi(l)  level  is 
split  into  four  by  a  second  boron  atom,  again  with  I  =  3/2.  This  results  in  16  levels 
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associated  with  ms  =  +1/2  and  16  levels  associated  with  ms=  -1/2.  These  32  levels  are 


labeled  as  mi(2).  The  16  total  energy  levels  possible  for  each  ms  are  displayed  in  Fig.  3.3 
in  the  section  describing  the  second  boron  interaction.  As  an  aid  to  understanding  Fig. 
3.3,  the  reader  can  follow  across  the  ms  =  +1/2  portion  of  the  diagram.  On  the  left 
side,  find  mi(l)  =3/2  underlined  in  blue  and  then  look  at  the  four  blue  lines  just  to  the 
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Figure  3.2.  EPR  spectrum  of  the  intrinsic  trapped  hole  center  in  an  IJ2B4O7  crystal. 
Stick  diagrams  show  the  1  'B  hyperfine  lines,  (a)  Magnetic  field  along  the  [001]  direction, 
(b)  Magnetic  field  along  the  [100]  direction,  (c)  Magnetic  field  along  the  [110]  direction. 
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Energy  level  diagram  for  an  S=l/2  interacting  with  2 
equivalent  “Boron  nuclei  of  1=3/2 


First  Boron  Second  Boron  interacting  and  subsequent  splitting 
interacting  w/  of  first  Boron  energy  levels  (lines  in  first  boron 
S=l/2  coordinate  with  lines  in  each  interaction) 


Final  Energy  levels 
(length  corresponds  to  #  of 
degenerate  levels) 


EPR  Transitions 


Transitions  (color  denotes  equivalence) 


Figure  3.3.  Energy  level  diagram  for  an  electron  spin  S  =  1/2  interacting  with  two 

equivalent  1  'B  nuclei  of  spin  I  =  3/2. 


right  in  the  “Second  boron  interacting. . .”  portion.  These  four  levels  represent  the  only 
energies  that  are  possible  when  ms  =  +1/2  and  mi(l)  =  +3/2.  These  levels  are  labeled 
using  the  scheme  E(ms,  mi(l),  mi(2)). 

E(+l/2,  +3/2,  +3/2) 

E(+l/2,  +3/2,  +1/2) 

E(+l/2,  +3/2, -1/2) 

E(+l/2,  +3/2,  -3/2) 

Notice  that  many  of  the  levels  in  Fig.  3.3  are  degenerate  with  each  other.  For  example, 
E(+l/2,  +3/2,  -1/2),  E(+l/2,  +1/2,  +1/2),  and  E(+l/2,  -1/2,  +3/2)  have  the  same  energy, 
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while  E(+l/2,  +3/2,  +3/2)  is  not  degenerate  with  any  other  levels.  Following  across  to 
the  right  side  of  Fig.  3.3,  the  final  energy  levels  are  drawn  using  line  lengths  that 
correspond  to  the  degree  of  degeneracy  for  each  level.  Fooking  down  on  the  “Final 
Energy  Fevels”  portion  corresponding  to  ms  =  +1/2,  the  bar  length  pattern  of  1,  2,  3,  4,  3, 
2,  1  in  terms  of  degeneracy  can  be  seen.  This  is  the  vertical  intensity  pattern  observed  in 
Fig.  3.2(a).  Recalling  that  the  selection  rules  for  EPR  transitions  are  Ams  =  ±1  and  Ami  = 
0,  the  dashed  arrows  on  the  right  side  of  Fig.  3.3  correspond  to  the  possible  EPR 
transitions,  of  which  there  are  seven.  The  small  arrows  between  the  final  energy  levels 
correspond  to  ENDOR  transitions  where  Ams  =  0  and  Ami  =  ±1. 

Although  it  is  not  shown,  the  final  energy  levels  on  the  right  side  of  Fig.  3.3  are 
not  equally  separated  in  energy.  The  1  'B  nucleus  has  a  nuclear  electric  quadrupole 
moment  and  this  causes  small  but  measurable  changes  in  the  separations  of  the  energy 
levels  within  the  ms  =  +1/2  and  ms  =  -1/2  sets.  The  arrows  in  this  figure  are  color  coded 
to  show  which  transitions  have  the  same  energy  differences.  Consider  the  ms  =  +1/2 
energy  level  portion  of  the  figure.  There  are  six  energy  levels,  but  there  are  only  three 
distinct  energy  differences  between  levels,  i.e.,  the  red,  blue,  and  green  arrows.  Because 
of  this  nuclear  electric  quadrupole  interaction,  an  ENDOR  spectrum  would  show  three 
transitions  associated  with  ms  =  +1/2.  The  sets  of  three  ENDOR  lines  are  shown  in  the 
ENDOR  spectra  discussed  later  in  this  chapter.  These  sets  of  three  ENDOR  lines  also 
emerge  from  the  analysis  of  the  complete  spin  Hamiltonian.  Equivalent  hyperfine  and 
quadrupole  matrices  for  the  two  boron  nuclei  interacting  with  the  unpaired  electron,  i.e., 
the  hole,  cause  the  sets  of  three  lines.  If  the  two  boron  nuclei  did  not  have  equivalent 
matrices,  there  would  be  six  different  energy  separations,  and  thus  six  discrete  ENDOR 
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transitions.  This  behavior  is  very  useful  in  determining  a  defect  model.  In  a  practical 
sense,  it  is  useful  when  aligning  the  magnetic  field  to  be  in  particular  crystal  planes  so 
that  useful  EPR  and  ENDOR  data  could  be  collected  and  interpreted. 

In  order  to  understand  the  spectra  presented  in  Fig.  3.2,  the  crystal  structure  of 
lithium  tetraborate  must  be  discussed.  The  basic  building  block  in  LTB  is  a  (B4O9)6 
unit.  Each  of  these  units  is  surrounded  by  6  Li+  ions,  thus  maintaining  charge  neutrality. 
These  (B4O9)6  units  are  stacked  and  rotated  by  90°  about  the  c  axis  of  the  crystal  to  give 
four  crystallographically  equivalent  orientations  of  the  basic  unit.  The  units  do  not  sit 
directly  above  and  below  each  other  because  they  rotate  along  a  screw  axis.  A  screw  axis 
is  a  rotation  followed  by  a  translation  along  the  direction  of  the  axis.  The  four  (B4C>9)6_ 
units  can  be  seen  in  Fig.  3.4(a)  where  the  lithium  atoms  have  been  removed  for  better 
viewing.  Notice  that  the  four  (B4O9)6-  units  are  identical  other  than  rotation.  The  change 
in  size  demonstrates  the  relative  depths  of  the  (B4C>9)6~  units.  All  other  (B409)6-  units  in 
the  crystal  will  be  crystallographically  equivalent  to  one  of  these  four  or  will  be  mirror 
images  of  these  four  (this  makes  a  total  of  eight  distinguishable  units).  A  defect  can  be 
located  in  any  of  these  crystallographically  equivalent  units.  A  defect,  therefore,  has 
eight  possible  orientations  of  the  principal  axes  directions  associated  with  its  g  matrix  and 
hyperfinc  matrices.  These  eight  orientations,  or  “sites,”  can  often  be  resolved  in  an  EPR 
experiment.  For  the  remainder  of  this  dissertation,  “sites”  will  be  used  to  refer  to 
crystallographically  equivalent  orientations  of  a  defect,  thus  there  are  eight  sites  in  LTB 
for  any  given  defect. 

When  the  magnetic  field  is  parallel  to  the  c  axis  of  an  LTB  crystal,  all  eight  of  the 
crystallographically  equivalent  sites  of  a  defect  are  magnetically  equivalent.  The  EPR 
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Figure  3.4.  Four  (B4O9)6"  units  comprising  the  LTB  crystal  lattice,  (a)  when  viewed 
along  the  c  axis  and  (b)  when  viewed  looking  along  the  a  axis. 

spectrum  in  Fig.  3.2(a)  represents  this  case  where  the  intrinsic  trapped  hole  center  has  the 
simplest  spectrum,  consisting  of  seven  lines,  because  all  of  the  sites  are  magnetically 
equivalent.  In  contrast,  when  the  magnetic  field  is  parallel  to  the  a  axis  of  the  crystal,  the 
eight  crystallographically  equivalent  sites  divide  into  two  distinct  groups  of  magnetically 
equivalent  sites.  Figure  3.4(b),  viewed  along  the  a  axis,  shows  four  (B4O9)6”  units.  Each 
of  these  four  also  has  a  mirror  to  make  a  total  of  eight.  The  upper  two  sites  in  Fig  3.4(b) 
are  magnetically  equivalent,  as  are  the  lower  two  sites.  The  equivalent  sites  are  mirror 
images  of  each  other.  When  EPR  is  performed  with  the  magnetic  field  parallel  to  the 
crystal  a  axis,  there  will  be  two  sets  of  EPR  lines  in  the  spectrum.  This  behavior  is 
illustrated  in  Fig.  3.2(b).  A  full  seven  peak  set  of  lines  similar  to  Fig.  3.2(a)  can  be  seen 
on  the  right  side  of  Fig.  3.2(b),  extending  from  3348  to  3384  G.  A  second  set  of  seven 
lines  underlay  the  larger  initial  set  of  seven  lines  and  has  two  components  that  extend  to 
lower  magnetic  field,  at  3336  and  3342  G.  This  is  shown  in  the  stick  diagram  above  the 
spectrum  in  Fig.  3.2(b).  Thus,  when  the  magnetic  field  is  along  the  a  axis,  the  two  sets  of 
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lines  representing  the  two  magnetically  inequivalent  sites  are  very  similar,  yet 
distinguishable,  in  the  EPR  spectrum. 

The  farthest  apart  that  the  sets  of  lines  ever  move  from  each  other  is  shown  in  Fig. 
3.2(c)  and  was  taken  with  the  magnetic  field  parallel  to  the  [110]  crystal  direction.  The 
fully  visible  spectrum  on  the  right  has  a  g  value  of  2.005  and  the  mostly  visible  spectrum 
on  the  left  has  a  g  value  very  close  to  2.020.  The  g  value  for  the  left  set  of  lines  is  an 
estimate  of  where  the  center  of  the  spectrum  occurs,  since  the  complete  seven-line 
spectrum  is  not  observed.  The  average  g  value  of  the  two  sets  of  lines  is  2.0125,  and  this 
is  also  the  value  of  g  along  the  c  axis  and  the  average  of  the  two  spectra  observed  when 
the  field  is  parallel  to  the  a  axis.  The  EPR  spectra  of  the  intrinsic  trapped  hole  center 
never  shift  far  from  this  2.0125  g  value.  This  lack  of  angular  dependence  means  that  the 
g  matrix  for  this  defect  is  nearly  isotropic. 

During  the  characterization  of  the  intrinsic  trapped  hole  center,  it  was  impossible 
to  observe  more  than  one  complete  seven  line  set  with  EPR.  This  meant  that  an  exact  g 
matrix  could  not  be  determined  from  the  experimental  data.  In  the  spin  Hamiltonian,  the 
g  matrix  is  described  by  three  principal  values,  gx,  gy,  gz,  and  three  Euler  angles  that 
specify  the  directions  of  the  principal  axes  relative  to  the  crystal  axes.  Each  observed 
seven-line  set  can  be  used  as  a  single  data  point  in  fitting  the  g  matrix.  Thus,  six 
independent  seven  line  sets  are  required  to  determine  the  complete  g  matrix.  Only  four 
independent  seven-line  sets  could  be  obtained  for  the  intrinsic  trapped  hole  center,  when 
the  magnetic  field  was  parallel  to  the  [001],  [100],  [110],  and  [101]  crystal  directions. 
Data  taken  with  the  field  parallel  to  the  [010]  crystal  direction  does  not  help  since  it  is 
equivalent  to  data  taken  when  the  field  is  parallel  to  the  [100]  crystal  direction  in  this 
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tetragonal  crystal.  Additional  independent  sets  of  lines  from  other  orientations  of  the 
magnetic  field  would  have  provided  the  extra  data  points  needed  for  a  complete 
determination  of  the  g  matrix. 

Although  an  experimental  g  matrix  could  not  be  precisely  determined,  it  is  still 
possible  to  identify  the  EPR  spectrum  as  “holelike”.  The  individual  g  values  that  were 
obtained  from  the  data  in  Fig.  3.2  ranged  from  2.0050  to  2.0125.  In  a  simple  ionic 
picture  (which  ignores  the  high  degree  of  covalency  within  the  borate  units  in  LEEEOy), 
the  O  ion  that  contains  the  hole  would  have  a  2p5  electron  configuration  (2px2,  2py2, 

2pz  )  with  L=l,  S=l/2.  The  threefold  orbital  degeneracy  of  this  “P  state  of  the  O  ion  is 
then  removed  by  the  crystalline  electric  field.  In  ascending  order,  the  three  energy  levels 
would  be  Ei,  E2,  and  E3.  Shifts  of  the  principal  values  from  the  free-spin  value  (ge  = 
2.0023)  are  caused  by  spin-orbit  interactions  that  admix  excited  states  to  the  ground  state. 
This  very  simple  analysis  (based  on  the  ionic  picture)  predicts  the  following  g  values. 

gl  =  ge  (3-1) 

21 

g2=ge~  ~ — ~  (3-2) 

g3=ge~  ^  (3-3) 

The  spin  orbit  coupling  constant,  X,  for  the  O  ion  is  approximately  -135  cm1.  Equations 
3.2  and  3.3  predict  positive  shifts  in  g  from  the  free-spin  value  for  the  O  ion.  In  general, 
a  hole  trapped  on  an  oxygen  ion  has  g  values  that  vary  from  2.003  to  2.050  [35-38]. 

These  theoretical  predictions  are  consistent  with  the  estimated  range  of  experimental  g 
values  (2.0050  to  2.0125),  thus  verifying  that  the  EPR  spectra  in  Fig.  3.2  are 
representative  of  a  trapped  hole. 
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3.4.  Electron-Nuclear  Double  Resonance  Results 


This  investigation  of  the  intrinsic  trapped  hole  center  in  LTB  will  proceed  with 
the  assumption  of  an  isotropic  g  matrix  (in  other  words,  the  small  differences  predicted  in 
Eqs.  3.1,  3.2,  and  3.3  will  be  disregarded  and  an  averaged  value  will  be  used).  This  is  not 
a  severe  limitation  in  the  analysis  of  the  ground  state  electronic  structure  of  the  trapped 
hole  center.  The  hyperfine  and  nuclear  electric  quadrupole  splittings  that  are  observed  in 
ENDOR  spectra  do  not  depend  on  the  g  matrix.  With  this  qualification,  focus  is  now 
shifted  to  determining  values  for  the  hyperfine  and  quadrupole  parameters. 

In  order  to  solve  the  spin  Hamiltonian  for  the  hyperfine  (A)  and  nuclear  electric 
quadrupole  ( Q )  matrices,  ENDOR  data  had  to  be  acquired.  To  completely  specify  the  A 
and  Q  matrices,  a  minimum  of  six  independent  data  sets  are  needed.  Independent  sets  of 
data  were  collected  on  different  days.  Included  were  sets  of  data  taken  with  the  magnetic 
field  parallel  to  several  high-symmetry  directions  of  the  crystal:  [001],  [100],  [110],  and 
[101].  These  data  sets  alone  yielded  enough  spectra  for  a  least  squares  fitting  routine  to 
be  used  to  determine  the  A  and  Q  matrices  in  the  spin  Hamiltonian.  Additional  data  were 
taken  at  angles  within  two  of  the  crystal  planes  to  increase  the  accuracy  of  the  least 
squares  fit  and  verify  the  results  of  the  fitting. 

In  the  ENDOR  experiment,  each  orientation  of  the  defect  in  the  crystal  will 
produce  a  spectrum  that  consists  of  six  peaks.  These  six  peaks  are  divided  into  two  sets 
of  three,  or  triplets.  One  set,  shown  in  Fig.  3.5,  consists  of  triplets  located  at  higher 
frequencies.  The  second  set  appears  at  lower  frequencies  near  4-8  MHz.  The  lower 
frequency  set  is  not  shown  nor  was  it  used  to  fit  the  spin  Hamiltonian.  The  lower 
frequency  ENDOR  spectra  overlap  other  spectra  that  appear  at  lower  frequencies.  These 
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overlapping  spectra  include  lines  due  to  the  nuclear  free  spins  of  lithium  and  boron. 

Since  each  defect  site  produces  a  triplet,  one  would  expect  to  see  all  four  sites  line  up 
when  the  magnetic  field  is  parallel  to  the  c  axis  of  the  crystal  and  there  to  be  two  sets  of 
triplets  when  the  magnetic  field  is  aligned  parallel  to  the  a  axis.  This  is  shown  by  the  top 
and  bottom  spectra  in  Fig.  3.5.  The  top  spectrum  has  two  sets  of  triplets  and,  though  the 
sets  labeled  sites  1  and  4  have  small  intensities,  their  positions  could  be  clearly  identified. 
The  remaining  spectra  in  Fig.  3.5  demonstrate  the  angular  dependence  of  the  site 
positions  in  the  crystal  as  the  direction  of  the  magnetic  field  is  varied.  The  site 
assignments  are  not  made  at  the  time  of  data  collection  but  occur  during  the  fitting 
process. 

The  spacing  within  the  triplets  is  the  result  of  nuclear  electric  quadrupole 
splitting.  Looking  at  one  of  these  triplets,  the  two  outer  peaks  are  separated  from  the 
center  peak  by  twice  the  effective  quadrupole  coupling  parameter  for  that  angle.  This 
will  be  further  discussed  later  in  this  chapter.  In  Fig.  3.5,  there  are  instances  where  the 
three  quadrupole  components  are  not  resolved,  i.e.,  a  single  line  is  observed  instead  of  a 
triplet.  For  example,  site  1  in  the  spectrum  labeled  “c  axis  +  80  deg  toward  a  axis”  has 
one  line  instead  of  three  lines.  A  single  peak  will  occur  when  the  effective  quadrupole 
moment  is  zero.  At  values  close  to  zero,  the  triplet  will  appear  as  a  single  broad  peak. 
This  occurs  when  the  linewidths  are  larger  than  the  separation  between  lines.  This  can  be 
seen  when  following  Site  2  in  Fig.  3.5  from  20  to  60  degrees  from  the  c  axis  toward  the  a 
axis.  During  the  initial  collection  of  data,  this  single  large  peak  was  puzzling  as  it  was 
not  yet  understood  that  it  was  actually  three  peaks  on  top  of  each  other. 


52 


Magnetic  Field 

Parallel  to:  Site  1  and  4  Site  2  and  3 


Site  1 


c  axis  +  60  deg  toward  a  axis 


Site  2  Site  3  and  4 


Site  j  Site  2  Site  3  and  4 


Site  1  Site  2 


Site  3  and  4 


10 


11 


t  1  t  E  l 

12 

Frequency  (MHz) 


F 

13 


14 


15 


Figure  3.5.  ENDOR  spectra  collected  as  the  magnetic  field  was  rotated  from  the  c  axis  to 
the  a  axis  in  20  degree  increments.  These  spectra  are  best  viewed  starting  from  the 
bottom  where  all  four  crystal  sites  are  equivalent  and  yield  a  single  spectrum. 


3.5.  Spin  Hamiltonian 

Swinney  et  al.  [31]  suggested  that  the  low  temperature  trapped  hole  center  in  LTB 
consisted  of  a  hole  with  spin  S  =  1/2  on  an  oxygen  ion  and  has  “slightly”  inequivalent 
hyperfine  interactions  with  two  neighboring  boron  nuclei.  The  following  spin 
Hamiltonian  describes  this  situation. 
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H  =  (3 S  ■  g  ■  B  +  Ix  A  ■  S  +  I2  ■  A  S  +  Ii  ■  Q  ■  Ix  +  I2  ■  Q  ■  I2 

~SnPn^i  '  B  —  §n^w^2  '  B  (3-4) 

In  this  spin  Hamiltonian,  each  boron  nucleus  interacting  has  its  own  hyperfine,  nuclear 
electric  quadrupole,  and  nuclear  Zeeman  term. 

During  the  present  investigation,  it  became  clear  that  the  slight  inequivalence 
reported  for  the  interactions  with  the  two  boron  nuclei  was  due  to  the  difficulty  of 
aligning  the  magnetic  field  along  the  crystal  axes.  Instead,  the  present  investigation 
shows  that  the  interactions  with  the  two  boron  nuclei  are  equivalent,  within  experimental 
error.  After  determining  that  the  two  boron  nuclei  are  equivalent,  the  spin  Hamiltonian 
used  to  describe  the  ENDOR  spectra  only  needs  to  be  written  for  one  nucleus. 

//  =  £S  g  B  +  I  A  S  +  I  Q  I  —  gN/?wI  ■  B  (3.5) 

With  an  isotropic  g  value  and  experimental  ENDOR  data,  Eq.  3.5  was  used  to  determine 
the  spin  Hamiltonian  parameters  describing  the  hyperfine  and  nuclear  electric  quadrupole 
matrices  for  the  low  temperature  trapped  hole  center  in  LTB. 

The  spin  Hamiltonian  is  “solved”  using  a  least-squares  fitting  routine.  The  fitting 
routine  diagonalizes  the  spin  Hamiltonian  matrix  and  iterates  through  each  parameter  that 
needs  to  be  determined.  The  parameters  are  incrementally  stepped  until  the  sum  of  the 
squares  is  minimized.  The  quantity  being  minimized  is  the  sum  of  the  differences 
(squared)  between  the  measured  and  calculated  ENDOR  frequencies.  The  fitting  routine 
determined  the  set  of  parameters  that  gave  a  “best  fit”  to  the  ENDOR  data  collected  when 
the  magnetic  field  was  aligned  parallel  to  the  [001],  [101],  [100],  and  [110]  crystal 
directions.  These  results  are  listed  in  Table  3.1.  Keep  in  mind  that  there  are  three 
principal  values  and  three  independent  Euler  angles.  The  Euler  angles  describe  the 
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rotations  needed  to  transform  the  given  parameter  from  the  lab  frame  (of  which  the  fitting 


parameter  is  solved)  to  the  parameter’s  principal  directions  where  its  representative 
matrix  is  diagonalized  with  the  principal  values  given. 


Table  3.1.  Principal  values  and  principal  angles  of  the  g,  A,  and  Q  matrices  for  the  hole 
center  in  LTB.  Values  for  g  are  assumed  isotropic.  There  are  no  principal  angles  listed 
for  the  g  matrix  since  it  is  assumed  isotropic. 


Parameter 

Principal  Value 

Principal  (Euler)  Angle 

g 

gx 

2.0125 

— 

gy 

2.0125 

— 

gz 

2.0125 

— 

Hyperfine 

Aj 

9.516 

23.6° 

A2 

14.549 

-14.6° 

a3 

20.094 

27.7° 

Quadrupole 

Qi 

-0.069 

46.5° 

q2 

-0.119 

71.8° 

Qb 

0.188 

0.1° 

Figure  3.6  shows  the  graphical  results  of  the  ENDOR  fitting  routine.  The  solid 
lines  represent  a  plot  of  the  spin  Hamiltonian  results  using  the  g.  A,  and  Q  matrix 
parameters  determined  by  the  fitting  routine.  The  blue  circles  represent  the  experimental 
data.  Four  of  the  eight  crystallographically  equivalent  sites  are  also  labeled  and  color 
coded  for  continuity.  A  consistent  set  of  site  assignments  were  determined  during  the 
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fitting  process.  Incorrect  assignments  result  in  large  values  for  the  sum  of  the  squares 


calculated  by  the  fitting  routine. 


Figure  3.6.  ENDOR  angular  dependence  plot.  Solid  lines  are  calculated  using  the  best- 
fit  parameter  values  and  open  circles  represent  experimental  data. 


3.6.  Discussion 

This  investigation  provided  several  interesting  results.  First,  an  explanation  has 
been  provided  for  the  single  peak  corresponding  to  Site  2  that  is  present  in  the  ENDOR 
spectra  as  the  magnetic  field  is  rotated  from  the  c  axis  to  the  a  axis.  The  quadrupole 
splitting  has  a  value  of  zero  for  two  different  orientations  of  the  magnetic  field  in  the  c  to 
a  plane.  This  causes  the  ENDOR  triplet  to  collapse  into  a  single  peak  at  these  points.  In 
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between  these  zero  points,  the  lines  comprising  the  triplet  are  close  together  and  remain 
unresolved.  Thus,  in  this  c  to  a  plane,  Site  2  appears  as  a  single  broad  line  for  60  degrees 
of  magnetic  field  rotation. 

Second,  the  two  boron  nuclei  adjacent  to  the  trapped  hole  have  nearly  identical 
hyperfine  matrices.  This  observation  was  initially  surprising,  as  there  is  no  fundamental 
reason  why  the  two  interactions  should  be  the  same.  The  evidence  that  these  two 
hyperfine  matrices  are  nearly  identical  is  quite  clear  when  the  EPR  and  ENDOR  spectra 
are  both  taken  into  consideration.  The  EPR  spectrum,  shown  in  Fig.  3.2(a),  contains 
seven  equally  spaced  lines  and  thus  demonstrates  that  two  boron  nuclei  are  involved.  On 
the  other  hand,  the  c-axis  ENDOR  spectrum,  shown  in  Fig.  3.6,  does  not  contain  resolved 
sets  of  lines  from  two  different  boron  nuclei.  Thus,  the  ENDOR  lines  from  the  two 
nuclei  must  coincide  for  this  orientation  of  magnetic  field.  Furthermore,  the  angular 
study  shows  that  the  ENDOR  lines  of  the  two  nuclei  overlap  for  all  angles.  This  proves 
that  the  hyperfine  matrices  for  the  two  boron  nuclei  are  equivalent  (i.e.,  the  sets  of 
ENDOR  lines  for  the  two  separate  nuclei  cannot  be  experimentally  resolved). 

The  nearly  equivalent  nature  of  the  two  boron  hyperfine  interactions  means  that 
the  principal  axes  associated  with  the  anisotropic  portion  of  each  matrix  must  have  the 
same  directions,  i.e.,  they  are  collinear.  In  terms  of  a  model,  recall  that  the  unpaired  spin 
resides  primarily  on  an  oxygen  ion  and  interacts  with  the  two  adjacent  boron  ions. 
Equivalent  boron  interactions  would  be  easily  explained  if  these  three  ions  were  to  lie 
collinear  in  the  crystal.  However,  in  the  LTB  lattice,  there  is  no  location  where  these 
three  ions  line  up.  This  suggests  that  lattice  relaxations  may  be  important  in  determining 
the  equilibrium  ground  state  structure  of  the  trapped  hole  center.  A  typical  lattice 
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relaxation  can  be  between  5  and  10  percent  of  a  lattice  spacing  and  this  movement  might 
be  sufficient  to  bring  an  oxygen  ion  more  in  line  with  its  two  boron  neighbors. 

The  next  step  is  to  suggest  a  specific  model  for  the  defect.  Originally,  Swinney  et 
al.  [31]  proposed  that  a  likely  model  had  the  hole  localized  on  an  0(4)  atom  (see  the 
basic  LTB  structure  in  Fig.  3.7(a)).  This  site  was  suggested  because  the  0(4)  atom  is  the 
only  oxygen  that  has  two  equivalent  boron  neighbors,  both  being  B(2)  ions.  Though  this 
model  could  be  correct,  there  are  several  concerns.  One  question  deals  with  the  boron- 
oxygen-boron  angle  which  is  108  degrees  for  the  0(4)  site  versus  the  approximately  180 
degree  angle  suggested  by  the  experimental  data.  A  second  questions  deals  with  the 
separation  distance  between  the  adjacent  stabilizing  lithium  vacancy  and  the  0(4)  oxygen 
ion  trapping  the  hole.  In  the  unrelaxed  lattice,  the  separation  distance  between  the  0(4) 
ion  and  the  lithium  vacancy  is  the  greatest  of  all  lithium-oxygen  separations. 

There  are  other  possible  oxygen  ions  where  the  hole  may  be  trapped.  It  could  be  trapped 
on  an  0(1)  ion  near  the  lithium  vacancy.  This  possible  hole  location  is  illustrated  in  Fig. 
3.7(b).  There  are  two  (B^O6-  units  shown  in  this  figure  along  with  an  adjacent  lithium 
ion.  The  figure  is  rotated  to  show  the  best  view.  This  0(1)  location  is  very  close  to  a 
lithium  vacancy.  This  is  the  second  closest  location  among  the  five  oxygen  ions  near 
each  lithium  vacancy.  Furthermore,  this  location  of  the  hole  has,  at  126  degrees,  the 
largest  angle  between  an  oxygen  ion  and  its  two  neighboring  boron  ions.  Although  this 
angle  is  still  much  less  than  180  degrees  in  the  unrelaxed  lattice,  it  does  not  necessarily 
eliminate  this  oxygen  ion  as  being  the  site  of  the  trapped  hole.  The  unique  principal  axes 
of  the  hyperfine  matrices  may  not  coincide  with  the  boron-oxygen  bond  directions.  The 
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Figure  3.7.  (a)  Schematic  representation  of  the  basic  (B4O9)6"  structural  building  block  in 
LFB4O7  crystals.  This  is  an  arbitrary  projection,  (b)  Possible  hole  trap  location  that  is 
different  from  the  model  suggested  by  Swinney  et  al.  [31]. 

spatial  distribution  of  the  /7-like  wave  function  of  the  hole  governs  the  principal  axes 
directions  of  the  hyperfine  matrices.  Without  a  full  ab  initio  quantum  chemistry 
calculation,  the  orientation  of  the  p  orbital  containing  the  hole  is  not  known.  The 
structure  of  LTB  makes  it  difficult  to  be  certain  of  the  trapped  hole  location. 

Similar  EPR  and  ENDOR  experiments  were  performed  in  2002  at  West  Virginia 
University  on  lithium  triborate  (LIB3O5)  crystals.  These  crystals  are  referred  to  as  LBO 
crystals.  The  investigation  was  part  of  the  PhD  dissertation  of  Wei  Hong  [39].  This 
study  Wei  Hong  identified  and  characterized  two  low  temperature  hole  centers  in  the 
LBO  crystals.  The  two  trapped  hole  centers  had  similar  g,  hyperfine,  and  nuclear  electric 
quadrupole  parameters.  The  thermal  stabilities  of  the  trapped  hole  centers  in  the  two 
materials  are  similar.  However,  the  trapped  hole  centers  in  the  two  materials  had  one 
major  difference  in  their  properties.  In  LBO,  only  one  boron  ion  interacted  strongly  with 
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the  trapped  hole  on  an  oxygen  ion,  whereas  in  LTB,  two  boron  ions  interact  strongly  with 
the  trapped  hole  on  an  oxygen  ion,  as  shown  in  this  dissertation.  A  representative 
spectrum  from  the  Wei  Hong  dissertation  is  provided  in  Fig.  3.8.  It  shows  the  four  EPR 
transitions  resulting  from  an  electron  spin  S  =  1/2  hole  interacting  with  a  1 1 B  nucleus  of 
nuclear  spin  I  =  3/2.  The  lithium  triborate  crystal,  with  orthorhombic  symmetry,  has  four 
crystallographically  equivalent  defect  orientations  compared  to  eight  orientations  found 
in  the  tetragonal  lithium  tetraborate  crystal.  In  LBO,  all  of  the  sites  are  equivalent  when 
the  magnetic  field  is  parallel  to  any  one  of  the  three  crystal  axes.  In  LTB,  site 
equivalence  only  occurs  when  the  magnetic  field  is  parallel  to  the  crystal’s  c  axis. 


3330  3355  3380  3405 

Magnetic  Field  (Gauss) 

Figure  3.8.  Representative  EPR  spectrum  from  a  trapped  hole  center  in  an  LBO  crystal. 
This  spectrum  was  taken  with  the  magnetic  field  parallel  to  the  c  axis  of  the  crystal. 
(Taken  from  the  PhD  Dissertation  of  Wei  Hong.)  [39] 

Similarities  between  the  LBO  and  LTB  crystals  are  found  in  their  structural  units. 
The  lithium  ions  are  crystallographically  equivalent  in  their  respective  lattices.  In  other 
words,  there  is  only  one  lithium  site  in  each  lattice.  The  BO4  and  BO3  structural  units  are 
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present  in  both  LBO  and  LTB.  In  an  LTB  crystal,  the  basic  building  block  is  a  (B4O9)6 
group  consisting  of  two  BO4  units  and  two  BO3  units.  In  contrast,  the  basic  building 
block  in  an  LBO  crystal  is  a  (B3O7)5  group,  consisting  of  one  BO4  unit  and  two  BO3 
units.  The  LTB  structure  is  shown  in  Fig.  3.7(a)  and  the  LBO  structure  is  shown  in  Fig. 
3.9  (a).  Both  structures  have  oxygen  ions  that  have  two  boron  neighbors  with  similar 

o  o 

oxygen-boron  bond  lengths,  between  1.35  A  and  1.50  A. 
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Figure  3.9.  (a)  Basic  structural  unit  of  an  LiBaOs  crystal.  Labels  identify  inequivalent  ion 
sites,  (b)  Proposed  model  for  the  hole  trap  in  LBO  (from  the  Wei  Hong  dissertation). 

In  the  LBO  crystal,  the  oxygen  atom  at  the  0(4)  site  was  determined  to  be  the 
location  of  the  trapped  hole.  The  sites  are  labeled  in  Fig.  3.9(a).  After  trapping  the  hole, 
the  neighboring  B(2)  ion  relaxes  significantly  toward  the  three  remaining  oxygen  atoms 
that  comprise  the  BO4  unit.  Figure  3.9(b)  shows  the  hole  location  with  one  bond 
removed  and  the  ions  at  their  relaxed  positions.  Gaussian  98  ab  initio  calculations  were 
also  performed  and  confirmed  the  EPR  experimental  data  [40].  Gaussian  is  an  ab  initio 
molecular  orbital  computer  program  that  uses  the  unrestricted  Hartree-Fock  method. 
When  using  this  program,  one  simulates  a  defect  by  constructing  a  large  molecule  that 
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represents  a  subset  of  ions  within  the  crystal  and  has  the  defect  near  the  center.  This 
simulated  molecule  contains  one  or  more  of  the  basic  building  blocks  of  the  crystal.  An 
appropriate  basis  set  is  chosen  and  the  program  then  calculates  the  electronic  energy  of 
the  ground  state.  A  geometry  optimization  provides  the  final  relaxed  positions  of  the  ions 
around  the  defect. 

Table  3.2  shows  the  g,  A,  and  Q  principal  values  obtained  from  the  trapped  hole 
centers  in  LBO.  These  were  obtained  from  angular  dependent  EPR  data  that  were  fit 
with  a  least  squares  routine.  There  are  two  items  to  note  in  the  principal  values  for  the 
LBO  hole  traps.  First,  the  g  values  in  LTB  are  similar  to  those  in  LBO  and  thus  support 
the  assignment  in  LTB  to  a  trapped  hole  center.  Second,  the  hyperfine  parameters  in 
LTB  are  approximately  half  of  those  determined  for  LBO. 

The  results  from  LBO  can  be  used  to  develop  a  plausible  model  for  LTB.  The 
borate  structures  are  expected  to  behave  in  a  similar  manner  in  the  two  materials.  In 
LBO,  the  oxygen  ion  traps  a  hole  and  causes  the  neighboring  boron  ion  to  relax  away 
from  the  hole  and  toward  the  remaining  three  oxygen  ions.  Applying  this  scenario  to 
LTB  suggests  that  the  0(1)  ion  (see  Fig.  3.7(b))  may  be  a  good  candidate  for  the  hole 
location.  If  the  0(1)  ion  contains  a  hole,  it  would  be  reasonable  to  have  the  neighboring 
boron  ions,  B(l)  and  B(2),  relax  away  from  the  hole  and  toward  their  remaining  oxygen 
ions.  This  may  allow  the  0(1)  ion  and  its  two  neighboring  boron  ions  to  form  a  linear 
chain.  This  would  explain  the  experimental  observation  that  LTB  has  equivalent 
hyperfine  interactions  with  two  boron  ions.  A  quantum  chemistry  molecular  orbital 
calculation  could  help  prove  or  disprove  this  suggested  model. 
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Table  3.2.  Principal  values  obtained  for  trapped  hole  centers  in  LBO.  Principal  axes 


directions  are  specified  by  polar  and  azimuthal  angles  9  and  cp.  (taken  from  the 
dissertation  of  Wei  Hong)  [39] 


Parameter 

Principal  Value 

Principal  Axis  Directions 

g  matrix 

0 

9 

gx 

2.0021 

50.5° 

259.1° 

gy 

2.0101 

62.4° 

143.6° 

gz 

2.0456 

51.9° 

29.4° 

Hyperfine  matrix 

Ai 

-34.49 

56.8° 

266.2° 

a2 

-38.53 

53.4° 

147.2° 

a3 

-19.84 

54.0° 

24.6° 

Quadrupole  matrix 

Qi 

-0.138 

O 

O 

'xf- 

350° 

q2 

-0.288 

82° 

250° 

Q3 

0.426 

51° 

153° 

This  research  project  has  provided  some  insight  to  the  use  of  lithium  tetraborate 
as  a  possible  neutron  detection  material.  Experiment  has  shown  that  there  are  large 
numbers  of  lithium  and  compensating  oxygen  vacancies  in  the  as-grown  crystals  [31]. 
These  vacancies  are  responsible  for  trapping  electrons  and  holes.  The  trapped  hole  center 
investigated  in  this  chapter  has  a  very  low  thermal  stability,  less  than  90  K,  and  is 
populated  during  exposure  to  ionizing  radiation.  These  trapped  holes  recombine  with 
trapped  electrons  and  release  visible  light  when  the  crystal  is  heated  above  90  K.  The 
hole  trapping  centers  could  have  both  a  positive  and  negative  effect  on  the  performance 
of  LTB  as  a  scintillating  material.  With  respect  to  a  positive  aspect,  the  oxygen  and 
lithium  vacancies  may  aid  in  scintillation  detection.  A  nuclear  event,  i.e.,  neutron 
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capture,  in  LTB  is  accompanied  by  ionizing  radiation.  This  ionizing  radiation  produces 
electrons  and  holes  that  will  then  “migrate”  through  the  crystal  lattice.  These  electrons 
and  holes  may  possibly  recombine  at  the  hole  trapping  site  investigated  in  this  chapter 
and  produce  visible  light  that  could  be  collected  in  a  photomultiplier  tube.  In  terms  of  a 
negative  effect,  ionizing  radiation  unrelated  to  neutron  capture  may  produce  electron  hole 
pairs  that  then  lead  to  a  “false”  count  in  a  neutron  detector  from  the  same  principles. 

Additional  work  is  needed  to  determine  the  feasibility  of  using  LTB  as  a  neutron 
detection  material.  EPR  has  provided  insight  into  the  electron  and  hole  centers  involving 
lithium  and  oxygen  vacancies.  Future  research  may  determine  if  a  post-growth  annealing 
process  could  be  utilized  to  reduce  the  lithium  and  oxygen  vacancy  concentrations.  Such 
a  process  has  been  developed  for  lithium  niobate  [41].  This  process  is  referred  to  as 
vapor  phase  equilibration  and  consists  of  heating  the  crystal  in  a  lithium-containing 
vapor.  Lithium  ions  enter  the  crystal  as  interstitials,  eliminate  the  lithium  vacancies,  and 
restore  the  crystal  to  its  stoichiometric  form.  If  a  post-growth  annealing  process  proved 
successful  in  LTB,  crystals  having  different  concentrations  of  vacancies  could  then  be 
exposed  to  a  neutron  source  and  the  resulting  detector  performance  evaluated. 

This  work,  characterizing  the  point  defects  in  lithium  tetraborate,  has  applications 
that  could  extend  beyond  neutron  detection.  Lithium  tetraborate  has  an  appreciable 
pyroelectric  coefficient  and  thus  could  be  used  as  a  thermal  sensor  [42-43].  LTB  also  has 
moderately  large  piezoelectric  constants  that  make  these  crystals  suitable  for  signal 
processing  and  frequency  control  devices  [44-46].  Acoustic  loss  peaks  associated  with 
point  defects  degrade  device  performance. 
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IV.  Perturbed  and  Unperturbed  Hole  Centers  (Ag2+  Ions)  in 
Silver-Doped  Lithium  Tetraborate  Crystals 
4.1.  Background 

The  introduction  to  this  dissertation  presents  the  idea  that  lithium  tetraborate 
doped  with  an  appropriate  impurity  could  be  utilized  in  radiation  detection  applications 
involving  scintillation.  Lithium  tetraborate  crystals  doped  with  Ag  are  excellent 
candidates  for  radiation-detection  applications  involving  light  emission  because  the  silver 
nuclei  trap  electrons  and  holes  that  recombine  above  room  temperature  and  radiate 
efficiently  [47-50].  Collaborative  research  conducted  by  the  Air  Force  Institute  of 
Technology  and  West  Virginia  University  identified  two  distinct  hole  centers  and  an 
electron  center  in  lithium  tetraborate  crystals  doped  with  silver  ions  [51].  These  centers 
were  identified  using  electron  paramagnetic  resonance,  after  the  crystals  had  been 
exposed  to  x-ray  radiation.  In  this  investigation,  thermoluminescence  showed  that  the 
electrons  and  holes  trapped  at  these  centers  recombine  above  room  temperature,  emitting 
visible  light  that  was  collected  using  a  photomultiplier  tube.  This  initial  work  was 
published  [51]  and  demonstrates  the  potential  for  using  silver-doped  lithium  tetraborate 
for  radiation  dosimetry  [52-54]  and  neutron  detection  involving  light  emission. 

However,  the  centers  participating  in  the  electron-hole  recombination  mechanisms  need 
to  be  more  completely  characterized  in  order  to  fully  develop  lithium  tetraborate  for  these 
optical  applications.  It  is  important  to  investigate  all  of  the  point  defects,  both  intrinsic 
and  extrinsic  (unintentional  or  intentional  impurities),  that  will  affect  the  performance  of 
LTB  as  a  radiation  detector. 
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The  research  of  Brant  et  al.  [51]  included  the  identification  of  two  Ag2+  hole 
centers  that  were  labeled  as  perturbed  and  unperturbed.  The  unperturbed  hole  center  was 
characterized  using  both  EPR  and  ENDOR  to  determine  the  g  and  hyperfine  matrix 
parameters.  This  center  was  determined  to  be  a  Ag2+  ion  occupying  a  lithium  site  in  the 
lattice  and  having  no  nearby  defects.  It  was  suggested  that  the  perturbed  center  was  a 
Ag2+  ion  sitting  at  a  lithium  site,  except  that  it  has  a  perturbing  ion  in  the  local 
environment.  For  example,  a  neighboring  sodium  or  silver  monovalent  ion  sitting  on  an 
adjacent  lithium  site  would  cause  the  slight  shifts  in  g  and  hyperfine  parameters 
compared  to  that  of  the  unperturbed  center.  On  the  other  hand,  the  electron  center  was 
identified  as  an  interstitial  silver  atom  (Ag°)  residing  near  a  boron  ion  at  a  regular  lattice 
position  [51].  The  unperturbed  Ag2+  hole  center  and  the  electron  center  were 
characterized  in  the  initial  investigation.  The  third  center,  the  perturbed  Ag2+  ion,  was 
not  characterized  at  that  time  because  its  EPR  spectrum  was  not  easily  resolved. 

Both  of  the  hole  centers  observed  in  silver-doped  LTB  are  the  result  of  a  hole 
trapped  on  a  silver  ion,  i.e.,  a  Ag2+  ion  is  formed.  Before  irradiation,  a  portion  of  the 
silver  ions  are  present  as  nonparamagnetic  Ag+  ions  (4d10)  replacing  Li+  ions.  After 
irradiation,  these  centers  trap  a  hole  and  convert  to  Ag2+  ions  (4d9)  with  a  spin  of  S  =  1/2. 
Silver  has  two  stable  isotopes,  107 Ag  and  109Ag.  Both  isotopes  are  approximately  50% 
abundant  (107Ag  =  5 1 .83%  and  109 Ag  =  48. 17%)  and  have  a  nuclear  spin  of  I  =  1/2.  They 
also  have  similar  nuclear  magnetic  moments,  -0.1  135  (107Ag)  and  -0.1305  (109Ag).  The 
magnetic  moments  here  are  expressed  in  “units”  of  nuclear  magnetons,  i.e.,  they  are  the 

ratio  ( — )  of  their  moments  to  that  of  the  nuclear  magneton.  For  each  silver  isotope,  the 
g-n 

interaction  of  an  electron  of  spin  1/2  with  a  nucleus  of  spin  I  =  1/2  produces  an  EPR 
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spectrum  consisting  of  two  lines.  Thus,  for  one  Ag2+  ion,  there  will  be  two  doublets,  one 
associated  with  each  silver  isotope,  with  slightly  different  separations  as  a  result  of  their 
different  magnetic  moments.  However,  the  nearly  equal  natural  abundances  and  the 
nearly  equal  magnetic  moments  of  these  two  isotopes  cause  the  EPR  hole-center 
spectrum  associated  with  107 Ag  and  109 Ag  to  appear  as  a  single  broadened  doublet.  In 
other  words,  the  individual  doublets  due  to  the  107Ag  and  109Ag  nuclei  are  not  resolved  in 
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Fig.  4.1.  Energy  level  diagram  for  the  Ag2+  hole  centers.  The  energy  levels  associated 


with  the  two  isotopes  are  illustrated  and  the  allowed  EPR  and  ENDOR  transitions  are 
shown  where  the  solid  red  lines  correspond  to  109Ag  and  the  dotted  orange  lines 

correspond  to  1<)7Ag. 
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EPR.  The  energy  level  diagram  shown  in  Fig.  4.1  illustrates  the  lack  of  resolution  in  the 
EPR  experiment  and  the  increased  resolution  in  the  ENDOR  experiment.  The  increased 
resolution  is  demonstrated  in  Fig.  4. 1  under  the  ENDOR  portion. 

As  shown  in  the  right  portion  of  Fig.  4.1,  electron-nuclear  double  resonance 
(ENDOR)  can  resolve  the  hyperline  lines.  Because  A/2  is  larger  than  Vn  for  the  Ag2+ 
ions,  each  silver  isotope  produces  an  ENDOR  spectrum  consisting  of  two  peaks  centered 
on  A/2  and  separated  by  2vn.  Each  Ag  isotope  has  a  different  value  of  Vn  (.606  Mhz  for 
l07Ag  and  .698  MHz  for  109 Ag);  this  leads  to  a  different  separation  between  the  ENDOR 
doublets  and  aids  in  identification  of  silver  as  the  responsible  nuclei.  In  Fig.  4.2,  the 
upper  blue  spectrum  consists  of  two  doublets.  The  doublet  between  32  and  34  MHz  has  a 
separation  close  to  1.4  MHz  and  can  be  assigned  to  1<)9Ag,  while  the  doublet  between  28 
and  30  MHz  has  a  separation  close  to  1.2  MHz  and  can  be  assigned  to  Ag.  A  second 
aid  in  identifying  silver  as  the  responsible  nuclei  resides  in  the  mid-points,  or  A/2  values, 
of  the  two  doublets  observed  in  the  ENDOR  spectrum.  Both  vn  and  A  are  proportional  to 
the  magnetic  moment.  The  ratio  obtained  by  dividing  the  A/2  measured  for  the  109Ag  by 
the  A/2  measured  for  the  Ag  nucleus  should  match  the  ratio  of  their  magnetic 
moments,  or  1.152.  The  ENDOR  spectrum  from  the  perturbed  Ag2+  hole  center, 
presented  in  the  lower  portion  of  Fig.  4.2,  shows  two  transitions  at  21.34  and  20.23  MHz 
for  109Ag  and  two  transitions  at  18.47  and  17.48  MHz  for  107 Ag.  The  mid-points  are  thus 
20.785  and  17.975  MHz,  respectively,  and  have  a  ratio  of  1.156.  This  ratio  is  in  close 
agreement  with  the  ratio  of  the  Ag  magnetic  moments. 

The  earlier  study  performed  by  Brant  et  al.  [51]  characterized  the  unperturbed 
Ag2+  hole  center.  Angular  dependence  EPR  data  was  collected  to  determine  the  g  matrix 
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and  ENDOR  data  was  collected  to  determine  the  hyperfine  matrix  for  each  isotope.  In 
this  study,  a  similar  approach  was  attempted  for  the  perturbed  Ag2+  hole  center. 

However,  there  were  several  complications  that  made  characterizing  the  perturbed  Ag2+ 
hole  center  more  difficult.  First,  the  intensity  of  the  perturbed  center  in  EPR  spectra  is 
less  than  half  of  the  intensity  of  the  unperturbed  center.  This  is  due  to  the  lower 
concentration  of  the  perturbed  defect  compared  to  that  of  the  unperturbed.  Also,  the 
lower  intensity  of  the  EPR  signal  correlates  with  a  lower  intensity  of  ENDOR  signals. 

An  even  greater  problem  arises  because  the  spin-lattice  relaxation  times  associated  with 
the  perturbed  Ag2+center  are  not  conducive  to  strong  ENDOR  transitions.  Another 
complication  encountered  while  characterizing  the  perturbed  Ag2+  hole  center  is  a  lack  of 
clearly  observable  EPR  lines  at  a  variety  of  orientations  of  the  magnetic  field,  i.e.,  the 
perturbed  center  EPR  lines  are  often  under  the  more  prominent  EPR  signals  from  the 
unperturbed  center.  There  are  eight  independent  sites  in  the  crystal  and  thus  eight 
orientations  of  a  given  defect.  When  the  magnetic  field  is  rotated  in  a  high  symmetry 
plane  of  the  crystal,  the  tetragonal  structure  of  LTB  reduces  the  number  of  sites  observed 
in  an  EPR  spectrum  to  four  or  less.  The  lines  of  the  perturbed  Ag2+  sites,  are  often  times 
under  the  unperturbed  Ag2+  lines.  Thus,  it  is  difficult  to  observe  a  large  number  of 
transitions  of  the  perturbed  Ag2+  hole  center  when  the  magnetic  field  is  parallel  to  many 
of  the  directions  in  the  crystal. 

Attempts  were  repeatedly  made  in  the  current  study  to  obtain  a  complete  set  of 
ENDOR  data  for  the  perturbed  Ag2+  center.  Success  was  not  achieved.  Much  of  this 
effort  involved  trying  to  determine  the  optimum  ENDOR  spectrometer  settings.  These 
included  temperature,  modulation  frequency,  modulation  amplitude,  RF  power,  and 
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microwave  power.  The  ENDOR  spectrum  from  the  perturbed  Ag2+  ions  was  always 
much  lower  in  intensity  than  that  from  the  unperturbed  Ag2+  ions.  The  difference  in 
intensities  of  the  ENDOR  spectra  from  the  unperturbed  and  perturbed  Ag2+  hole  center  is 
shown  in  Fig.  4.2.  The  lower  frequency  signals  of  the  upper  spectra  are  from  the 


Fig.  4.2.  ENDOR  spectra  collected  for  the  unperturbed  (blue  and  upper  right)  and 
perturbed  (black  and  upper  left)  Ag  hole  centers.  These  spectra  were  collected  during 
the  same  experiment,  but  while  the  magnetic  field  was  “parked”  on  each  respective  EPR 
transition,  i.e.,  two  separate  data  collections.  The  lower  trace  shows  the  magnified 
spectrum  of  the  perturbed  Ag“  hole  center. 
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perturbed  Ag2+  center  and  the  higher  frequency  signals  in  this  spectrum  (shown  in  blue) 
are  from  the  unperturbed  Ag2+  center.  In  this  figure,  the  ENDOR  spectrum  for  the 
unperturbed  Ag2+  hole  center  is  eight  times  greater  in  magnitude  compared  to  the 
perturbed  Ag2+  hole  center.  The  lower  portion  of  Fig.  4.2  shows  an  expanded  and 
magnified  view  of  the  ENDOR  spectrum  of  the  perturbed  Ag2+  center.  Here,  all  four 
ENDOR  transitions  associated  with  the  two  silver  isotopes  are  easily  identified.  During 
this  investigation,  most  of  the  ENDOR  peaks  could  be  seen  above  the  noise  when  the 
magnetic  field  was  rotated  between  the  c  axis  and  the  a  axis  of  the  crystal.  However, 
several  attempts  failed  to  find  ENDOR  lines  above  the  noise  when  the  magnetic  field  was 
rotated  in  the  basal  plane  from  an  a  axis  to  an  a  axis  of  the  crystal.  Most  likely,  this  is  a 
result  of  the  dependence  of  the  spin-lattice  relaxation  time  on  the  direction  of  the  external 
static  magnetic  field  in  the  crystal.  After  several  attempts,  it  became  clear  that  the 
collection  of  ENDOR  spectra  in  the  basal  plane  of  the  crystal  would  not  be  possible  using 
the  crystals  that  were  available.  If  a  full  set  of  ENDOR  were  collected,  it  would  have 
been  possible  to  determine  the  hyperfine  parameters  for  each  silver  isotope  just  as  it  was 
for  the  unperturbed  Ag2+  hole  center. 

An  angular  dependence  study  using  only  the  EPR  spectra  collected  in  the  high 
symmetry  crystal  planes  is  sufficient  to  determine  the  g  and  hyperfine  parameters.  The 
hyperfine  parameters  determined  from  the  EPR  spectra  will  be  an  average  of  the 
parameters  for  the  separate  silver  isotopes.  By  using  only  EPR  spectra  for  data 
collection,  it  is  still  possible  to  get  a  full  set  of  data  and  calculate  all  of  the  parameters  of 
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Fig.  4.3.  EPR  spectra  of  perturbed  (P)  and  unperturbed  (U)  Ag2+  ions  in  LTB. 
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the  spin  Hamiltonian.  Some  accuracy  may  be  lost  due  to  the  higher  resolution  of  the 
ENDOR  experiment. 

Another  complication  encountered  in  characterizing  the  perturbed  Ag2+  hole 
center  was  having  enough  independent  data  points  for  a  fit  of  the  spin  Hamiltonian. 

Figure  4.3  shows  EPR  spectra  collected  when  the  magnetic  field  is  parallel  to  high 
symmetry  directions  of  the  crystal,  as  well  as  the  [101]  crystal  direction.  In  each 
spectrum,  transitions  are  labeled  with  a  U  if  the  transition  corresponds  to  an  unperturbed 
Ag2+  hole  center  transition  or  a  P  if  the  transition  corresponds  to  a  perturbed  Ag2+  hole 
center  transition.  The  g  and  “average”  hyperfinc  parameters  can  be  found  for  the 
unperturbed  Ag2+  hole  center  using  the  four  spectra  shown  in  Fig.  4.3.  However,  the  g 
and  hyperfine  parameters  cannot  be  precisely  determined  using  the  same  spectra  for  the 
perturbed  Ag2+  hole  center.  When  the  magnetic  field  is  parallel  to  the  crystal  [101] 
direction,  several  of  the  perturbed  Ag2+  hole  center  peaks  are  obscured  by  the 
unperturbed  Ag2+  hole  center  peaks.  For  example,  in  Fig.  4.3  (magnetic  field  along  [101] 
direction),  there  is  a  doublet  centered  on  2850  Gauss  that  corresponds  to  unperturbed 
Ag2+  hole  center  transitions.  Upon  close  observation  this  doublet  is  distorted  compared 
to  the  unperturbed  double  centered  at  3150  Gauss.  This  distortion  of  the  low-field 
doublet  is  due  to  an  underlying  doublet  associated  with  the  perturbed  Ag2+  hole  center. 
Because  the  spectra  of  Fig.  4.3  do  not  provide  enough  independent  data  for  a  fitting  of  the 
spin  Hamiltonian  of  the  perurbed  Ag2+  hole  center,  a  full  angular  study  was  conducted, 
and  provided  sufficient  independent  data  to  allow  a  good  fit  of  the  spin  Hamiltonian 
describing  the  perturbed  Ag2+  hole  center. 
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4.2.  EPR  Results 


Since  ENDOR  was  unable  to  provide  the  hyperfine  parameters  for  the  perturbed 
center,  it  was  decided  to  first  use  EPR  spectra  alone  to  re-determine  the  parameters  of  the 
spin  Hamiltonian  describing  the  unperturbed  Ag2+  hole  center.  Then,  if  the  parameters 
obtained  using  EPR  spectra  alone  for  the  unperturbed  Ag2+  hole  center  are  in  good 
agreement  with  the  previous  parameters  obtained  using  both  EPR  and  ENDOR  [51], 
there  is  confidence  that  the  perturbed  Ag2+  hole  center  parameters  obtained  from  the  EPR 
spectra  will  be  accurate.  Data  were  taken  every  10  degrees  as  the  magnetic  field  was 
rotated  from  the  c  to  the  a  axis  of  the  crystal.  Also,  data  were  taken  in  this  plane  for  the 
[101]  direction.  Representative  spectra  taken  in  the  c  to  a  plane  are  shown  in  Fig.  4.4. 
The  larger  peaks  in  this  figure  correspond  to  the  unperturbed  Ag2+  hole  center  and  can  be 
clearly  seen.  The  smaller  peaks,  corresponding  to  the  perturbed  hole  center,  are  present 
yet  harder  to  identify  and  are  marked  with  arrows.  In  some  cases,  the  spectra  from  the 
two  centers  overlap  thus  making  identification  difficult.  For  example,  the  spectrum 
labeled  “c  axis  +  10  deg  twd  a  axis”  in  Fig.  4.4  contains  closely  space  and  overlapping 
lines.  In  order  to  obtain  a  complete  set  of  spin  Hamiltonian  parameters  for  the  perturbed 
Ag2+  hole  center,  there  needs  to  be  enough  independent  data  to  allow  a  fitting  that  solves 
for  the  six  parameters  that  describe  each  term  in  the  spin  Hamiltonian. 

After  taking  EPR  spectra  in  more  than  20  independent  directions,  four  sets  of 
spectra  were  identified  where  a  full  set  of  transitions  related  to  the  perturbed  Ag2+  hole 
center  could  be  identified.  These  directions  are  shown  in  Fig.  4.5.  When  the  magnetic 
field  is  parallel  to  the  c  axis  of  the  crystal,  the  eight  orientations  of  the  defect  that  exist  in 
the  unit  cell  of  the  crystal  are  all  crystallographically  equivalent  and  a  single  pair  of  EPR 
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Fig.  4.4.  Angular  dependence  data  collected  while  the  magnetic  field  was  rotated 
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Fig.  4.5.  Four  independent  sets  of  data  where  all  of  the  EPR  transitions  for  the  Ag  hole 
center  can  be  identified.  Transitions  for  the  perturbed  Ag2+  hole  center  are  identified 

with  arrows. 
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transitions  is  observed.  When  the  magnetic  field  is  along  the  a  axis  and  [110]  directions 
of  the  crystal,  the  eight  orientations  have  two  equivalences  thus  producing  four  EPR 
transitions.  Finally,  when  the  magnetic  field  is  in  a-c  plane  of  the  crystal  80  degrees  from 
c  axis  toward  the  a  axis,  there  were  four  equivalences  and  eight  identifiable  EPR 
transitions.  Arrows  are  drawn  in  Fig.  4.5  to  show  the  identified  transitions.  When  the 
magnetic  field  was  along  the  [110]  direction  there  is  a  transition  that  is  both  a  perturbed 
and  unperturbed  Ag2+  hole  center  transition  at  3175  G.  The  data  sets  of  Fig.  4.5  can  be 
used  to  fit  the  spin  Hamiltonian.  To  increase  the  accuracy  of  the  fit,  106  transitions  that 
could  be  identified  for  each  angle  of  magnetic  field  were  used  compared  to  the  1 8  that 
could  have  been  used. 

With  enough  data  to  solve  for  both  the  perturbed  and  unperturbed  Ag2+  hole 
center,  the  unperturbed  center  was  solved  first  in  order  to  compare  to  the  published 
results.  The  parameters  determined  from  the  EPR  spectra  collected  for  this  dissertation 
are  shown  in  Table  4.1.  These  parameters  were  determined  using  a  least-squares  fitting 
routine  to  find  the  best  “fit”  to  the  spin  Hamiltonian  that  models  the  unperturbed  Ag2+ 
hole  center.  The  published  results  of  Brant  et  al.  [51]  are  provided  in  Table  4.2.  The 
published  data  was  solved  using  the  same  least  squares  fitting  routine  with  the  exception 
that  the  hyperfine  parameters  for  each  isotope  were  calculated  individually;  in  other 
words,  the  fitting  routine  was  run  separately  to  solve  for  107 Ag  and  109Ag.  In  the  current 
study,  an  average  hyperfine  was  determined  due  to  the  lack  of  ENDOR  data.  The 
hyperfine  principal  values  are  the  only  data  that  are  different  and  since  the  two  isotopes 
are  52%  and  48%  abundant,  respectively,  an  average  of  the  values  is  provided  for 
reference  to  Table  4.1.  The  angular  dependence  associated  with  the  two  sets  of 
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parameters  are  plotted  in  Fig.  4.6.  The  blue  circles  indicate  data  points  taken  from  EPR 
spectra  acquired  in  the  current  study.  The  blue  dashed  line  represents  a  plot  of  the 
spin  Hamiltonian  using  the  parameters  that  were  calculated  from  the  EPR  alone  (i.e.,  the 
Table  4.1  parameters).  The  solid  red  line  represents  a  plot  of  the  spin  Hamiltonian  using 
the  parameters  determined  by  Brant  et  al.  from  EPR  and  ENDOR  (i.e.,  the  Table  4.2 
parameters).  The  averaged  values  of  the  hyperfine  parameters  were  used  when 
generating  these  solid  red  lines.  Both  sets  of  parameters  follow  the  data  closely  and  the 
largest  difference  between  the  two  was  found  to  be  six  Gauss.  The  widths  of  the 
individual  lines  in  the  EPR  spectra  for  the  unperturbed  center  were  at  least  that  value  and 
tended  to  be  closer  to  8.5  Gauss.  It  can  be  concluded  that  the  two  sets  of  data  are 
accurate  within  a  line  width  of  an  EPR  trace,  which  is  excellent.  As  further  illustration  of 
experimental  accuracy,  an  error  of  one  degree  of  angular  rotation  of  the  magnetic  field 
causes  a  shift  of  six  to  eight  Gauss  when  conducting  the  EPR  experiment  for  silver-doped 
LTB. 


Table  4.1.  Spin-Hamiltonian  parameters  calculated  using  EPR  spectra  alone  for  the 
unperturbed  Ag2+  hole  center. 

Parameter  Principal  Value  Principal  (Euler)  Angles 


8 


gx 

2.0427 

40. 0‘ 

gy 

2.1263 

18. 7‘ 

gz 

2.3677 

65. 0‘ 

Hyperfine 

Ai 

45.05  MHz 

29.  r 

a2 

65.89  MHz 

13.4‘ 

a3 

82.88  MHz 

84.  r 
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Table  4.2.  Spin-Hamiltonian  parameters  calculated  using  EPR  and  ENDOR  for  the 
unperturbed  Ag2+  hole  center,  (from  Brant  et  al.)  [51] 


Parameter 

Principal 

Value 

107Ag 

Principal 

Value 

109Ag 

Principal 
(Euler)  Angles 

g 

gx 

2.0439 

40.2 

gy 

2.1294 

19.0 

gz 

2.3708 

64.4 

Hyperfine 

(AVERAGE) 

Ai 

45.75  MHz 

42.30 

49.19 

31.0 

A2 

67.13  MHz 

62.20 

72.07 

15.6 

a3 

83.09  MHz 

77.05 

89.12 

76.1 

As  a  result  of  the  agreement  between  the  two  methods  (EPR  alone  versus 
combined  EPR  and  ENDOR)  for  determining  the  parameters  of  the  unperturbed  Ag2+ 
hole  center,  the  perturbed  Ag2+  hole  center  was  fit  using  data  collected  from  the  same 
EPR  spectra.  EPR  transitions  (or  peaks)  taken  from  the  full  angular  dependence  data 
were  used  as  the  experimental  points  in  the  spin-Hamiltonian  fitting  routine.  The  “best- 
fit”  parameters  of  the  spin  Hamiltonian  are  presented  in  Table  4.3.  The  perturbed  Ag2+ 
hole  center  angular  dependence  is  shown  graphically  in  Fig.  4.7.  The  solid  lines  are 
calculated  using  the  best-fit  parameter  values  and  the  open  circles  represent  the 
experimental  data.  The  lower  pair  of  lines,  between  the  [001]  and  [100]  direction,  has 
comparatively  fewer  data  points  than  the  other  three  pairs  of  lines.  Transitions  due  to 
these  orientations  of  the  defect  could  not  be  determined  due  to  an  overlapping  set  of 
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Table  4.3.  Spin-Hamiltonian  parameters  calculated  using  only  EPR  for  the  perturbed 
Ag2+  hole  center. 

Parameter  Principal  Value  Principal  (Euler)  Angles 


g 


gx 

2.0319 

35.5° 

gy 

2.1212 

9.7° 

gz 

2.3604 

-9.6° 

Hyperline 

Ai 

28.18  MHz 

35.2° 

A2 

35.81  MHz 

25.9° 

A3 

65.12  MHz 

— 38.7C 

Fig.  4.6.  Angular  dependence  comparison  of  the  unperturbed  Ag2+  hole  center  where  the 
dashed  blue  lines  represent  fit  parameters  determined  from  EPR  spectra  alone  are  plotted 
along  with  solid  red  lines  that  represent  parameters  detennined  using  EPR  and  ENDOR. 
Blue  dots  represent  experimental  data  points. 


80 


peaks  from  the  unperturbed  Ag2+  hole  center.  Squares  are  used  in  Fig.  4.7  to  represent 
the  unperturbed  peak  positions  that  covered  the  perturbed  Ag2+  hole  center  transitions. 

There  is  one  last  important  item  discovered  during  the  research  for  this  chapter. 
In  order  to  fit  the  spin  Hamiltonian  using  the  least- squares  fitting  routine,  site 
assignments  need  to  be  made  relative  to  the  data.  This  takes  place  when  the  data  is 
entered  into  the  fitting  routine.  In  the  first  plane  of  the  fitting  routine,  or  for  the  [001]  to 
the  [100]  direction,  site  assignments  can  be  randomly  assigned.  Each  color  in  Fig.  4.7 
corresponds  to  an  individual  site  assignment,  hence  eight  different  colors.  The  sites  can 


Fig.  4.7.  EPR  angular  dependence  plot.  Solid  lines  are  calculated  using  the  best-fit 
parameter  values  and  open  circles  represent  experimental  data.  Each  color  corresponds 
to  a  specific  site  assignment,  one  through  eight. 
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then  only  take  on  certain  values  between  the  [100]  and  [010]  directions.  For  example  the 
red  line  in  Fig.  4.7  may  correspond  to  site  1  in  the  crystal.  The  red  line  goes  from  the 
[001]  plane  to  the  upper  set  of  transitions  along  the  [100]  direction.  From  the  [100] 
direction  the  red  line  can  only  follow  one  of  the  four  “paths”  leading  away  from  the  upper 
set.  While  some  choices  reduced  the  fitting  accuracy  greatly,  there  was  a  second  set  of 
site  assignments  that  gave  results  close  to  those  that  gave  the  lowest  sum  in  the  least- 
squares  fit.  There  were  thus  two  sets  of  parameters  that  closely  fit  the  two  planes  of  data 
being  used  in  the  least-squares  fitting  routine.  For  some  crystal  symmetries,  it  is 
impossible  to  create  site  assignments  that,  in  turn,  allow  for  two  sets  of  parameters  to  fit 
the  data  for  two  planes  of  the  crystal.  In  this  study,  it  was  mistakenly  believed  to  be  the 
case  for  tetragonal  symmetry. 

When  the  two  different  sets  of  parameters  were  plotted,  the  lines  overlapped  each 
other  for  the  two  planes  in  Fig.  4.7.  However,  when  a  third  plane  was  plotted  the  two  sets 
of  data  diverged  and  showed  EPR  transitions  at  different  magnetic  fields.  It  was 
therefore  decided  that  data  must  be  taken  in  this  third  plane  to  confirm  that  the  correct 
site  assignments  were  being  made.  Data  were  then  taken  by  rotating  the  magnetic  field 
between  the  [110]  direction  to  the  [001]  direction  and  confirmed,  in  this  case,  that  the 
lowest  sum  in  the  least-squares  fit  yielded  the  correct  parameters.  The  fact  that  two  sets 
of  parameters  can  yield  the  same  fitting  for  two  planes  has  shown  that  a  third  set  of  data 
must  be  taken  for  future  cases  involving  crystals  with  tetragonal  symmetry  when  an 
angular  dependence  study  is  going  to  be  performed.  The  third  plane  of  data  is  not 
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presented  in  this  dissertation  as  it  does  not  improve  the  accuracy  of  the  spin-Hamiltonian 
parameters,  but  does  confirm  correct  site  assignments. 

4.3.  Discussion 

In  order  for  lithium  tetraborate  to  become  a  viable  candidate  for  dosimetry  and 
nuclear  detection  applications  [55-56],  it  is  important  to  characterize  all  of  the  defects  in 
the  material.  The  g  and  hyperfine  values  for  the  perturbed  Ag2+  hole  center  have  not 
been  previously  reported.  This  chapter  has  shown  that  it  is  possible  to  find  a  good  set  of 
g  and  hyperfine  parameters  using  EPR  alone  for  a  defect  that  contains  two  isotopes  of 
similar  nuclear  magnetic  moments  that  could  not  otherwise  have  been  characterized. 

The  g  matrix  principal  values  determined  in  this  research  are  in  good  agreement 
with  results  reported  in  earlier  studies  of  Ag2+  ions  in  other  crystals  [57-59].  The  values 
of  the  g  and  hyperfine  matrices  do  not,  however,  provide  enough  information  for  an 
improved  model  that  clearly  establishes  the  difference  between  the  perturbed  and 
unperturbed  Ag2+  hole  center.  The  values  do  support  the  assignment  of  “perturbed”  and 
“unperturbed”  labels.  The  unperturbed  center  is  posited  to  be  a  substitutional  ion,  i.e.,  a 
silver  ion  that  sits  at  a  lithium  site  in  the  crystal.  Since  silver  is  a  dopant  and  similar  to 
lithium  in  that  it  has  a  single  unpaired  electron  in  its  outer  shell,  it  would  appear 
statistically  favorable  to  be  a  substitutional  ion  with  no  neighboring  defect.  This  is 
explained  by  the  higher  signal  intensity  of  the  unperturbed  center  relative  to  the  perturbed 
center  in  the  EPR  spectra.  The  perturbed  center,  with  its  lower  signal  intensity  and 
population,  would  then  be  likely  to  have  a  neighboring  defect.  This  perspective  is 
reinforced  with  the  hyperfine  parameters  determined  for  the  perturbed  Ag2+  hole  center, 
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since  they  are  lower  than  the  hyperfine  parameters  of  the  unperturbed  Ag2+  hole  center. 

If  a  hole  is  trapped  on  a  silver  ion  with  a  neighboring  monovalent  silver  ion,  smaller 
hyperfine  values  appear  appropriate.  It  has  also  been  suggested  [51]  that  this  perturbation 
is  due  to  Na+  ions  occupying  adjacent  substitutional  lithium  sites.  This  suggestion  is 
based  upon  neutron  activation  results  that  showed  concentrations  of  sodium  present  in 
LTB  crystals  (those  used  in  the  present  study)  with  intentionally  large  concentrations  of 
silver  [31].  Sodium  ions  could  cause  a  slight  perturbation  to  the  g  and  hyperfine  values, 
perhaps  with  no  other  effect  on  the  spectra,  i.e.,  further  splitting  due  to  nuclear  spin 
interactions  from  sodium.  Also,  although  unlikely,  this  nearest  silver  could  be  an 
interstitial  monovalent  silver  ion  that  is  responsible  for  the  Ag°  electron  center  seen  in  the 
EPR  spectra  of  x-ray  irradiated  crystals. 

Future  research  can  help  determine  the  impact  of  these  defects,  and  techniques 
can  be  developed  to  reduce  or  magnify  the  defect  density  dependent  upon  its  usefulness 
in  a  given  application.  Currently,  copper  appears  in  the  literature  as  the  preferred  dopant 
candidate  in  lithium  tetraborate  but  the  work  done  by  Brant  et  al.  [51]  and  this  latest 
research  effort  has  shown  silver-doped  lithium  tetraborate  to  be  a  strong  candidate  for 
dosimetry.  The  perturbed  and  unperturbed  Ag2+  hole  centers  appear  to  recombine  with 
electrons  from  the  electron  center  that  was  previously  identified.  This  recombination 
explains  the  strong  thermoluminescent  peak  that  is  produced  when  x-ray  irradiated  silver- 
doped  lithium  tetraborate  is  heated  above  room  temperature.  This  research  provides  the 
explanation  for  the  thermoluminescent  curve.  The  advantage  silver  appears  to  have  as  a 
dopant  versus  copper  in  lithium  tetraborate  is  that  silver-doped  lithium  tetraborate  has 
only  one  thermoluminescent  peak,  whereas  copper-doped  lithium  tetraborate  has  two 
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peaks.  This  chapter  has  been  the  final  set  of  characterization  needed  to  explain  the  defect 
sites  responsible  for  electron-hole  recombination  in  lithium  tetraborate  doped  with  silver. 
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V.  Identification  of  Radiation-Induced  Point  Defects  in  Copper- 
Doped  Lithium  Tetraborate  Crystals 

5.1.  Introduction 

Transition  metals  ions  have  been  doped  into  lithium  tetraborate  crystals  for  their 
light  emitting  properties  [60-62].  Specifically,  copper,  manganese,  and  silver  have  been 
the  subject  of  numerous  studies  [63-70].  These  controlled  dopants  have  been  shown  to 
be  highly  efficient  with  respect  to  their  radiative  emissions.  Copper-doped  LTB  has  been 
the  subject  of  an  EPR  investigation  in  the  literature.  Corradi  et  al.  [71]  fully 
characterized  a  Cu"  trapped  hole  center  by  performing  a  complete  angular  dependence 
study  from  which  the  parameters  for  the  g  and  hyperfine  matrices  were  determined.  This 
hole  center  was  modeled  as  an  S  =  1/2  hole  trapped  on  a  copper  ion  substituting  for  a 
lithium  ion  and  stabilized  by  a  neighboring  lithium  vacancy,  thus  creating  a  Cu  hole 
center.  This  center  investigated  by  Corradi  et  al.  will  be  referred  to  as  the  Cu  -Vy  hole 
center. 

2d- 

Unlike  silver-doped  crystals  where  there  are  no  Ag“  centers  prior  to  x-ray 

2d_ 

irradiation,  copper-doped  LTB  has  the  Cu'  -Vy  hole  center  always  present,  i.e.,  room 
temperature  x  rays  are  not  required  to  populate  the  Cu“  -Vy  hole  center.  The  EPR  signal 

2d- 

from  the  Cu"  -Vy  hole  center  has  a  significant  magnitude  without  x-ray  exposure,  and  is 
comparable  to  the  magnitude  of  the  Ag"  hole  centers  in  silver-doped  LTB  after  x-ray 
exposure.  An  initial  goal  of  this  study  was  to  verify  the  spin  Hamiltonian  parameters 
detennined  by  Corradi  et  al.  [71]  for  the  Cu'  -Vy  hole  center.  However,  before  data 
collection  began,  the  LTB  crystal  was  exposed  to  x  rays  to  further  magnify  and  populate 
the  Cu'  -Vy  hole  center  and  a  new  effect  was  observed  that  changed  the  scope  of  this 
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investigation.  Surprisingly,  x-ray  exposure  did  not  increase  the  signal  intensity  of  the 
Cu  -Vu  hole  center,  but  rather  created  several  new  spectra  that  have  not  been  previously 
reported  in  the  literature.  Thus,  this  investigation  has  focused  on  two  new  trapped 
electron  centers  and  two  new  trapped  hole  centers.  Their  identification  and 
characterization  are  the  subject  of  this  chapter.  The  two  electron  centers  have  isotropic 
EPR  spectra.  The  spectrum  from  one  hole  center  continuously  overlaps  the  Cu2+-Vij  hole 
center  and  a  second  (low  concentration)  “hole”  center  has  little  angular  dependence.  This 
second  low  concentration  hole  center  was  not  investigated  in  depth  in  this  dissertation. 
After  an  analysis  of  the  copper-related  EPR  spectra,  this  study  continued  with  an 
investigation  of  the  thermal  stability  of  the  copper  defects.  It  is  shown  that  the  x-ray- 
induced  electron  and  hole  centers  recombine  at  temperatures  that  correspond  to 
thermoluminescence  glow  curve  peak  positions  reported  in  the  literature  [68,  70,  72]. 

5.2.  X-Ray-Induced  Electron  Centers  in  Copper-Doped  LTB 

The  spectra  of  the  two  electron  centers  partially  overlap  making  identification 
difficult  but  not  impossible.  These  electron  centers  are  the  result  of  Cu+  interstitial  ions 
trapping  an  electron  during  irradiation.  The  trapped  electron  converts  the  interstitial  Cu+ 
(3d10)  ions  into  Cu°  (Sd’^s1)  atoms.  The  two  x-ray-induced  electron  centers  have  large 
hyperfine  parameters.  The  large  hyperfine  is  due  to  the  interaction  between  the  S  =  1/2 
unpaired  electron  and  the  copper  nucleus.  Calculations  published  by  Morton  and  Preston 
[74-75]  predict  a  large  hyperfine  constant  for  copper  nuclei.  The  electron  center  spectra 
are  described  in  this  section  along  with  support  for  this  defect  model. 

Prior  to  x-ray  irradiation,  the  EPR  spectrum  of  copper-doped  lithium  tetraborate 
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Cu2+liole  centers 


Fig.  5.1.  EPR  spectrum  of  x-ray  irradiated  copper-doped  lithium  tetraborate.  This 
spectrum  was  taken  at  45  K  with  the  magnetic  field  along  the  [001]  direction. 

consists  of  a  set  of  EPR  transitions  in  the  2700  to  2850  Gauss  range  representing  the 
Cu2+-VLi  hole  center  reported  by  Corradi  et  al.  [71]  After  x  ray  exposure,  two  new  EPR 
spectra  are  present.  These  two  spectra  consist  of  four  widely  separated  sets  of  lines  and 
represent  the  two  electron  centers.  These  spectra  with  their  multiple  lines  are  due  to  two 
similar,  yet  different,  Cu°  electron  centers  in  the  LTB  crystal.  The  spectrum  shown  in 
Fig.  5.1  was  taken  after  a  one  hour  exposure  to  x-ray  irradiation  produced  by  a  rotating 
anode  x-ray  tube.  The  Cu2+-VLi  hole  center  and  a  second  hole  center,  discussed  in  more 
detail  later  in  this  chapter,  extend  off-scale  in  the  figure  due  to  the  EPR  spectrometer 
settings  used  to  optimize  the  signal  of  the  electron  centers.  The  two  Cu°  electron  centers 
are  identified  and  labeled  as  Cu°  electron  centers  #1  and  #2.  These  two  centers, 
representing  two  different  defects,  can  be  grouped  into  sets  of  four  transitions  because  of 
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the  two  isotopes  of  copper.  The  stick  diagrams  for  each  center  designate  the  mid-points 
of  the  transition  sets  for  each  isotope.  Copper  has  two  isotopes,  65Cu  and  63Cu,  that  have 
relative  abundances  of  30.8%  and  69.2%,  respectively.  They  possess  similar  magnetic 
moments  with  a  ratio  of  65Cu  to  63Cu  being  1.070  and  both  have  a  nuclear  spin  of  I  =  3/2. 
As  mentioned  earlier,  copper  atoms  with  their  3d !  ()4s 1  configuration  are  expected  to  have 
very  large  hyperfine  interactions.  These  properties  will  be  used  to  explain  the  spectra 
resulting  from  the  Cu°  electron  centers. 

The  energy  level  diagrams  presented  in  this  dissertation,  so  far,  have  been  made 
with  the  assumption  that  the  electron  Zeeman  term  of  the  spin  Hamiltonian  has  a  much 
larger  contribution  than  the  hyperfine  term,  or  any  other  term.  From  this  assumption  the 
effective  g,  or  geJf,  can  be  calculated  by  using  the  mid-point  of  a  set  of  EPR  transitions. 
For  example,  in  Chapter  Three,  a  gejf  could  be  calculated  from  the  center  of  the  seven 
peak  spectrum  presented  in  Fig.  3.2(a).  This  is  not  the  case  for  the  Cu°  electron  centers; 
their  effective  g  value  corresponds  to  a  lower  magnetic  field  value  than  the  mid-point  of 
the  spectrum.  The  spin  Hamiltonian  for  the  copper  electron  centers  are  the  same  as  any 
S  =  1/2  electron  interacting  with  an  I  =  3/2  nucleus,  with  the  exception  that  the  hyperfine 
parameter,  A,  will  have  a  significant  contribution  to  the  energy  levels.  This  behavior  is 
described  in  Appendix  C  of  the  book  “ Electron  Paramagnetic  Resonance ”  by  John  Weil 
and  James  Bolton.  [25]  The  resulting  equations  for  the  energies  of  a  spin  Hamiltonian 
are  dependent  upon  the  magnitude  of  the  applied  magnetic  field.  The  equations  for  these 
energy  levels  are  called  the  ‘Breit-Rabi’  formulas.  A  Breit-Rabi  energy  level  diagram  of 
an  electron  of  spin  S  =  1/2  interacting  with  a  nucleus  of  spin  I  =  3/2  is  illustrated  in 
Fig. 5. 2.  The  Breit-Rabi  formulas  predict  the  changes  in  an  EPR  spectrum  when  the 
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S  I  2. 1  =  3/2  small  hyperfine  Interaction 


g  value  and  middle  of  spectmm 


Fig  5.2.  Breit-Rabi  energy  level  diagrams  illustrating  the  interaction  of  an  S  =  1/2 
electron  spin  interacting  with  an  I  =  3/2  nuclear  spin.  (Left)  Energy  level  diagram 
demonstrating  a  small  hyperfine  interaction.  (Right)  Energy  level  diagram  demonstrating 

a  large  hyperfine  interaction. 


hyperfine  interaction  becomes  progressively  larger.  A  symmetric  EPR  spectrum  is 
produced  when 

the  electron  Zeeman  term  is  considerably  larger  than  the  hyperfine  term.  Conversely,  an 
asymmetric  EPR  spectrum  is  produced  when  the  hyperfine  term  approaches  the  electron 
Zeeman  term.  These  effects  are  shown  in  Fig.  5.2.  The  red  double  headed  arrows 
represent  the  constant  energy  (AE  =  hv)  provided  by  the  microwaves  in  an  EPR 
experiment.  It  may  be  deceptive,  but  the  arrows  are  all  the  same  length.  The  left  portion 
of  Fig.  5.2  illustrates  the  positions  of  EPR  transitions,  and  the  resulting  EPR  spectrum 
that  occurs,  when  there  is  only  a  small  contribution  from  the  hyperfine  term  of  the  spin 
Hamiltonian.  The  right  hand  portion  of  Fig.  5.2  illustrates  the  positions  of  EPR 
transitions,  and  the  resulting  EPR  spectrum  that  occurs,  when  there  is  a  large  contribution 
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from  the  hyperfine  term  of  the  spin  Hamiltonian.  The  large  hyperfine  value  of  copper 
and  the  Breit-Rabi  energy  level  diagram  of  Fig.  5.2  thus  explain  the  relative  spacings  of 
the  sets  of  transitions  in  the  EPR  spectrum  in  Fig.  5.1  for  the  Cu°  electron  centers. 

The  Cu°  electron  centers  are  complex  and  thus  difficult  to  describe.  Each  center 
is  the  result  of  an  unpaired  s  shell  electron  interacting  with  its  copper  nucleus  of  spin 
I  =  3/2.  This  causes  the  spectrum  to  split  into  four  groups  of  spectra  that  have  a  spacing 
described  by  the  Breit-Rabi  formulas.  If  there  were  no  other  interactions  than  this,  a 
resulting  four  peak  spectrum  would  be  seen.  This  is  not  the  case.  Each  of  the  four 
splittings  from  the  copper  interaction  is  further  split  by  an  interaction  with  a  boron 
nucleus.  Boron  has  two  isotopes  (10B  and  1  B)  that  have  different  nuclear  spins  (1  =  3  and 
I  =  3/2,  respectively)  and  different  abundances  (19.8%  and  80.2%,  also  respectively). 
Since  hyperfine-split  signals  associated  with  the  1  'B  nuclei  will  be  more  intense  because 
of  the  higher  abundance  of  this  isotope,  these  1  'B  lines  will  be  used  to  describe  the  four 
groups  of  EPR  spectra  of  the  Cu°  electron  centers.  Each  Cu°  electron  center  “group”  is 
split  by  an  adjacent  interacting  boron  ion,  which  is  called  a  superhyperfine  interaction. 

1  'B  (with  I  =  3/2)  will  cause  a  Cu  EPR  transition  to  separate  into  four  EPR  lines.  As 
mentioned  earlier,  copper  has  two  isotopes  with  similar  magnetic  moments.  65Cu  has  a 
larger  magnetic  moment  and  thus  has  four  wider  spaced  Cu°  electron  center  groups  than 
those  groups  associated  with  63Cu.  With  their  similar  magnetic  moments,  the  groups 
from  the  two  isotopes  overlap. 

A  diagram  illustrating  the  spectral  splittings  of  Cu°  electron  center  #1  that  takes 
both  copper  isotopes  into  account  is  shown  in  Fig.  5.3.  Starting  from  the  top  of  the 
figure,  an  S  =  1/2  spin  system,  is  represented  interacting  with  each  copper  isotope.  Each 
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I  65Cu  Isotope  ■  63Cu  Isotope 

Fig  5.3.  Diagram  illustrating  the  spectral  splittings  of  Cu°  electron  center  #1  taking  both 
isotopes  into  account  and  the  superhyperfine  splitting  from  a  1  *B  nucleus. 

isotope  (both  with  I  =  3/2)  splits  the  transition  into  four  sets  (four  for  each  isotope)  shown 
with  short  brown  dashes  for  65Cu  and  long  green  dashes  for  63Cu.  Each  of  these  sets  is 
shown  interacting  with  a  1  B  nucleus  that  splits  each  set  into  four  transitions,  shown  again 
in  solid  brown  for  65Cu  and  solid  green  for  63Cu.  The  four  transitions  are  then  spaced 
appropriately  for  their  magnetic  moments.  When  two  transitions  are  close  together,  they 
tend  to  “add”  up  to  a  larger  single  transition  in  an  actual  EPR  spectrum.  The  resulting 
experimental  spectra,  shown  in  the  lower  portion  of  the  figure,  represent  the  groups  of 
lines  associated  with  Cu°  electron  center  #1  and  were  taken  from  Fig.  5.1.  The  transition 
labeled  with  an  “x”  is  not  shown  because  it  is  covered  by  the  spectra  from  the  hole 
centers.  The  experimental  EPR  spectra  in  Fig.  5.3  have  more  transitions  than  the  ones 
described  above  and  labeled  in  the  figure.  These  additional  transitions  arise  from  10B 
interactions  as  well  as  the  Cu°  electron  center  #2. 

There  are  several  ways  to  verify  that  the  correct  assignments  have  been  made  for 
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Fig.  5.4.  Expanded  view  of  the  lowest  “group”  of  lines  for  the  two  Cu°  electron  centers. 
The  stick  diagrams  indicate  the  four  1  B  superhyperfine  lines  for  the  two  Cu  isotopes  in 
each  center.  A  single  set  of  10B  superhyperfine  lines  is  shown  for  electron  center  #1  and 

its  associated  63Cu  isotope. 


the  Cu°  electron  centers.  Again,  focus  will  be  given  to  the  lower  group  of  lines  in  the 
1600  to  2200  Gauss  magnetic  field  range  for  the  two  Cu°  electron  centers,  shown  in  Fig. 
5.4.  The  two  peaks  on  the  far  right  hand  side  of  the  figure  at  2096  and  2120  Gauss  have 
intensities  (taken  from  the  experimental  spectrum)  of  2.53  and  17.34  units,  respectively, 
in  height.  These  two  peaks  represent  10B  and  1  'B  ion  interactions  for  the  63Cu  isotope. 
Multiplying  17.34  units  by  4  ('  'B  with  I  =  3/2  results  in  four  spectral  peaks)  equals  69.36 
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units.  This  represents  the  total  magnitude  of  transitions  within  this  set  that  are  due  to  1  'B. 
Next,  multiplying  2.53  units  by  7  (l0B  with  1  =  3  results  in  seven  spectral  peaks)  equals 
17.71  units  and  represents  the  total  magnitude  of  transitions  within  this  set  that  are  due  to 
l0B.  Dividing  the  10B  natural  abundance  (19.8%)  by  the  1 1 B  natural  abundance  (80.2%) 
gives  a  ratio  of  0.2469.  Finally,  dividing  17.71  units  by  69.36  units  gives  a  ratio  of 
0.2553  and  is  very  close  to  boron’s  natural  abundance  ratio.  Carrying  out  a  similar 
comparison  for  a  single  63Cu  peak  intensity  to  that  of  a  single  65Cu  peak  intensity  yielded 
a  ratio  of  0.4266  that  is  also  in  good  agreement  with  the  0.445 1  ratio  of  the  two  natural 
abundances  for  copper. 

Recall  from  Chapter  Two  that  EPR  transitions  are  separated  by  their  hyperfine 
parameters  and  that  those  hyperfine  values  are  also  directly  proportional  to  their  magnetic 
moments.  This  allows  for  a  comparison  of  the  spacing  between  the  seven  10B  transitions 
and  the  four  1  B  hyperfine  transitions  associated  with  the  63Cu°  portion  of  electron  center 
#1.  The  spacing  between  adjacent  lines,  or  the  A  value,  for  the  1 1 B  interaction  is 
measured  to  be  52.33  Gauss.  Likewise,  the  spacing  between  adjacent  lines,  or  the  A 
value,  for  the  l0B  hyperfine  interaction  is  measured  to  be  17.67  Gauss.  The  ratio  of  1()B 
hyperfine  (or  A  value)  to  1  'B  hyperfine  is  0.3376  and  is  in  excellent  agreement  with  the 
ratio  of  the  two  magnetic  moments  of  boron,  or  0.3347.  Comparing  the  hyperfine  in  this 
manner  works  as  long  as  the  hyperfine  parameter  does  not  significantly  affect  the  overall 
energy  of  the  spin  Hamiltonian.  The  hyperfine  of  copper  has  a  significant  effect  on  the 
overall  energy  of  the  spin  Hamiltonian. 

Because  the  Cu°  electron  centers  are  isotropic,  data  taken  with  the  magnetic  field 
aligned  parallel  to  the  c  axis  of  the  crystal  was  used  to  determine  the  g  parameters.  A 
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least-squares  fitting  routine,  similar  to  those  used  in  earlier  chapters  of  this  dissertation, 
was  used  to  solve  for  the  g  value.  Only  one  set  of  data  (magnetic  field  values  and 
microwave  frequencies)  is  needed.  Since  the  defect  model  places  the  electron  in  an  s 
orbital,  the  experimentally  measured  g  should  be  very  close  to  ge,  the  g  value  of  a  free 
electron.  Indeed,  this  is  the  case,  since  both  electron  centers  (#1  and  #2)  have  a  g  value 
of  2.0023  ±  0.0005,  as  determined  in  the  least-squares  fitting  routine.  The  error  was 
estimated  by  using  the  line  widths  of  the  EPR  spectra. 

At  the  same  time,  the  hyperfine  parameters  were  determined  for  the  electron 
centers.  This  is  not  a  straight  forward  calculation  as  in  earlier  chapters  due  to  the  uneven 
spacings  of  the  EPR  transitions  that  are  due  to  the  large  hyperfine  interaction.  John  Weil 
has  a  paper  titled  “The  Analysis  of  Large  Hyperfine  Splitting  in  Paramagnetic  Resonance 
Spectroscopy”  [73].  This  paper  starts  with  an  isotropic  spin  Hamiltonian  and  can  be  used 
with  hyperfine  values  that  are  not  “too  large.”  In  his  derivation,  Weil  arrives  at  the 
following  equation: 


^-mj  ^mj 

2  irq 


5.1 


In  the  above  equation,  B_mi  —  Bnl[,  is  the  magnetic  field  position  corresponding  to  a 
specific  nuclear  transition  (mi).  It  may  help  to  look  back  at  the  four  EPR  transitions 
shown  for  the  right  side  of  Fig.  5.2.  These  four  transitions  go  with  specific  values  of  mi. 
Also  in  Eq.  5.1,  Aiso  represents  the  isotropic  hyperfine  parameter,  and  giso  is  the  isotropic 
value  calculated  in  the  preceding  paragraph.  Using  the  magnetic  field  positions  of  the 
midpoints  of  the  sets  of  transitions  for  each  Cu°  center  and  each  isotope  of  copper,  Aiso 
could  be  determined  for  each  isotope.  For  clarity,  the  magnetic  field  positions  that  were 
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used  to  solve  for  the  hyperfine  values  of  each  isotope  are  shown  in  Fig.  5.1  (i.e.,  the 
vertical  lines  in  the  stick  diagrams  above  and  below  the  spectra  for  electron  centers  #1 
and  #2).  The  values  of  the  hyperfine  parameters  for  the  two  electron  centers  are  listed  in 
Table  5.1. 

Table  5.1.  Isotropic  spin-Hamiltonian  parameters  for  Cu°  electron  centers  #1  and  #2. 
These  values  were  obtained  from  the  EPR  spectra  in  Fig.  5.1.  Estimated  error  limits  are 
±  0.0005  for  the  g  values  and  ±15  MHz  for  the  A  values. 

Defect  g  A(63Cu)  A(65Cu) 

Cu°  #1  2.0023  2480  MHz  2663  MHz 

Cu°  #2  2.0023  2178  MHz  2333  MHz 


The  two  defects  are  very  similar  and  defect  models  to  support  the  differences  in 
the  centers  would  be  difficult  to  establish.  The  superhyperfine  interaction  with  boron 
strongly  suggests  a  model  in  which  these  electron  centers  are  interstitial  copper  atoms. 
One  of  these  two  defects  could  be  the  result  of  an  interstitial  copper  that  is  isolated  from 
any  other  defects,  while  the  second  defect  is  the  result  of  an  interstitial  copper  with  a 
nearby  defect  that  “slightly”  changes  the  hyperfine  parameter.  The  g  parameter  is  the 
same  for  both  centers  because  the  unpaired  spin  is  in  an  s  orbital  and  thus  has  no  orbital 
angular  momentum.  There  is  a  table  provided  by  Morton  and  Preston  [74-75]  that 
predicts  the  isotropic  hyperfine  constant  for  each  element  in  the  periodic  table.  Their 
prediction  for  65Cu  is  6410  MHz.  The  average  hyperfine  for  the  two  electron  centers 
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containing  65Cu  is  2498  MHZ.  This  means  that  39%  of  the  unpaired  spin  is  in  a  4s  orbital 
and  the  remaining  spin  density  resides  elsewhere,  such  as  on  the  neighboring  boron 
nucleus,  as  illustrated  in  Fig.  5.1  where  the  spin  of  the  electron  is  split  by  an  interacting 
boron  nucleus.  This  large  spin  density  residing  on  the  boron  ion  is  further  evidence  to 
support  the  model  that  these  two  electron  centers  are  interstitial  copper  atoms.  If  a 
copper  atom  were  occupying  a  lithium  site,  a  smaller  interaction  with  the  neighboring 
substitutional  boron  would  be  expected  because  of  the  large  separation  distance  of  about 

o 

2.7  A  from  a  lithium  site  to  the  nearest  boron  lattice  site. 

5.3.  X-Ray-Induced  Hole  Centers  in  Copper-Doped  LTB 

As-grown  copper-doped  lithium  tetraborate  crystals  contain  a  large  signal  due  to 
Cu2+  ions  that  has  been  reported  in  the  literature,  labeled  Cu2+-Vij.  This  was  the  only 
EPR  signal  observed  in  as-grown  crystals.  After  exposure  to  x-rays,  a  radiation-induced 
Cu2+  hole  center  signal  is  present  and  its  EPR  spectrum  overlaps  the  spectrum  assigned  to 
the  Cu2+-Vu  center.  Significant  overlap  occurs  for  every  orientation  of  magnetic  field. 
This  radiation-induced  Cu2+  hole  center,  henceforth  referred  to  as  the  Cu2+-active  center, 
is  best  observed  when  the  magnetic  field  is  oriented  parallel  to  the  crystal  [110]  direction. 
Recall  that  the  basic  B4O96  unit  in  lithium  tetraborate  has  eight  crystallographic 
orientations  in  the  LTB  crystal  lattice  and  a  paramagnetic  defect  would  have  an 
equivalent  set  of  eight  crystallographic  orientations.  When  the  magnetic  field  is  aligned 
parallel  to  the  c  axis  of  the  crystal  all  eight  defect  sites  are  magnetically  equivalent  and 
produce  a  single  EPR  spectrum.  However,  when  the  magnetic  field  is  aligned  parallel  to 


97 


Cu2+-V 


Li 


Fig.  5.5.  Upper  EPR  spectrum  showing  the  Cu2+-VLi  characterized  by  Corradi  et  al.  [71]. 
Lower  EPR  spectrum  showing  both  the  Cu2+-Vij  center  and  the  radiation  induced  Cu2+- 
active  center  after  irradiation.  These  spectra  were  taken  at  25  K  with  the  magnetic  field 

parallel  to  the  crystal  [110]  direction. 
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the  a  axis  or  the  crystal  [110]  direction,  the  eight  sites  separate  into  two  magnetically 
inequivalent  sets,  each  with  fourfold  degeneracy.  The  two  magnetically  inequivalent  sets 
are  identified  and  labeled  in  Fig.  5.5  for  both  Cu2+-Vij  and  Cu2+-active  centers. 

The  EPR  spectra  of  the  Cu2+  hole  centers,  with  their  sets  of  four  equally  spaced 
lines,  are  the  result  of  an  S  =  1/2  spin  interacting  with  a  copper  (I  =  3/2)  nucleus.  This 
interaction  results  in  a  four-peak  spectrum  that  is  best  seen  in  the  upper  portion  of  Fig. 

5.5  for  the  four  peaks  between  2800  and  3000  Gauss.  Since  65Cu  has  a  larger  magnetic 
moment  than  63Cu,  the  outer  lines  associated  with  the  different  isotopes  within  a  set  of 
four  lines  often  do  not  overlap  at  higher  fields.  This  is  best  seen  in  the  upper  portion  of 
Fig.  5.5  for  the  peaks  between  3100  and  3500  Gauss.  The  outer  most  peaks  in  this  upper 
set  are  due  solely  to  65 Cu  and  the  peaks  next  to  each  of  these  are  due  solely  to  63Cu.  The 
middle  two  peaks  in  this  set  are  overlapping  peaks  due  to  the  presence  of  both  copper 
isotopes  and  this  explains  their  misshapen  nature,  i.e.,  there  are  actually  two  overlapping 
lines  at  each  location.  In  the  lower  set  of  peaks  between  2800  and  3000  Gauss,  the  lines 
due  to  both  isotopes  are  all  overlapping.  This  overlap  results  in  broadened  peaks  and  the 
individual  isotopes  are  not  resolved.  The  lower  portion  of  Fig.  5.5  contains  the  lines  of 
the  Cu2+-Vu  center  plus  additional  lines  from  the  Cu2+-active  center.  This  “active”  center 
has  been  given  this  name  because  it  appears  during  irradiation  and  disappears  after 
irradiation  when  the  temperature  is  increased  (when  this  hole  center  recombines  with  the 
electron  centers  #1  and  #2).  The  EPR  signal  intensity  (i.e.,  concentration)  of  the  Cu2+-Vij 
center  does  not  change  during  an  annealing  process  and  thus  do  not  “actively”  participate 
in  electron-hole  recombination.  This  will  be  discussed  further  in  the  thermal  stability 
section  of  this  chapter. 
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Fig.  5.6.  EPR  spectra  of  the  Cu^-Vu  and  Cu2+-active  centers.  Spectra  taken  before  and 
after  irradiation  are  compared.  These  data  were  taken  at  25  K.  The  non-irradiated 
spectra  contain  the  Cu2+-Vu  center  alone.  The  Cu2+- active  center  is  present  in  the 
irradiated  spectra,  but  difficult  to  identify. 
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A  set  of  spin-Hamiltonian  parameters  describing  the  Cu2+-active  center  was  not 
obtained  in  the  current  study  due  to  the  strong  overlap  with  the  Cu2+-Vij  center.  The  EPR 
spectrum  of  the  “active”  center  is  best  observed  when  the  magnetic  field  is  oriented 
parallel  to  the  crystal  [110]  direction.  The  overlap  between  the  two  spectra  is  even  more 
extreme  for  other  high  symmetry  directions.  This  overlapping  nature  is  illustrated  in  Fig. 
5.6,  where  spectra  are  shown  for  the  magnetic  field  oriented  along  the  crystal  a  and  c 
directions.  In  both  cases,  a  full  set  of  Cu2+-active  center  peaks  cannot  be  identified. 
Though  not  shown,  the  spectra  were  observed  at  many  orientations  of  magnetic  field 
relative  to  the  crystal,  and  overlapped  in  all  cases.  This  leads  to  the  conclusion  that  this 
radiation  induced  Cu2+-active  center  has  a  set  of  spin  Hamiltonian  parameters  that  are 
close  to  those  published  by  Corradi  et  al.  [71]  for  the  Cu2+-Vij  center,  and  thus  is  a 
similar  defect.  The  Cu2+-active  center  is  likely  a  Cu2+  ion  sitting  at  a  Li+  substitutional 
site  that  does  not  have  a  neighboring  Li+  vacancy  stabilizing  the  hole  as  put  forward  for 
the  Cu2+-VLi  center. 

There  are  two  additional  items  concerning  hole  centers  in  copper-doped  lithium 
tetraborate  that  need  elaboration.  First,  as  shown  in  Fig.  5.1,  there  is  a  less  intense  hole 
center  located  between  2950  and  3300  Gauss  that  is  produced  by  x-ray  radiation.  It  is 
clearly  associated  with  copper  as  the  outer  lines  show  the  typical  copper  63Cu/65Cu 
splitting.  However,  its  concentration  is  much  less  than  that  of  the  Cu2+-Vij  and  Cu2+- 
active  centers  and  is  not  investigated  in  any  detail.  This  “small”  Cu-related  center  did  not 
appear  to  contribute  additional  information  to  the  present  study.  Second,  while 
performing  the  EPR  experiments  on  copper-doped  lithium  tetraborate  and  viewing  the 
spectra  along  many  magnetic  field  orientations,  it  was  observed  that  the  hyperfine 


101 


splittings  for  the  Cu2+-Vu  center  reduced  to  zero  for  a  few  orientations  of  magnetic  field. 
Corradi  et  al.  [71]  published  hyperfine  parameters  of-16.1,  -57.5,  and  97.6  MHz  for  the 
Cu2+-Vu  center.  The  observation  that  the  hyperfine  decreases  to  zero  at  some  angles 
provides  evidence  that  the  hyperfine  parameters  do  not  all  have  the  same  sign.  When 
hyperfine  parameters  are  either  all  positive  or  all  negative,  the  hyperfine  splitting  cannot 
reduce  to  zero.  This  is  the  case  for  the  hyperfine  in  silver-doped  lithium  tetraborate 
crystals  where  there  is  always  an  observable  hyperfine  splitting  and  it  has  been  found  that 
the  hyperfine  parameters  are  all  positive  in  sign. 

5.4.  Thermal  Stability  of  the  Radiation-Induced  Centers 

After  identification  of  the  radiation-induced  electron  and  hole  centers,  their 
thermal  stabilities  were  determined.  The  copper-doped  lithium  tetraborate  crystals  were 
subject  to  a  pulsed  isochronal  anneal.  Crystals  were  first  exposed  for  one  hour  with  x 
rays  generated  by  a  rotating  anode  x-ray  tube.  The  subsequent  anneal  was  accomplished 
by  setting  an  oven  to  a  set  temperature,  placing  the  sample  in  a  ceramic  boat,  and  sliding 
the  boat  into  the  center  of  the  oven  within  a  quartz  tube  that  runs  through  the  oven.  The 
sample  remained  in  the  center  of  the  oven  for  two  minutes  at  each  temperature.  After 
heating,  the  sample  was  quickly  removed  from  the  oven.  Electron  paramagnetic 
resonance  (at  25  K  for  the  hole  center  and  45  K  for  the  electron  centers)  was  used  to 
determine  the  defect  concentration  after  each  annealing  step.  The  defect  concentration 
was  determined  by  measuring  the  intensities  and  line  widths  of  EPR  peaks  that 
corresponded  to  each  defect  center.  The  monitored  defect  centers  included  the  two 
electron  centers,  the  Cu2+-Vu  center,  and  the  Cu2+-active  center.  The  EPR  isochronal 
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pulsed  anneal  and  the  thermoluminescence  experiments  (discussed  in  the  next  paragraph) 
were  repeated  multiple  times  to  ensure  reproducibility. 

Thermoluminescence  (TL)  data  were  obtained  using  a  Harshaw  TLD  reader 
where  the  resulting  glow  curves  represent  total  light  output  versus  temperature.  The 
heating  rate  was  1  °C/s;  this  is  the  slowest  heating  rate  available  on  the  instrument.  The 
thermoluminescent  glow  curve  for  the  copper-doped  LTB  crystals  showed  two  peaks,  one 
near  100  °C  and  the  other  near  200  °C.  Thermoluminescent  glow  curves  containing  two 
peaks  have  been  reported  for  copper-doped  lithium  tetraborate  in  the  literature  [68,  70, 

72]  and  some  authors  have  reported  similar  two  peak  TL  glow  curves  for  manganese- 
doped  LTB  [76-78].  Only  one  above-room-temperature  TL  peak  has  been  reported  in 
silver-doped  LTB  [51]. 

The  results  of  the  EPR  isochronal  anneal  are  presented  in  Fig.  5.7  and  are  overlaid 
with  the  TL  glow  curve.  In  the  Cu2+-active  center  decay  data,  two  decreasing  steps  are 
observed  near  100  and  200  °C,  respectively.  The  disappearance  of  this  hole  appears  to  be 
related  to  the  TL  glow  peaks.  The  electron  center  is  different.  During  the  isochronal 
annealing  experiment,  the  Cu°-electron  centers  increase  in  intensity  between  25  and  125 
°C  and  then  decay  entirely  between  125  and  225  °C.  These  data,  which  were  collected  on 
several  occasions,  suggests  that  there  must  be  an  unidentified  electron  center  that 
participates  in  recombination  between  25  and  125  °C,  at  temperatures  lower  than  the 
decay  of  the  Cu°-electron  centers.  As  this  unidentified  center  releases  electrons,  some 
recombine  with  the  Cu2+-active  center  causing  its  decay  and  others  become  trapped  at 
unpopulated  Cu°-electron  centers  thus  increasing  their  signal  intensity.  In  previous  work 
[31],  LTB  has  been  shown  to  have  significant  concentrations  of  lithium  and  oxygen 
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vacancies.  In  this  earlier  work,  the  oxygen  vacancies  trapped  electrons  and  could  be  a 


possible  candidate  for  the  unidentified  electron  centers. 


Fig.  5.7.  Thermal  stabilities  of  the  Cu2+-active  hole  center  (blue)  and  the  two  Cu° 
electron  centers  (red),  as  determined  from  the  magnitudes  of  their  EPR  spectra.  The  data 
points  were  obtained  by  collecting  EPR  signal  intensities  at  25  and  45  K,  respectively. 
The  effective  heating  rate  during  the  EPR  pulsed  anneal  experiment  was  0.2  °C/s.  The 
thermoluminescence  glow  curve  (black)  was  obtained  by  collecting  the  total  light  output 

with  a  heating  rate  of  1  °C/s. 
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5.5.  Discussion 


In  the  present  study,  two  new  trapped  electron  centers  and  two  new  trapped  hole  centers 
have  been  identified  after  exposing  copper-doped  lithium  tetraborate  crystals  to  x-ray 
radiation.  The  electron  centers  are  isotropic  and  have  similar  hyperfine  parameters.  The 
radiation-induced  hole  center  has  g  and  hyperfine  parameters  that  are  similar  to  a 
previously  identified  hole  center  [71].  These  radiation-induced  centers  are  stable  at  room 
temperature  and  release  their  trapped  electrons  and  holes  as  they  are  warmed  above  room 
temperature.  The  two  peaks  seen  in  thermoluminescent  glow  curves  correspond  to  the 
two  decay  steps  of  the  radiation-induced  hole  center.  This  research  has  identified  the 
defects  responsible  for  the  TL  glow  curves  and  can  be  beneficial  in  utilizing  copper- 
doped  lithium  tetraborate  as  an  x-ray  radiation  dosimeter.  Future  work  could  be  done  to 
determine  the  sensitivity  of  these  crystals  to  lower  magnitude  x-ray  exposures  (i.e., 
medical  x-ray  dose  levels)  furthering  their  potential  in  dosimetry  applications. 
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VI.  Optical  Properties  of  Silver-Doped  Lithium  Tetraborate  Crystals 


6.1.  Optical  Experiment 

Studies  concerning  the  optical  properties  of  lithium  tetraborate  crystals  and 
glasses  can  be  found  in  the  literature  [79-83].  Many  of  these  studies  contain  absorption, 
excitation,  and  emission  spectra  and  offer  theoretical  explanations.  In  this  chapter, 
electron  paramagnetic  resonance  (EPR),  optical  absorption  (OA),  photoluminescence 
(PL),  and  photoluminescence  excitation  (PLE)  data  have  been  collected  in  an  effort  to 
provide  a  more  convincing  argument  for  the  optical  features  observed  from  silver-doped 
lithium  tetraborate  crystals. 

Ag-doped  LTB  crystals  are  transparent  prior  to  being  exposed  to  x-ray  radiation. 
During  exposure,  the  crystals  turn  green  in  color.  Crystals  will  maintain  this  coloration 
for  about  a  month,  if  kept  at  room  temperature.  Electron  paramagnetic  resonance  spectra 
are  not  present  in  colorless  crystals.  This  lack  of  spectra  occurs  in  silver  doped  crystals 
because  all  of  the  silver  impurities  are  in  the  form  of  Ag+  (3d10)  ions  in  as-grown  crystals 
(including  both  interstitial  and  substitutional  Ag+  ions).  In  the  singly  ionized  charge 
state,  there  is  no  unpaired  spin  residing  at  the  silver  ion  and  thus  no  signal  in  an  EPR 
spectrum.  This  behavior  has  been  reported  by  Brant  et  al.  [51]  During  exposure  to  x-ray 
radiation,  the  interstitial  Ag+  ions  trap  electrons  to  form  Ag°  ( 3d 1  ()4s 1 )  atoms  and  the 
substitutional  Ag+  ions  trap  holes  to  form  Ag2+  (3d9)  ions.  Both  of  these  charge  states  are 
paramagnetic  and  results  in  EPR  spectra  with  very  large  intensities.  EPR  spectra 
identifying  the  radiation-induced  trapped  electron  and  trapped  hole  centers  are  shown  in 
Fig.  6.1. 
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An  optical  absorption  experiment  was  conducted  using  a  Varian-Cary 
spectrophotometer.  The  results,  obtained  at  room  temperature,  are  shown  in  Fig.  6.2. 
The  red  curve  shows  the  absorption  spectrum  through  the  visible  and  near- ultraviolet 
region  for  an  as-grown  LTB  crystal. 


Unperturbed  Hole  Center 
10£kg  rn 
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Fig.  6.1.  EPR  Spectrum  illustrating  the  x-ray  radiation  induced  electron  and  hole  centers. 

It  contains  one  band  peaking  near  205  nm,  which  has  been  assigned  in  the  literature  to 
Ag+  ions  [80].  The  black  curve  in  Fig.  6.2  was  taken  after  the  crystal  was  exposed  to  x 
rays  at  room  temperature.  Several  large  optical  absorption  peaks  were  produced  as  a 
result  of  the  x  rays.  The  dominant  peak,  at  390  nm,  has  been  assigned  to  Ag°  atoms  [82- 
83].  Additional  radiation-induced  peaks  occur  at  261,  298,  and  650  nm.  Interestingly, 
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the  peak  at  204  nm  decreased  in  intensity  as  a  result  of  the  x  rays.  This  is  in  agreement 
with  EPR  data  where  Ag+  ions  are  converted  to  Ag°  and  Ag2+  centers. 

Fluorolog  experiments  (PL  and  PLE)  were  conducted  on  pre-  and  post-  x-ray 
irradiated  Ag-doped  LTB  crystals.  Prior  to  irradiation,  a  weak  emission  peak  is  observed 
near  400  nm  when  using  298  nm  light  as  an  excitation  source.  After  irradiation,  a  much 
larger  (more  than  20  times  greater)  emission  peak  is  observed  at  471  nm  when  using  298 
nm  light  as  an  excitation  source.  This  latter  peak  must  be  due  to  a  defect  that  was 


Fig.  6.2.  Optical  absorption  of  pre-  (red)  and  post-  (black)  x-ray  irradiated  LTB:Ag 


crytals. 


produced  by  the  radiation.  The  two  candidates  for  this  defect  are  Ag°  atoms  or  Ag2+  ions 
(see  the  EPR  spectrum  in  Fig.  6.1).  To  determine  which  defect  is  responsible  for  the 
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intense  emission  at  471  nm,  a  photoluminescence  excitation  experiment  was  performed. 
Excitation  spectra  were  collected  on  pre-  and  post-x-ray  irradiated  crystals.  Prior  to 
irradiation,  there  is  a  small  excitation  peak  near  340  nm  when  the  emission 
monochromator  is  set  to  monitor  47 1  nm.  After  the  crystals  were  exposed  to  x-ray 
radiation,  a  very  large  (more  than  20  times  greater)  excitation  peak  is  observed  at  298  nm 
when  the  emission  monochromator  is  set  to  monitor  471  nm.  This  clearly  shows  that  the 
298  nm  excitation  peak  correlates  with  the  471  nm  emission  peak.  The  pre-  and  post¬ 
irradiation  emission  spectra  are  shown  in  Fig.  6.3,  and  the  pre-  and  post-irradiation 
excitation  spectra  are  shown  in  Fig.  6.4.  Based  on  the  PFE  data,  the  471  nm  emission 
peak  can  now  be  conclusively  assigned  to  Ag2+  ions.  The  Ag°  atom  is  eliminated  as  a 
possibility  because  its  broad  absorption  band  at  390  nm  does  not  give  rise  to  the  471  nm 
emission  band.  In  Fig.  6.5,  the  optical  absorption  data  from  the  Varian-Cary 
spectrophotometer  and  the  PFE  data  from  the  Fluorolog  are  plotted  together. 

A  thermal  anneal  study  was  performed  for  the  radiation  induced  optical 
absorption  spectra.  These  results  are  shown  in  Fig.  6.6.  The  isochronal  annealing 
procedure  was  the  same  as  the  one  used  for  EPR  spectra  in  Chapter  Five.  Each  step  in 
the  anneal  consisted  of  holding  the  sample  at  a  fixed  temperature  in  the  furnace  for  two 
minutes.  The  sample  was  then  cooled  and  placed  in  the  spectrophotometer.  The  decrease 
in  the  intensity  as  a  function  of  the  annealing  temperature  correlates  with  the  EPR  data 
previously  published  by  Brant  et  al.  [51] 
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Fig.  6.3.  Emission  (PL)  spectra  of  pre-  and  post-  x-ray  irradiated  crystals. 


Fig.  6.4.  Excitation  (PLE)  spectra  of  pre-  and  post-  x-ray  irradiated  crystals. 
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Fig.  6.5.  Absorption  spectra  (black  and  red)  combined  with  PLE  spectra  (blue). 


Fig.  6.6.  Absorption  spectra  taken  during  isochronal  annealing  experiment.  The  inset  is 

a  magnification  of  the  peaks  at  205  nm. 


Ill 


6.2.  Discussion 


The  research  of  this  chapter  has  demonstrated  that  EPR  spectra  along  with  optical 
absorption  and  fluorospectroscopy  data  (i.e.,  PL  and  PLE)  can  be  combined  to 
conclusively  show  that  the  absorption  peak  near  298  nm  is  due  to  Ag2+  ions  in  x-ray 
irradiated  Ag-doped  LTB  crystals.  This  study  provides  insight  about  the  absorption 
spectrum  of  Ag-doped  LTB  as  well  as  further  confirmation  for  the  defect  models 
presented  in  Chapter  Four  concerning  the  silver  valence  states  both  prior  to  and  after 
exposure  to  x-ray  radiation.  Understanding  the  photoluminescence  and  absorption 
characteristics  of  LTB  are  also  important  in  utilizing  LTB  crystals  for  dosimetry 
applications. 

Absorption  spectra  are  commonly  plotted  as  energy  versus  optical  density  which 

he 

is  a  direct  conversion  using  E  =  — .  A  slight  correction  for  wavelength  dependence  to 

A 

spectral  response  is  also  needed  in  the  emission  and  excitation  data  presented.  This 
correction  is  illustrated  in  appendix  A  of  Luminescent  Materials  [84]. 
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VII.  Conclusions 


7.1.  Summary  of  Findings 

Lithium  tetraborate  (LLB4O7),  commonly  referred  to  as  LTB,  is  a  versatile  insulating 
crystal  with  potential  use  in  neutron  detection.  LTB  has  been  grown  using  enriched  6Li 
and  10B  to  increase  its  neutron  opacity  [4].  When  doped  with  an  appropriate  optically 
active  element,  LTB  could  be  used  as  an  effective  scintillation  material  and  thus  is  an 
important  candidate  for  high- sensitivity  thermoluminescence  dosimetry  applications  [52- 
54].  Radiative  processes  associated  with  these  dopants  are  highly  efficient  in  this 
material,  and  bright,  above-room-temperature  thermoluminescence  peaks  are  observed. 
Lithium  tetraborate  has  been  grown  pure  and  doped  with  many  different  elements 
including  transition  metals,  actinides,  and  rare  earth  elements.  The  research  in  this 
dissertation  addressed  undoped  lithium  tetraborate  as  well  as  LTB  doped  with  copper  and 
silver.  Specifically,  this  dissertation  has  focused  on  identifying  and  characterizing  point 
defects  in  lithium  tetraborate  crystals  as  well  as  explaining  their  role  in  electron  and  hole 
recombination  processes.  These  electron-hole  recombination  processes  have  been  shown 
to  explain  the  thennoluminescent  glow  curves  in  both  Cu-  and  Ag-doped  LTB  crystals. 

In  Chapter  Three,  the  research  has  provided  insight  to  the  use  of  lithium 
tetraborate  as  a  possible  neutron  detection  material.  Experiment  has  shown  that  there  are 
large  numbers  of  lithium  and  compensating  oxygen  vacancies  in  the  as-grown  crystals 
[31].  These  vacancies  are  responsible  for  trapping  electrons  and  holes.  The  trapped  hole 
center  investigated  in  this  chapter  has  a  very  low  thermal  stability,  less  than  90  K,  and  is 
populated  during  exposure  to  ionizing  radiation.  These  trapped  holes  recombine  with 
trapped  electrons  and  release  visible  light  when  the  crystal  is  heated  above  90  K.  This  is 
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demonstrated  with  the  low  temperature  thermoluminescent  data  presented.  A  full  set  of 
spin-Hamiltonian  parameters  were  determined  for  the  low  temperature  hole  center.  The 
hole  trapping  centers  may  have  different  effects  on  the  performance  of  LTB  as  a 
scintillating  material.  A  nuclear  event,  i.e.,  neutron  capture,  in  LTB  is  accompanied  by 
ionizing  radiation.  This  ionizing  radiation  produces  electrons  and  holes  that  will  then 
migrate  through  the  crystal  lattice.  These  electrons  and  holes  may  recombine  at  the  hole 
trapping  site  investigated  in  Chapter  Three  and  produce  visible  light  that  could  be 
collected  in  a  photomultiplier  tube.  This  would  be  a  benefitting  defect  in  neutron 
detection.  However,  ionizing  radiation  unrelated  to  neutron  capture  may  also  produce 
electron  hole  pairs  that  again  recombine  at  the  hole  trapping  site  investigated  and  lead  to 
a  “false”  count  in  a  neutron  detector.  This  would  be  a  negative  effect  in  determining 
neutron  capture. 

Concerning  Chapter  Four,  there  are  several  significant  achievements.  A  full  set  of 
spin  Hamiltonian  parameters  were  determined  for  a  defect  that  had  been  previously 
reported  [31]  but  not  characterized.  The  method  of  characterization,  i.e.,  using  EPR 
alone  to  solve  for  the  g  and  hyperfine  parameters,  was  validated  by  comparing  the 
research  of  Chapter  Four  to  published  results  for  the  unperturbed  hole  center.  This 
chapter  has  shown  that  it  is  possible  to  find  a  good  set  of  g  and  hyperfine  parameters 
using  EPR  alone  for  a  defect  that  contains  two  isotopes  of  similar  nuclear  magnetic 
moments  that  could  not  otherwise  have  been  characterized.  Now,  all  of  the  defects  have 
been  characterized  that  participate  in  electron-hole  recombination  and  contribute  to  the 
thermoluminescent  glow  curve  of  x-ray  exposed  Ag-doped  LTB.  In  order  for  lithium 
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tetraborate  to  become  a  viable  candidate  for  dosimetry  and  nuclear  detection  applications 
[55-56],  it  is  important  to  characterize  all  of  the  defects  in  the  material. 

While  performing  research  for  Chapter  Five,  an  attempt  to  increase  signal 
intensity  of  a  previously  identified  trapped  hole  center  in  Cu-doped  lithium  tetraborate 
led  to  the  discovery  of  several  new  spectra  that  have  not  been  reported  in  the  literature. 
Two  new  similar,  yet  different,  electron  centers  are  identified  and  characterized.  Two 
new  hole  centers  are  also  identified.  After  establishing  models  for  the  copper-related 
EPR  spectra,  the  study  continued  with  an  investigation  of  the  thermal  stability  of  the 
copper  defects.  It  is  shown  that  the  x-ray-induced  electron  and  hole  centers  recombine  at 
temperatures  that  correspond  to  thermoluminescence  peak  temperatures  reported  in  the 
literature  [68,  70,  72].  This  research  has  identified  the  defects  responsible  for  the  TL 
glow  curves  and  can  be  beneficial  in  utilizing  copper-doped  lithium  tetraborate  as  an  x- 
ray  radiation  dosimeter. 

In  Chapter  Six,  electron  paramagnetic  resonance  (EPR),  optical  absorption  (OA), 
photoluminescence  (PL),  and  photoluminescence  excitation  (PLE)  data  have  been 
collected  in  an  effort  to  provide  a  comprehensive  description  of  the  optical  features 
observed  from  silver-doped  lithium  tetraborate  crystals.  The  results  in  this  chapter  show 
that  EPR  spectra  along  with  optical  absorption  and  fluorospectroscopy  data  (i.e.,  PL  and 
PLE)  can  be  combined  to  conclusively  establish  that  the  absorption  peak  near  298  nm  is 
due  to  Ag2+  ions  in  x-ray  irradiated  Ag-doped  LTB  crystals.  This  study  provides  insight 
about  the  absorption  spectrum  of  Ag-doped  LTB  as  well  as  further  confirmation  for  the 
defect  models  presented  in  Chapter  Four  concerning  the  silver  valence  states  both  prior  to 
and  after  exposure  to  x-ray  radiation.  Understanding  the  photoluminescence  and 
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absorption  characteristics  of  LTB  are  also  important  in  utilizing  LTB  crystals  for 
dosimetry  applications. 

7.2.  Future  Work 

If  LTB  is  to  be  incorporated  in  neutron  detectors  and  sensors,  much  research 
remains  to  be  done.  This  dissertation  has  characterized  the  intrinsic  hole  trap  in  LTB,  but 
has  not  addressed  the  intrinsic  electron  trap  (i.e.,  the  oxygen  vacancies).  The 
paramagnetic  defects  in  silver-doped  LTB  have  been  investigated  at  length  and  the  spin- 
Hamiltonian  parameters  have  been  determined  for  both  electron  and  hole  centers.  Silver- 
doped  LTB  has  also  been  optically  characterized.  Initial  work  has  been  performed  to 
optically  characterize  copper-doped  LTB  but  work  remains  and  the  preliminary  results 
were  not  included  in  this  dissertation. 

Additional  investigations  should  be  considered  if  a  wider  variety  of  samples  were 
available  that  complemented  those  used  in  the  present  work.  Due  to  their  light-emitting 
properties,  copper-  and  silver-doped  LTB  are  strong  candidates  for  scintillation  and 
neutron  detection  applications.  Similar  research  should  be  performed  on  other  dopants 
that  are  known  to  be  optically  active.  On  the  basis  of  size  and  charge,  aluminum  and 
silicon  will  substitute  for  boron  in  the  LTB  lattice.  Nitrogen  may  substitute  for  oxygen. 

If  crystals  were  grown  with  these  dopants,  surprisingly  new  behaviors  may  occur  and 
their  paramagnetic  defect  and  optical  properties  should  be  investigated  in  a  similar 
manner  to  those  experiments  described  in  this  dissertation. 

A  final  interesting  potential  project  relates  to  the  work  performed  in  Chapter 
Three  concerning  undoped  lithium  tetraborate.  The  work  of  Chapter  Three  showed  that 
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the  electron  and  hole  traps  due  to  oxygen  and  lithium  vacancies  produce  luminescence. 
Vapor  phase  equilibration  is  a  process  that  reduces  the  concentration  of  defects  by 
heating  a  crystal  in  a  material  that  will  vaporize  and  cause  ions  to  enter  the  crystal  where 
they  restore  the  perfect  crystallinity.  Undoped  LTB  crystals  could  be  submerged  in  a 
lithium-containing  powder  in  a  ceramic  boat.  The  crystal  and  powder  could  be  heated 
using  a  furnace  set  to  a  temperature  that  is  below  LTB’s  melting  point  but  high  enough 
that  lithium  ions  from  the  powder  may  fill  the  lithium  vacancies  in  the  lattice.  If  the 
lithium  vacancies  are  filled,  the  oxygen  vacancies  will  be  eliminated  at  the  same  time  so 
as  to  maintain  charge  neutrality.  This  has  been  demonstrated  in  both  lithium  niobate  and 
lithium  tantalite  crystals.  After  vapor  phase  equilibration,  LTB  could  be  exposed  to  x 
rays  at  77  K.  The  EPR  and  thermoluminescent  experiment  outlined  in  Chapter  Three 
could  then  be  used  to  determine  if  any  lithium  and  oxygen  vacancies  remain.  If  the 
vacancies  are  completely  removed  and  the  crystal  fails  to  luminesce,  the  crystal  could 
then  be  placed  in  a  thermal  neutron  environment.  If  the  lithium  and  boron  nuclei  in  the 
lattice  capture  neutrons,  defects  would  be  produced  in  the  crystal.  Post  neutron- 
irradiation  measurements  of  EPR  and  thermoluminescence  could  be  performed.  Signals 
produced  in  these  measurements  could  be  assigned  to  neutron  capture  thus  proving  that 
LTB  could  be  used  an  a  neutron  detecting  material.  These  newly  detected  signals  may 
not  match  previously  identified  EPR  spectra. 
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Appendix  A 


Experimental  EPR  and  ENDOR  Data 


Table  A.  1  Chapter  III:  LTB  trapped  hole  center  data. 


Magnetic  Field 
Direction 

EPR  Resonance 
Magnetic  Field 

ENDOR  Transition 
Frequency  (MHz) 

[001] 

3367.60 

12.356 

[001] 

3367.60 

12.460 

[001] 

3367.60 

12.601 

[101] 

3376.64 

9.070 

[101] 

3376.64 

9.516 

[101] 

3376.64 

9.980 

[101] 

3376.64 

13.170 

[101] 

3376.64 

13.460 

[101] 

3376.64 

13.760 

[101] 

3376.64 

11.600 

[101] 

3376.64 

11.650 

[101] 

3376.64 

11.700 

[100] 

3379.18 

11.226 

[100] 

3379.18 

11.351 

[100] 

3379.18 

11.445 

[100] 

3379.18 

12.459 

[100] 

3379.18 

12.662 

[100] 

3379.18 

12.858 

[110] 

3361.10 

10.720 

[110] 

3361.10 

10.920 

[110] 

3361.10 

11.170 

[110] 

3386.74 

12.849 

[110] 

3386.74 

12.973 

[110] 

3386.74 

13.081 
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Table  A. 2  Chapter  IV:  LTB:Ag  unperturbed  trapped  hole  center  data. 


Magnetic  Field 
Direction  /  Angle  (°) 

EPR  Resonance 
Magnetic  Field 

Data  for  a  to  a  plane 

[001] 

3087.55 

[001] 

3108.098 

[001] 

3136.258 

[001] 

3158.562 

15 

3059.343 

15 

3079.281 

15 

3100.412 

15 

3137.968 

15 

3162.058 

15 

3184.938 

25 

3048.214 

25 

3099.262 

25 

3141.506 

25 

3161.461 

25 

3175.549 

25 

3197.849 

35 

3044.109 

35 

3061.719 

35 

3081.046 

35 

3161.461 

35 

3180.849 

35 

3201.985 

45 

3048.193 

45 

3068.729 

45 

3176.141 

45 

3197.27 
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Table  A. 2  (cont’d)  Chapter  IV:LTB:Ag  unperturbed  trapped  hole  center  data. 


Magnetic  Field  EPR  Resonance 

Direction  /  Angle  (°)  Magnetic  Field 

Data  for  c  to  a  plane 


Magnetic  Field  EPR  Resonance 

Direction  /  Angle  (°)  Magnetic  Field 

Data  for  c  to  a  plane  (cont'd) 


[001] 

2977.315 

60 

2889.682 

[001] 

3001.566 

60 

2912.370 

60 

3058.719 

10 

2919.952 

60 

3082.192 

10 

2945.757 

60 

3137.751 

10 

2969.622 

60 

3160.452 

10 

2995.052 

60 

3269.977 

10 

3019.676 

60 

3285.624 

10 

3046.700 

10 

3070.161 

70 

2941.180 

70 

2963.867 

20 

2872.812 

70 

3087.515 

20 

2897.551 

70 

3111.756 

20 

2969.104 

70 

3149.313 

20 

2993.309 

70 

3171.223 

20 

3021.815 

70 

3230.675 

20 

3044.942 

70 

3247.896 

20 

3125.644 

20 

3148.230 

80 

3014.719 

80 

3035.839 

30 

2841.827 

80 

3116.812 

30 

2867.058 

80 

3140.197 

30 

2980.335 

80 

3147.324 

30 

3002.628 

80 

3170.805 

30 

3053.110 

80 

3162.000 

30 

3077.184 

80 

3179.009 

30 

3199.845 

30 

3219.778 

[001] 

3087.550 

[001] 

3108.098 

45 

2842.557 

[001] 

3136.258 

45 

2866.801 

[001] 

3158.562 

45 

3012.331 

45 

3035.806 

45 

3099.145 

45 

3122.617 

45 

3267.361 

45 

3285.328 

T 


Table  A. 3  Chapter  IV:  LTB:Ag  perturbed  trapped  hole  center. 


Magnetic  Field 
Direction  /  Angle  (°) 

EPR  Resonance 
Magnetic  Field 

Data  for  a  to  a  plane 

[001] 

3042.965 

[001] 

3058.819 

[001] 

3294.133 

[001] 

3304.680 

15 

3042.308 

15 

3261.255 

15 

3272.969 

15 

3288.831 

15 

3299.390 

25 

3110.997 

25 

3224.281 

25 

3237.190 

25 

3266.529 

25 

3277.682 

35 

3091.055 

35 

3105.149 

35 

3136.245 

35 

3150.346 

35 

3193.166 

35 

3231.323 

35 

3243.043 

45  [110] 

3132.114 

45  [110] 

3146.786 

45  [110] 

3185.534 

45  [110] 

3197.270 
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Table  A. 3  (cont’d)  Chapter  IV:  LTB:Ag  perturbed  trapped  hole  center. 


Magnetic  Field 
Direction  /  Angle  (°) 


EPR  Resonance 
Magnetic  Field 


Data  for  c  to  a  plane 


Magnetic  Field  EPR  Resonance 

Direction  /  Angle  (°)  Magnetic  Field 

Data  for  c  to  a  plane 


[001] 

2931.543 

50 

3089.195 

[001] 

2951.301 

50 

3105.625 

50 

3159.608 

10 

2889.426 

50 

3172.122 

10 

2911.328 

50 

3186.214 

10 

2931.281 

10 

2981.737 

60 

3172.171 

60 

3231.667 

20 

2858.818 

60 

3244.939 

20 

2945.435 

20 

3008.362 

70 

2932.574 

20 

3057.329 

70 

2952.130 

70 

3139.145 

30 

3094.174 

70 

3215.042 

30 

3110.016 

70 

3276.065 

70 

3287.795 

40 

3047.808 

40 

3120.576 

80 

2988.303 

40 

3134.660 

80 

3005.893 

40 

3149.529 

80 

3093.348 

80 

3108.603 

[101] 

3056.118 

80 

3267.043 

[101] 

3073.326 

80 

3279.369 

[101] 

3135.903 

80 

3298.719 

[101] 

3150.763 

80 

3309.865 

[101] 

3163.281 

[100] 

3042.965 

[100] 

3058.819 

[100] 

3294.133 

[100] 

3304.68 
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Sample  EPR  and  ENDOR  Computer  Routines 


%  EnergyLevels_g_and_A 

% 

%  This  subroutine  is  used  in  conjunction  with  Cntr_A_EPR_linepositions  to 
%  calculate  the  line  positions  as  a  function  of  angle  for  Center  A  in  LiB305. 

% 

%  It  calculates  the  eigenvalues  and  returns  them  to  the  main  program, 
function  EE  =  EnergyLevels_g_and_A(TG,TH,TQ,HH,P,B,gbnj 


W1  =B*HH*(P(1)*TG(1,1)*TG(1,3)+P(2)*TG(2,1)*TG(2,3)+P(3)*TG(3,1)*TG(3,3)); 

W2  =  B*HH*(P(  1)*TG(  1 ,2)*TG(  1 ,3)+P(2)*TG(2,2)*TG(2,3)+P(3)*TG(3.2)*TG(3,3)); 

W3  =  B*HH*(P(lj*TG(l,3)*TG(l,3)+P(2)*TG(2,3)*TG(2,3)+P(3)*TG(3,3)*TG(3,3)); 

W4  =  P(7)*TH(  1 , 1  )*TH(  1 . 1  )+P(8)*TH(2, 1  )*TH(2,1)+P(9)*TH(3, 1  )*TH(3, 1); 

W5  =  P(7)*TH(  1 ,2)*TH(  1 . 1  )+P(8)*TH(2,2)*TH(2, 1  )+P(9)*TH(3,2)*TH(3, 1 ); 

W6  =  P(7)*TH(  1 ,3)*TH(  1 , 1  )+P(8)*TH(2,3)*TH(2, 1  )+P(9)*TH(3,3)*TH(3, 1 ); 

W7  =  P(7)*TH(  1 ,2)*TH(  1 ,2)+P(8)*TH(2,2)*TH(2,2)+P(9)*TH(3.2)*TH(3.2); 

W8  =  P(7)*TH(  1 ,3)*TH(  1 ,2)+P(8)*TH(2,3)*TH(2,2)+P(9)*TH(3,3)*TH(3,2); 

W9  =  P(7)*TH(1,3)*TH(1,3)+P(8)*TH(2,3}*TH(2,3)+P(9)*TH(3,3)*TH(3,3); 

QZZ  =  -(P(13)+P(14j); 

W10  =  P(T3)*TQ(l,l)*TQ(l,l)+P(14)*TQ(2,lj*TQ(2,l)+QZZ*TQ(3,l)*TQ(3,l); 

W1 1  =  P(  13)*TQ(  1 ,2)*TQ(  1 , 1  )+P(  14)*TQ(2,2)*TQ(2,1  )+QZZ*TQ(3,2)*TQ(3, 1 ); 

W12  =  P(13)*TQ(1,3)*TQ(1,1)+P(14)*TQ(2,3)*TQ(2,1)+QZZ*TQ(3,3)*TQ(3,1); 

W13  =  P(13)*TQ(1.2)*TQ(1,2)+P(14)*TQ(2,2)*TQ(2,2)+QZZ*TQ(3,2)*TQ(3,2); 

W14  =  P(13)*TQ(  1 ,3)*TQ(  1 ,2)+P(  14)*TQ(2,3)*TQ(2,2)+QZZ*TQ(3,3)*TQ(3,2); 

W15  =  P(13)*TQ(1,3)*TQ(1,3)+P(14)*TQ(2,3)*TQ(2,3)+QZZ*TQ(3,3)*TQ(3,3); 


Q1  =0.5*(Wl+i*W2); 

Q2  =  0.25*(W4-W7)+0.5*i*W5; 

Q3  =  0.25*(W4+W7); 

Q4  =  0.5*(W6+i*W8); 

Q5  =  0.25*(W10-W13)+0.5*i*Wl  1; 
Q6  =  0.25*(W10+W13); 

Q7  =  0.5*(W12+i*W14); 


%  HAM  is  the  matrix  representing  the  spin-Hamiltonian 


Ham  =  zeros(8); 

Ham(l,l)  =  0.5*W3  +  0.75*W9  +  2.25*W15  -  1.5*gbn*HH  +  3*Q6; 
Ham(2,2)  =  0.5*W3  +  0.25*W9  +  0.25*W15  -  0.5*gbn*HH  +  7*Q6; 
Ham(3.3)  =  0.5*W3  -  0.25*W9  +  0.25*W15  +  0.5*gbn*HH  +  7*Q6; 
Ham(4,4)  =  0.5*W3  -  0.75*W9  +  2.25*W15  +  1.5*gbn*HH  +  3*Q6; 
Ham(5.5)  =  -0.5*W3  -  0.75*W9  +  2.25*W15  -  1.5*gbn*HH  +  3*Q6; 
Ham(6.6)  =  -0.5*W3  -  0.25*W9  +  0.25*W15  -  0.5*gbn*HH  +  7*Q6; 
Ham(7,7)  =  -0.5*W3  +  0.25*W9  +  0.25*W15  +  0.5*gbn*HH  +  7*Q6; 
Ham(8,8)  =  -0.5*W3  +  0.75*W9  +  2.25*W15  +  1.5*gbn*HH  +  3*Q6; 
Ham(2,l)  =  0.5*sqrt(3)*Q4  +  2*sqrt(3)*Q7; 

Ham(3.1)  =  2*sqrt(3)*Q5; 

Ham(3,2)  =  Q4; 

Ham(4,2)  =  2*sqrt(3)*Q5; 

Ham(4,3)  =  0.5*sqrt(3)*Q4  -  2*sqrt(3)*Q7; 

Ham(5.1)  =  Ql  +  1.5*Q4; 

Ham(5.2)  =  sqrt(3)*Q3; 

Ham(6,l)  =  sqrt(3)*Q2; 

Ham(6,2)  =  Ql  +0.5*Q4; 

Ham(6.3)  =  2*Q3; 

Ham(6,5)  =  -0.5*sqrt(3)*Q4  +  2*sqrt(3)*Q7; 

Ham(7,2)  =  2*Q2; 

Ham(7.31  =  Ql  -  0.5*Q4; 

Ham(7,4)  =  sqrt(3)*Q3; 

Ham(7,5)  =  2*sqrt(3)*Q5; 

Ham(7.6)  =  -Q4; 

Ham(8,3)  =  sqrt(3)*Q2; 

Ham(8,4)  =  Ql  -  1.5*Q4; 
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Ham(8,6) 
Ham(8,7) 
Ham(  1 ,2) 
Ham(l,3) 
Ham(2,3) 
Ham(2,4) 
Ham(3,4) 
Ham(l,5) 
Ham(2,5) 
Ham(  1 .6) 
I  lam(2.6) 
I  lam(3.6) 
Ham(5.6) 
Ham(2,7) 
Ham(3,7) 
Ham(4,7) 
Ham(5.7) 
Ham(6,7) 
Ham(3,8) 
Ham(4,8) 
Ham(6,8) 
Ham(7,8) 


2*sqrt(3)*Q5; 

-0.5*sqrt(3)*Q4  -2*sqrt(3)*Q7; 

conj(Ham(2,l)); 

conj(Ham(3,l)); 

conj(Ham(3,2)); 

conj(Ham(4,2)); 

conj(Ham(4,3)); 

conj(Ham(5,l)); 

conj(Ham(5,2)); 

conj(Ham(6,l)); 

conj(Ham(6,2)); 

conj(Ham(6,3)); 

conj(Ham(6,5)); 

conj(Ham(7,2)); 

conj(Ham(7,3)); 

conj(Ham(7,4)); 

conj(Ham(7,5)); 

conj(Ham(7,6)); 

conj(Ham(8,3)); 

conj(Ham(8,4)); 

conj(Ham(8,6)); 

conj(Ham(8,7)); 


EE  =  sort(real(eig(Ham))); 


%  Lithium  Tetraborate:  Subroutine  for  hyperfine  and 

%  quadrupole  parameters  fitting. 

% 

%  This  subroutine  is  used  in  conjunction  with  Li_vacancy_ENDOR_fitting 
%  to  determine  the  best  set  of  1  IB  hyperfine  parameters  for  the 
%  Li+  vacancy  hole  center  in  Li2B407. 

%  It  calculates  a  sum  of  the  frequency  differences  squared  and 
%  returns  the  value  to  the  main  program.  The  input  data  are  the 
%  measured  magnetic  fields  and  ENDOR  frequencies. 

function  summ  =  Li_vacancy_ENDOR_fitting_sub(PP,P,B,gbn) 

CTR  =  pi/180; 

%  G  is  the  3x3  rotation  matrix  which  takes  the  principal 
%  axes  of  the  g  matrix  into  the  crystal  coordinate  system. 

%  H  is  the  3x3  rotation  matrix  which  takes  the  principal 
%  axes  of  the  A  matrix  into  the  crystal  coordinate  system. 

%  Q  is  the  3x3  rotation  matrix  which  takes  the  principal 
%  axes  of  the  Q  matrix  into  the  crystal  coordinate  system. 

%  R  is  the  3x3  rotation  matrix  which  takes  the  crystal  coordinate 
%  system  into  the  magnetic  field  coordinate  system. 

G(l,l)  =  cos(PP(6))*cos(PP(5))  -  cos(PP(4))*sin(PP(5))*sin(PP(6)); 

G(l,2)  =  cos(PP(6))*sin(PP(5))  +  cos(PP(4))*cos(PP(5))*sin(PP(6)); 

G(  1,3)  =  sin(PP(6))*sin(PP(4)); 

G(2,l)  =  -sin(PP(6))*cos(PP(5))  -  cos(PP(4))*sin(PP(5))*cos(PP(6)); 
G(2,2)  =  -sin(  PP(6))*sin(PP(5 ))  +  cos(PP(4))*cos(PP(5))*cos(PP(6)); 
G(2,3)  =  cos(PP(6))*sin(PP(4)); 

G(3,l)  =  sin(PP(4))*sin(PP(5)); 

G(3,2)  =  -sin(PP(4))*cos(PP(5)); 

G(3,3)  =  cos(PP(4)); 

Hf  1,1)  =  cos(P(6) )*cos(P(5))  -  cos(P(4))*sin(P(5))*sin(P(6)); 

Hf  1,2)  =  cos(P(6))*sin(P(5))  +  cos(P(4))*cos(P(5))*sin(P(6)); 

H(l,3)  =  sin(P(6))*sin(P(4)); 

H(2,l)  =  -sin(P(6 ))*cos(P(5))  -  cos(P(4))*sin(P(5))*cos(P(6)); 

H(2,2)  =  -sin(P(6))*sin(P(5))  +  cos(P(4))*cos(P(5))*cos(P(6)); 

H(2,3)  =  cos(P(6))*sin(P(4)); 
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H(3,l)  =  sin(P(4))*sin(P(5 )); 

H(3,2)  =  -sin(P(4))*cos(P(5)); 

H(3,3)  =  cos(P(4)); 

Q(l,l)  =  cos(P(  1 1  ))*cos(P(  1 0))  -  cos(P(9))*sin(P(10)]*sin(P(ll)); 

Q(  1 ,2)  =  cos(P(  1 1  ))*sin(P(  1 0))  +  cos(P(9))*cos(P(10))*sin(P(ll)); 

Q(1 ,3)  =  sin(P(  1 1  ))*sin(P(9)); 

Q(2,l)  =  -sin(P(  1 1  ))*cos(P(  10))  -  cos(P(9))*sin(P(10))*cos(P(ll)); 

Q(2,2)  =  -sin(P(l  l))*sin(P(10))  +  cos(P(9))*cos(P(10))*cos(P(l  1)); 

Q(2,3)  =  cos(P(l  1  ))*sin(P(9)); 

Q(3,l)  =  sin(P(9))*sin(P(10)); 

Q(3,2)  =  -sin(P(9 )  )*cos(P(  10)); 

Q(3,3)  =  cos(P(9)); 

%  Data  for  [001] 

h(l)=3367.6;FRQ(l)=12.356;Line(l)=3;Alpha(l)=0*CTR;Beta(l)=0*CTR; 

h(2)=3367.6;FRQ(2)=12.46;Line(2)=2;Alpha(2)=0*CTR;Beta(2)=0*CTR; 

h(3)=3367.6;FRQ(3)=12.601;Line(3)=l;Alpha(3)=0*CTR;Beta(3)=0*CTR; 

%  Data  for  [101] 

h(4)=3376.64;FRQ(4)=9.07;Line(4)=l;Alpha(4)=45*CTR;Beta(4)=0*CTR; 

h(5)=3376.64;FRQ(5)=9.516:Line(5)=2;Alpha(5)=45*CTR;Beta(5)=0*CTR; 

h(6)=3376.64;FRQ(6)=9.98;Line(6)=3;Alpha(6)=45*CTR;Beta(6)=0*CTR; 

h(7)=3376.64;FRQ(7)=13.17;Line(7)=3;Alpha(7)=45*CTR;Beta(7)=90*CTR; 

h(8)=3376.64;FRQ(8)=13.46;Line(8)=2;Alpha(8)=45*CTR;Beta(8)=90*CTR; 

h(91=3376.64;FRQ(9)=13.76;Line(9)=l;Alpha(9)=45*CTR;Beta(9)=90*CTR; 

h(10)=3376.64;FRQ(10)=11.6;Line(10)=3;Alpha(10)=45*CTR;Beta(10)=270*CTR; 

h(  1 1  )=3376.64;FRQ(  1 1  )=1 1 .65  ;Line(  1 1  )=2 ;  Alpha(  1 1  )=45  *CTR;Beta(  1 1  )=270*CTR; 

h(12)=3376.64;FRQ(12)=11.7;Line(12)=l;Alpha(12)=45*CTR;Beta(12)=270*CTR; 

%  Data  for  [100] 

h(13)=3379.18;FRQ(13)=11.226;Line(13)=3;Alpha(13)=90*CTR;Beta(13]=0*CTR; 

h(14)=3379.18;FRQ(14)=11.351;Line(14]=2;Alpha(14)=90*CTR;Beta(14)=0*CTR; 

h(15)=3379.18;FRQ(15)=11.445;Line(15)=l;Alpha(15)=90*CTR;Beta(15)=0*CTR; 

h(16)=3379.18;FRQ(16)=12.459;Line(16]=3;Alpha(16)=90*CTR;Beta(16)=90*CTR; 

h(17)=3379.18;FRQ(17)=12.662;Line(17)=2;Alpha(17)=90*CTR;Betaa7)=90*CTR; 

h(18)=3379. 1 8;FRQ(  1 8)=12.858;Line(  1 8)=1  ;Alpha(  1 8)=90*CTR;Beta(  1 8)=90*CTR; 

%  Data  for  [110] 

h(19)=3361.10;FRQ(19]=10.72;Line(19)=3;Alpha(19)=90*CTR;Beta(19)=45*CTR; 
h(20)=3361.10;FRQ(20)=10.92;Line(20)=2;Alpha(20)=90*CTR;Beta(20)=45*CTR; 
h(2 1  )=336 1 . 1 0;FRQ(2 1  )=1 1 . 1 7;Line(2 1  )=  1 ;  Alpha(2 1  )=90*CTR;Beta(2 1  )=45  *CTR; 
h(22)=3386.74;FRQ(22)=12.849;Line(22)=3;Alpha(22)=90*CTR;Beta(22)=135*CTR; 
h(23)=3386.74;FRQ(23)=12.973;Line(23)=2;Alpha(23)=90*CTR;Beta(231=135*CTR; 
h(24)=3386.74;FRQ(24 )=13.08  l;Line(24]=l;Alpha(24)=90*CTR;Beta(24)=135*CTR; 

%Data  for  [100+10] 

%h(25)=3373.0;FRQ(25)=10.88;Line(25)=3;Alpha(25)=90*CTR;Beta(25)=10*CTR; 
%h(26)=3373.0;FRQ(26)=11.02;Line(26)=2;Alpha(26)=90*CTR;Beta(26)=10*CTR; 
%h(27)=3373.0;FRQ(27)=  11.17  ;Line(27)= 1  ;Alpha(27)=90*CTR;Beta(27)=  1 0*CTR; 
%h(28)=3373.0;FRQ(28)=12.79;line(28)=3;Alpha(28)=90*CTR;Beta(28)=100*CTR; 
%h(29)=3373.0;FRQ(29)=12.99;Line(29)=2;Alpha(29)=90*CTR;Beta(29)=100*CTR; 
%h(30)=3373.0;FRQ(30)=  13.15  ;Une(30)=  1  ;Alpha(30)=90*CTR;Beta(30)=  1 00*CTR; 

datapoints  =  length(h); 

for  nn=l:  datapoints 
HH  =  h(nn); 
line  =  Line(nn); 
alpha  =  Alpha(nn); 
beta  =  Beta(nn); 

Rf  1,1)  =  cos(alpha)*cos(beta); 

R(l,2)  =  -sin(beta); 

R(l,3)  =  sin(alpha)*cos(beta); 

R(2,l)  =  cos(alpha)*sin(beta); 

R(2,2)  =  cos(beta); 

R(2,3)  =  sin(alpha)*sin(beta); 

R(3,l)  =  -sin(alpha]; 
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R(3,2)  =  0; 

R(3,3)  =  cos(alpha); 

TG  =  G  *  R; 

TH  =  H  *  R; 

TQ  =  Q  *  R; 


W1  =  B*HH*(PP(1  )*TG(  1 , 1  )*TG(  1 ,3)+PP(2)*TG(2, 1  )*TG(2,3  )+PP(3)*TG(3, 1)*TG(3,3)); 
W2  =  B*HH*(PP(1)*TG(  1 ,2)*TG(  1 ,3  )+PP(2)*TG(2,2)*TG(2,3  )+PP(3)*TG(3,2)*TG(3,3)); 
W3  =  B*HH*(PP(1)*TG(1,3)*TG(1,3)+PP(2)*TG(2,3)*TG(2,3)+PP(3)*TG(3,3)*TG(3.3)); 
W4  =  P(  1)*TH(  1 , 1)*TH(1 , 1  )+P(2)*TH(2, 1)*TH(2, 1  )+P(3)*TH(3. 1  )*TH(3. 1 ); 

W5  =  P(  1  )*TH(  1 , 1)*TH(1 ,2)+P(2)*TH(2, 1)*TH(2,2)+P(3)*TH(3. 1  )*TH(3.2); 

W6  =  P(  1  )*TH(  1 , 1  )*TH(  1 ,3)+P(2)*TH(2, 1)*TH(2,3)+P(3)*TH(3, 1  )*TH(3,3); 

W7  =  P(  1  )*TH(  1 ,2)*TH(1,2)+P(2)*TH(2,2)*TH(2,2)+P(3)*TH(3,2)*TH(3,2); 

W8  =  P(  1  )*TH(  1 ,2)*TH(  1 ,3)+P(2)*TH(2,2)*TH(2,3)+P(3)*TH(3,2)*TH(3,3); 

W9  =  P(  1  )*TH(  1 ,3)*TH(1 ,3)+P(2)*TH(2,3)*TH(2,3)+P(3)*TH(3,3)*TH(3.3); 

QZZ  =  -(P(7)+P(8)); 

W10  =  P(7)*TQ(1,1)*TQ(1,1)+P(8)*TQ(2,1)*TQ(2,1)+QZZ*TQ(3,1)*TQ(3,1); 

W1 1  =  P(7)*TQ(1 ,2)*TQ(  1 , 1  )+P(8)*TQ(2,2)*TQ(2, 1  )+QZZ*TQ(3,2)*TQ(3, 1 ); 

W12  =  P(7)*TQ(1 ,3)*TQ(  1 , 1  )+P(8)*TQ(2,3)*TQ(2, 1  )+QZZ*TQ(3,3)*TQ(3, 1 ); 

W13  =  P(7)*TQ(1,2)*TQ(1,2)+P(8)*TQ(2,2)*TQ(2,2)+QZZ*TQ(3,2)*TQ(3,2); 

W14  =  P(7)*TQ(1,3)*TQ(1,2)+P(8)*TQ(2,3)*TQ(2,2)+QZZ*TQ(3,3)*TQ(3,2); 

W15  =  P(7)*TQ(1,3)*TQ(1,3)+P(B)*TQ(2,3)*TQ(2,3)+QZZ*TQ(3,3)*TQ(3,3); 

Q1  =  0.5*(Wl+i*W2); 

Q2  =  0.25*(W4-W7)+0.5*i*W5; 

Q3  =  0.25*(W4+W7 ); 

Q4  =  0.5*(W6+i*W8); 

Q5  =  0.25*(W10-W13)+0.5*i*Wl  1; 

Q6  =  0.25*(W10+W13); 

Q7  =  0.5*(W12+i*W14); 

%  Ham  is  the  matrix  representing  the  spin-Hamiltonian: 


Ham  =  xeros(8); 

Ham(l.l)  =  0.5*W3  +  0.75*W9  +  2.25*W15  -  1.5*gbn*HH  +  3*Q6; 
Ham(2,2)  =  0.5*W3  +  0.25*W9  +  0.25*W15  -  0.5*gbn*HH  +  7*Q6; 
Ham(3,3)  =  0.5*W3  -  0.25*W9  +  0.25*W15  +  0.5*gbn*HH  +  7*Q6; 
Ham(4,4)  =  0.5*W3  -  0.75*W9  +  2.25*W15  +  1.5*gbn*HH  +  3*Q6; 
Ham(5,5)  =  -0.5*W3  -  0.75*W9  +  2.25*W15  -  1.5*gbn*HH  +  3*Q6; 
Ham(6,6)  =  -0.5*W3  -  0.25*W9  +  0.25*W15  -  0.5*gbn*HH  +  7*Q6; 
Ham(7,7)  =  -0.5*W3  +  0.25*W9  +  0.25*W15  +  0.5*gbn*HH  +  7*Q6; 
Ham(8,8)  =  -0.5*W3  +  0.75*W9  +  2.25*W15  +  1.5*gbn*HH  +  3*Q6; 
Ham(2,l)  =  0.5*sqrt(3)*Q4  +  2*sqrt(3)*Q7; 

Ham(3,l)  =  2*sqrt(3)*Q5; 

Ham(3.2)  =  Q4; 

Ham(4,2)  =  2*sqrt(3)*Q5; 

Ham(4,3)  =  0.5*sqrt(3)*Q4  -  2*sqrt(3)*Q7; 

Ham(5,l)  =  Ql  +  1.5*Q4; 

Ham(5,2)  =  sqrt(3)*Q3; 

Ham(6,l)  =  sqrt(3)*Q2; 

Ham(6,2)  =  Ql  +0.5*Q4; 

Ham(6,3)  =  2*Q3; 

Ham(6,5)  =  -0.5*sqrt(3)*Q4  +  2*sqrt(3)*Q7; 

Ham(7,2)  =  2*Q2; 

Ham(7,3)  =  Ql  -  0.5*Q4; 

Ham(7,4)  =  sqrt(3)*Q3; 

Ham(7,5)  =  2*sqrt(3)*Q5; 

Ham(7,6)  =  -Q4; 

Ham(8,3)  =  sqrt(3)*Q2; 

Ham(8,4)  =  Ql  -  1.5*Q4; 

Ham(8,6)  =  2*sqrt(3)*Q5; 

Ham(8,7)  =  -0.5*sqrt(3)*Q4  -2*sqrt(3)*Q7; 

Ham(l,2)  =  conj(Ham(2,l)); 

Ham(l,3)  =  conj(Ham(3,l)); 

Ham(2,3)  =  conj(Ham(3,2)); 

Ham(2,4)  =  conj(Ham(4,2)); 

Ham(3,4)  =  conj(Ham(4,3)); 

Ham(l,5)  =  conj(Ham(5,l)); 

Ham(2,5)  =  conj(Ham(5,2)); 
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Ham(l,6)  =  conj(Ham(6,l)); 

Ham(2,6)  =  conj(Ham(6,2)); 

Ham(3,6)  =  conj(Ham(6,3)); 

Ham(5,6)  =  conj(Ham(6,5)); 

Ham(2,7)  =  conj(Ham(7,2)); 

Ham(3,7)  =  conj(Ham(7,3)); 

Ham(4,7)  =  conj(Ham(7,4)); 

Ham(5,7)  =  conj(Ham(7,5)); 

Ham(6,7)  =  conj(Ham(7,6)); 

Ham(3,8)  =  conj(Ham(8,3)); 

Ham(4,8)  =  conj(Ham(8,4)); 

Ham(6,8)  =  conj(Ham(8,6)); 

Ham(7,8)  =  conj(Ham(8,7)); 

EE  =  sort(real(eig(Ham))); 
if  line  ==  1 

Freq(nn)  =  abs(EE(2)  -  EE(1)); 
elseif  line  ==  2 

Freq(nn)  =  abs(EE(3)  -  EE(2)); 
elseif  line  ==  3 

Freq(nn)  =  abs(EE(4)  -  EE(3)); 
elseif  line  ==  4 

Freq(nn)  =  abs(EE(6)  -  EE(5)); 
elseif  line  ==  5 

Freq(nn)  =  abs(EE(7)  -  EE(6)); 
elseif  line  ==  6 

Freq(nn)  =  abs(EE(8)  -  EE(7)); 
end 
end 

summ=0; 

for  ii=l  :datapoints 

summ  =  summ  +  (FRQ(ii)-Freq(ii))A2; 
end 
end 


%  EPR_fitting_gandA_P 

%  This  program  determines  the  "best"  set  of  g  and  A  parameters  for  the 
%  perturbed  Ag2+  hole  center  in  Li2B407.  Because  lines  from  the  two 
%  isotopes  of  Ag  (107  and  109)  are  not  resolved  in  the  EPR  spectra,  an 
%  "averaged"  A  matrix  representing  the  two  isotopes  together  is  used  in 
%  this  fitting  program. 

%  Input  data  consists  of  74  EPR  magnetic  field  values  and  their  corresponding 
%  microwave  frequencies.  The  output  is  12  parameters  (three  principal  values 
%  and  three  Euler  angles  for  both  g  and  A). 

%  This  program  uses  a  subroutine  named  SUMM_EPR_fitting_gandA_P. 
clear  all 
%  Constants: 

h  =  6.626069;  %  Planck's  constant 

B  =  9.274009/h;  %  Bohr  magneton  divided  by  Planck's  constant 

gbn  =  -0.00018637;  %  averaged  value  of  gn*bn  for  107 Ag  and  109 Ag 

CTR  =  pi/180;  %  Conversion  constant,  degrees  to  radians 

%  Initial  values  of  the  spin-Hamiltonian  parameters: 

%  Six  for  the  g  matrix  (three  principal  values  and  three  angles). 

%  Six  for  the  A  matrix  (three  principal  values  and  three  angles). 

P(l)  =  2.0329; 

P(2)  =  2.1215; 

P(3)  =  2.3607; 

P(4)  =  35.7*CTR; 

P(5)  =  9.7*CTR; 

P(6)  =  -9.9*CTR; 
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P(7)  =  27.31; 

P(8)  =  34.13; 

P(9)  =  66.04; 

P(  10)  =  33.4*CTR; 

P(1 1)  =  26.0*CTR; 

P(12)  =  -40.6*CTR; 

%  Step  sizes  for  the  parameters: 

gg  =  0.0001 ;  %  step  size  for  the  g  principal  values 

aa  =  0.01 ;  %  step  size  for  the  A  principal  values 

delta  =  0.1  *CTR;  %  step  size  for  the  angles 

step(l)  =  gg; 
step(2)  =  gg; 
step(3)  =  gg; 
step(4)  =  delta; 
step(5)  =  delta; 
step(6)  =  delta; 
step(7)  =  aa; 
step(8)  =  aa; 
step(9)  =  aa; 
step(10)  =  delta; 
step(ll)  =  delta; 
step(12)  =  delta; 

sum2  =  0; 

suml  =  SUMM_EPR_fitting_gandA_P(P,B,gbn); 

while  sum2<suml 
for  n  =  1:12 

summ  =  SUMM_EPR_fitting_gandA_P(P,B,gbn); 
sum2  =  summ; 
if  n=l; 

suml  =  summ; 
end 

P(n)  =  P(n)  +  step(n); 

summ  =  SUMM_EPR_fitting_gandA_P(P,B,gbn); 
if  summ  >=  sum2; 

P(n)  =  P(n)  -  2*step(n); 

summ  =  SUMM_EPR_fitting_gandA_P(P,B,gbn); 
if  summ  >=  sum2; 

P(n)  =  P(n)  +  step(n); 
end 
end 
end 

if  summ<sum2; 

sum2  =  summ; 
end 
sum2 
end 

P(4)  =  P(4)/CTR; 

P(5)  =  P(5)/CTR; 

P(6)  =  P(6)/CTR; 

P(10)  =  P(10)/CTR; 

P(1 1)  =  P(1 1)/CTR; 

P(12)  =  P(12)/CTR; 

P  %  Display  final  set  of  parameters. 

sum2  %  Display  final  value  of  sum2. 

%  End  of  program. 
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SUMM_EPR_fitting_gandA_P 


%  This  subroutine  is  used  in  conjunction  with  EPR_fitting_gandA  to 
%  determine  the  best  set  of  g  and  "averaged"  A  hyperfine  parameters 
%  for  the  perturbed  Ag2+  trapped  hole  center  in  Li2B407. 

%  It  calculates  a  sum  of  the  frequency  differences  squared  and 
%  returns  the  value  to  the  main  program.  The  input  data  are  the 
%  measured  magnetic  fields  and  microwave  frequencies. 

function  summ  =  SUMM_EPR_fitting_gandA_P(P,B,gbn) 

CTR  =  pi/180; 

%  G  is  the  3x3  rotation  matrix  which  takes  the  principal 
%  axes  of  the  g  matrix  into  the  crystal  coordinate  system. 

%  Taken  from  "Classical  Mechanics,  2nd  ed."  by  Goldstein,  pp.  146-147. 

%  H  is  the  3x3  rotation  matrix  which  takes  the  principal 
%  axes  of  the  A  matrix  into  the  crystal  coordinate  system. 

%  R1  is  the  3x3  rotation  matrix  which  takes  sites  1  through  8 
%  into  the  crystal  coordinate  system. 

%  R  is  the  3x3  rotation  matrix  which  takes  the  crystal  coordinate 
%  system  into  the  magnetic  field  coordinate  system. 

G(  1,1)  =  cos(P(6))*cos(P(5))  -  cos(P(41)*sin(P(5))*sin(P(6)); 

G(l,2)  =  cos(P(6)  l*sin(P(5))  +  cos(P(4))*cos(P(5))*sin(P(6)); 

G(l,3)  =  sin(P(6))*sin(P(4)); 

G(2,l)  =  -sin(P(6 l)*cos(P(5)l  -  cos(P(4))*sin(P(5))*cos(P(6)); 

G(2,2)  =  -sin(P(6))*sin(P(5))  +  cos(P(4))*cos(P(5))*cos(P(6)); 

G(2,3)  =  cos(P(6))*sin(P(4)); 

G(3,l)  =  sin(P(4))*sin(P(5)l; 

G(3,2)  =  -sin(P(4 )  l*cos(P(51); 

G(3,3)  =  cos(P(4)); 

H(l,l)  =  cos(P(  1 2))*cos(P(  1 1 ))  -  cos(P(  1 0))*sin(P(T  1  ))*sin(P(  12)); 

H(l,2)  =  cos(P(  1 2  ))*sin(P(  1 1 ))  +  cos(P(10))*cos(P(l  l))*sin(P(12)); 

Hf  1,3)  =  sin(P(12))*sin(P(10)); 

H(2,l)  =  -sin(P(12))*cos(P(l  1))  -  cos(P(10))*sin(P(l  1  ))*cos(P(12)l; 

H(2,2)  =  -sin(P(12))*sin(P(l  1 ))  +  cos(P(10))*cos(P(l  l))*cos(P(12)); 

H(2,3)  =  cos(P(  1 2))*sin(P(  10)); 

H(3,l )  =  sin(P(10))*sin(P(l  1)); 

H(3,2)  =  -sin(P(10))*cos(P(l  1)); 

H(3,3)  =  cos(P(10j); 

%  Rotation  from  [001]  to  [100],  These  data  were  taken  on  July  24  and  25,  2012. 

h(  1  )=293 1 .543  ;FRQ(  1  )=9399.820;  Alpha(  1  )=0*CTR;Beta(  1  )=0*CTR;K(  1  )=  1  ;IC(  1  )=  1 ; 
h(2)=295 1 .301  ;FRQ(2)=9399.820;Alpha(2)=0*CTR;Beta(2  )=0*CTR;K(2)=  1  ;IC(2)=2; 

h(3)=2889.426;FRQ(3)=9399.298;Alpha(3)=10*CTR;Beta(3)=0*CTR;K(3)=2;IC(3)=l; 
h(4J=2911.328;FRQ(4)=9399.298;Alpha(4)=10*CTR;Beta(4)=0*CTR;K(4)=2;IC(4)=2; 
h(5)=293 1 .28 1  ;FRQ(5)=9399.298;Alpha(5)=10*CTR;Beta(5)=0*CTR;K(5)=l  ;IC(5)=1 ; 
h(6)=2981.737;FRQ(6)=9399.298;Alpha(6)=10*CTR;Beta(6)=0*CTR;K(6)=4;IC(6)=l; 

h(7)=2858.818;FRQ(7)=9399.706;Alpha(7)=20*CTR;Beta(7)=0*CTR;K(7)=2;IC(7)=l; 
h(8)=2945.435;FRQ(8)=9399.706;Alpha(8)=20*CTR;Beta(8)=0*CTR;K(8)=l;IC(8)=l; 
h(9)=3008.362;FRQ(9)=9399.706;Alpha(9)=20*CTR;Beta(9)=0*CTR;K(9)=3;IC(9)=l; 
h(  10)=3057.329;FRQ(  10)=9399.706;Alpha(10)=20*CTR;Beta(10)=0*CTR;K(  10)=4;IC(  1 0)=2; 

h(l  1)=3094.174;FRQ(1 1  )=9400.099 ;Alpha(l  l)=30*CTR;Beta(l  1)=0*CTR;K(1 1)=4;IC(1 1)=1; 
h(12)=3110.016;FRQ(12)=9400.099;Alpha(12]=30*CTR;Beta(12)=0*CTR;K(12)=4;IC(12)=2; 

h(13)=3047.808;FRQ(13)=9400.255;Alpha(13)=40*CTR;Beta(13)=0*CTR;K(13)=l;IC(  131=2; 
h(14)=3120.576;FRQ(14)=9400.255;Alpha(14)=40*CTR;Beta(14)=0*CTR;K(14)=3;IC(  14)=2; 
h(15)=3134.660;FRQ(15)=9400.255;Alpha(15)=40*CTR;Beta(15)=0*CTR;K(15l=4;IC(  151=1; 
h(16)=3149.529;FRQ(16)=9400.255;Alphan6l=40*CTR;Beta(16)=0*CTR;K(16)=4;IC(16)=2; 
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h(17)=3056.118;FRQ(17> 
h(18)=3073.326;FRQ(18)= 
h(  19)=3 135 .903  ;FRQ(  19)= 
h(20)=3 150.763  ;FRQ(20  )= 
h(21)=3150.763;FRQ(21)= 
h(22)=3 1 63 .28 1  ;FRQ(22)= 

h(23)=3089.195;FRQ(23  )= 
h(24)=3 105.625  ;FRQ(24)= 
h(25)=3159.608;FRQ(25)= 
h(26)=3 172. 122;FRQ(26)= 
h(27)=3 172. 122;FRQ(27> 
h(28)=3 186.2 14;FRQ(28)= 


9400. 1 92;  Alpha(  1 7)=45*CTR;Beta(  1 7> 
9400. 1 92;  Alpha!  1 8)=45  *CTR;Beta(  1 8> 
9400. 1 92;  Alpha!  1 9)=45  *CTR;Beta(  1 9> 
9400. 192;Alpha(20)=45  *CTR;Beta(20)= 
9400. 1 92;  Alpha(2 1  )=45  *CTR;Beta(2 1  > 
9400.192;Alpha(22)=45*CTR;Beta(22)= 


9400.445  ;Alpha(23> 
9400.445 ;  Alpha!24  )= 
9400.445  ;Alpha(25  )= 
9400.445;Alpha(26)= 
9400.445 ;  Alpha(27  )= 
9400.445  ;Alpha(28)= 


:50*CTR;Beta(23)= 
50*CTR;Beta(24)= 
50*CTR;Beta(25)= 
50*CTR;Beta(26  )= 
50*CTR;Beta(27)= 
50*CTR;Beta(28)= 


:0*CTR;K(17)= 

:0*CTR;K(18)= 

:0*CTR;K(19)= 

:0*CTR;K(20)= 

:0*CTR;K(21)= 

:0*CTR;K(22)= 

:0*CTR;K(23)= 

0*CTR;K(24)= 

0*CTR;K(25)= 

:0*CTR;K!26)= 

:0*CTR;K(27)= 

:0*CTR;K(28)= 


=  1  ;IC(  1 73=1 
=  1  ;IC(  1 8)=2; 
=3;IC(19)=1; 
=3  ;IC(20)=2; 
4;IC(21)=1; 
4;IC(22)=2; 

4;IC!23)=1; 

4  ;IC(24)=2; 

4;IC(25)=1; 

4;IC(26)=2; 

4;IC(27)=1 

4;IC(28)=2; 


h(29)=3172.171;FRQ(29)=9400.521;Alpha(29)=60*CTR;Beta(29)=0*CTR;K(29)=l;IC(29)=2; 
h(30)=323 1 .667  ;FRQ( 30)=9400.52 1  ;Alpha(30.)=60*CTR;Beta(30)=0*CTR;K(30)=3  ;IC(30)=1 ; 
h!3 1  )=3244.939;FRQ(31)=9400.521  ;Alpha(3  l)=60*CTR;Beta(3 1)=0*CTR;K(3 1  )=3  ;IC(3 1)=2; 


h!32)=2932.574;FRQ(32)= 

h!33)=2952.130;FRQ(33)= 

h(34)=3139.145;FRQ(34)= 

h!35)=3215.042;FRQ(35)= 

h!36)=3276.065;FRQ(36)= 

h!37)=3287.795;FRQ(37> 


9400.614;Alpha(32)= 

9400.614;Alpha(33)= 

9400.614;Alpha(34)= 

9400.614;Alpha(35)= 

9400.614;Alpha(36)= 

9400.614;Alpha(37)= 


70*CTR;Beta(32)= 

70*CTR;Beta(33)= 

70*CTR;Beta(34)= 

70*CTR;Beta(35)= 

70*CTR;Beta(36)= 

70*CTR;Beta(37)= 


0*CTR;K(32)= 

0*CTR;K(33)= 

0*CTR;K(34)= 

0*CTR;K(35)= 

0*CTR;K(36)= 

0*CTR;K(37)= 


h!38)= 

h(39)= 

h(40)= 

h(41)= 

h!42)= 

h!43)= 

h!44)= 

h!45)= 


:2988.303: 

4005.893: 

4093.348: 

4108.603: 

4267.043: 

4279.369: 

4298.719: 

4309.865: 


FRQ(38)= 

FRQ(39)= 

FRQ(40)= 

FRQ(41)= 

FRQ(42)= 

FRQ(43)= 

FRQ(44)= 

FRQ(45)= 


9400.569: 

9400.569: 

9400.569: 

9400.569: 

9400.569: 

9400.569: 

9400.569: 

9400.569: 


,Alpha(38)= 
;Alpha(39)= 
;Alpha(40)= 
;Alpha(41)= 
,Alpha(42)= 
Alpha(43)= 
;Alpha(44)= 
;Alpha(45  )= 


80*CTR: 

80*CTR: 

80*CTR: 

80*CTR: 

:80*CTR; 

:80*CTR; 

80*CTR: 

80*CTR: 


Beta(38)= 

Beta(39)= 

Beta(40)= 

;Beta(41)= 

Beta(42)= 

Beta(43)= 

;Beta(44)= 

Beta(45)= 


0*CTR; 

0*CTR; 

0*CTR; 

0*CTR; 

0*CTR; 

0*CTR; 

0*CTR: 

0*CTR; 


K(38)= 
K(39)= 
K(40  )= 
K(41  )= 
K(42  )= 
K(43)= 
K!44)= 
K(45)= 


=2;IC(32)=1 

:2;IC(.33')=2; 

4;IC(34)=1; 

d;IC(35)=l; 

4;IC(36)=1; 

4;IC(37)=2; 

:2;IC!381=1 

:2;IC(39)=2; 

4;IC(40)=1; 

4;IC(41)=2; 

d;IC(42)=l 

4;IC(431=2; 

4;IC(44)=1; 

4;IC(45)=2: 


h(46)=3042.965;FRQ(46)=9400.045;Alpha(46)=90*CTR;Beta(46)=0*CTR;K(46)=2;IC(46)=l 

h!47)=3058.819;FRQ(47)=9400.045;Alpha(47)=90*CTR;Beta(47)=0*CTR;K(47)=2;IC(47)=2 

h!48)=3294.133;FRQ(48)=9400.045;Alpha(48)=90*CTR;Beta(48)=0*CTR;K(48)=l;IC(48)=l 

h(49)=3304.680;FRQ(49)=9400.045;Alpha(49)=90*CTR;Beta(49)=0*CTR;K(49)=l;IC(49)=2 


%  Rotation  from  [100]  to  [110],  These  data  were  taken  on  July  24  and  25,  2012. 


h(50)=3042.965;FRQ(50)=9400.045;Alpha(50)=90*CTR;Beta(50)=0*CTR;K(50)=2;IC(50)=l; 

h!51)=3058.819;FRQ(51)=9400.045;Alpha(51)=90*CTR;Beta(51)=0*CTR;K(51)=2;IC(51)=2; 

h!52)=3294.133;FRQ(52)=9400.045;Alpha(52)=90*CTR;Beta(52)=0*CTR;K(52)=l;IC(52)=l; 

h!53)=3304.680;FRQ(53)=9400.045;Alpha(53)=90*CTR;Beta(53)=0*CTR;K(53)=l;IC(53)=2; 

h(54)=3042.308;FRQ(54)=9398.754;Alpha(54)=90*CTR;Beta(54)=15*CTR;K!54)=2;IC(54)=l; 

h(55)=3261.255;FRQ(55)=9398.754;Alpha(55)=90*CTR;Beta[55)=15*CTR;K(55)=5;IC(55)=l; 

h(56)=3272.969;FRQ(56)=9398.754;Alpha(56)=90*CTR;Beta(56)=15*CTR;K(56)=5;IC(56)=2; 

h!57)=3288.831;FRQ(57)=9398.754;Alpha(57)=90*CTR;Beta(57)=15*CTR;K!57)=l;IC(57)=l; 

h!58)=3299.390;FRQ(58)=9398.754;Alpha(58)=90*CTR;Beta(58)=15*CTR;K(58)=l;IC(58)=2; 


h!59)=31 10.997 ;FRQ(59)=9398.742;Alpha(59]=90*CTR;Beta(59)=25*CTR;K(59)=6;IC(59)=2; 
h(60)=3224.281;FRQ(60)=9398.742;Alpha(60)=90*CTR;Beta(60)=25*CTR;K!60)=5;IC(60)=l 
h(61)=3237.190;FRQ(61)=9398.742;Alpha(61)=90*CTR;Beta(61)=25*CTR;K!61)=5;IC(61)=2 
h!62)=3266.529;FRQ(62)=9398.742;Alpha(62)=90*CTR;Beta(62)=25*CTR;K!62)=l;IC(62)=l: 
h(63)=3277.682;FRQ(63)=9398.742;Alpha(63]=90*CTR;Beta(63)=25*CTR;K(63)=l;IC(63)=2; 


h(64)= 

h!65)= 

h(66)= 

h!67)= 

h!68)= 

h(69)= 

h!70)= 


4091.055: 
4105.149: 
4 136.245: 
4150.346: 
4193.166: 
4231.323: 
4243.043: 


FRQ(64)= 
FRQ(65)= 
FRQ(66)= 
FRQ(67  )= 
FRQ(68)= 
FRQ(69)= 
FRQ(70)= 


9398.832: 

9398.832: 

9398.832: 

9398.832: 

9398.832: 

9398.832: 

9398.832: 


;Alpha(64)= 
Alpha(65)= 
Alpha(66)= 
,Alpha(67)= 
;Alpha(68  )= 
;  Alpha(69  )= 
;Alpha(70)= 


90*CTR: 

90*CTR: 

90*CTR: 

90*CTR: 

90*CTR: 

90*CTR: 

90*CTR: 


Beta(64)= 

Beta(65)= 

Beta(66)= 

Beta(67)= 

Beta(68)= 

Beta(69)= 

;Beta(70)= 


35*CTR: 

35*CTR: 

35*CTR: 

35*CTR: 

35*CTR: 

35*CTR: 

35*CTR: 


K(64)= 
K!65)= 
K!66)= 
K(67)= 
K(68  )= 
K(69)= 
K(70)= 


2;IC(64)=1 

2;IC(65)=2; 

6;IC(66)=1: 

6;IC(67)=2 

5;IC(68)=2; 

1;IC(69)=1 

1;IC(70)=2 


h(71)=3132.114;FRQ(71)=9399.142;Alpha(71)=90*CTR;Beta(71)=45*CTR;K!71)=2;IC(71)=l: 

h!72)=3146.786;FRQ(72)=9399.142;Alpha(72)=90*CTR;Beta(72)=45*CTR;K!72)=2;IC(72)=2; 

h(73)=3185.534;FRQ(73)=9399.142;Alpha(73)=90*CTR;Beta(73)=45*CTR;K(73)=l;IC(73)=l; 
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h(74)=3197.270;FRQ(74)=9399.142;Alpha(74)=90*CTR;Beta(74)=45*CTR;K(74)=l;IC(74)=2; 


datapoints  =  length(h); 

for  nn=l:  datapoints 
HH  =  h(nn); 
alpha  =  Alpha(nn); 
beta  =  Beta(nn); 
k  =  K(nn); 
ic  =  IC(nn); 

R(l,l)  =  cos(alpha)*cos(beta); 
R(l,2)  -  -sin(beta); 

R(l,3)  =  sin(alpha)*cos(beta); 
R(2,l)  =  cos(alpha)*sin(beta); 
R(2,2)  =  cos(beta); 

R(2,3)  =  sin(alpha)*sin(beta); 
R(3,l)  =  -sin(alpha); 

R(3,2)  =  0; 

R(3,3)  =  cos(alpha); 


if  k==  1 


Rif  1,1)  =  1;R1(1,2)  =  0;R1(1,3)  =  0; 
Rl(2,l)  =  0;R1(2,2)  =  1;R1(2,3)  =  0; 
Rl(3,l)  =  0:R1(3,2)  =  0;R1(3,3)  =  1; 


elseif  k==2 


Rif  1,1)  =  0:R1(1,2)  =  1;R1(  1,3)  =  0; 
Rl(2,l)  =  -1  ;R1(2,2)  =  0;R1(2,3)  =  0; 
Rl(3,l)  =  0:R1(3,2)  =  0;R1(3,3)  =  1; 


elseif  k==3 


Rif  1,1)  =  -1;R1(1,2)  =  0;R1(1,3)  =  0; 
Rl(2,l)  =  0:R1(2,2)  =  -1;R1(2,3)  =  0; 
Rl(3,l)  =  0:R1(3,2)  =  0;R1(3,3)  =  1; 


elseif  k==4 


Rif  1,1)  =  0:R1(1,2)  =  -1;R1(1,3)  =  0; 
Rl(2,l)  =  1:R1(2,2)  =  0;R1(2,3)  =  0; 
R1(3,1)  =  0:R1(3,2)  =  0;R1(3,3)=  1: 


elseif  k==5 


Rif  1,1)  =  1;R1(1,2)  =  0;R1(1,3)  =  0; 
Rl(2,l)  =  0:R1(2,2)  =  -1;R1(2,3)  =  0; 
Rl(3,l)  =  0:R1(3,2)  =  0;R1(3,3)  =  -1; 


elseif  k==6 


Rif  1,1)  =  0:R1(1,2)  =  1;R1(  1,3)  =  0; 
Rl(2,l)=  1:R1(2,21  =  0;R1(2,3)  =  0; 
Rl(3,l)  =  0:R1(3,2)  =  0;R1(3,3)  =  -1; 


elseif  k==7 


Rif  1,1)  =  -1;R1(1,2)  =  0;R1(1,3)  =  0; 
Rl(2,l)  =  0:R1(2,21  =  1;R1(2,3)  =  0; 
Rl(3,l)  =  0:R1(3,2)  =  0;R1(3,3)  =  -1; 


elseif  k==8 


Rif  1,1)  =  0;R1(1,2)  =  -1;R1(1,3)  =  0; 
Rl(2,l)  =  -1;R1(2,2)  =  0;R1(2,3)  =  0; 
Rl(3,l)  =  0:R1(3,2)  =  0;R1(3,3)  =  -1; 

end 
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RT  =  R1  *  R; 
TG  =  G  *  RT; 
TH  =  H  *  RT; 


W1=B*HH*(P(1)*TG(1,1)*TG(1,3)+P(2)*TG(2,1)*TG(2,3)+P(3)*TG(3,1)*TG(3,3)); 

W2=B*HH*(P(1)*TG(  1 ,2)*TG(  1,3)+P(2)*TG(2,2)*TG(2>3)+P(3)*TG(3,2)*TG(3,3)); 

W3=B*HH*(P(1)*TG(1,3)*TG(13)+P(2)*TG(2,3)*TG(2,3)+P(3)*TG(3,3)*TG(3,3)); 

W4=P(7)*TH(1,1)*TH(1,1)+P(8)*TH(2,1)*TH(2,1)+P(9)*TH(3,1)*TH(3,1); 

W5=P(7)*TH(1,1)*TH(1,2)+P(8)*TH(2,1)*TH(2,2)+P(9)*TH(3,1)*TH(3,2); 

W6=P(7)*TH(1 , 1  )*TH(  1,3)+P(8)*TH(2, 1  )*TH(2,3)+P(9)*TH(3, 1)*TH(3,3); 

W7=P(7)*TH(1,2)*TH(1,2)+P(8)*TH(2,2)*TH(2,2)+P(9)*TH(3,2)*TH(3,2); 

W8=P(7)*TH(  1 ,2)*TH(  1 3)+P(8)*TH(2,2)*TH(2,3)+P(9)*TH(3,2)*TH(3,3); 

W9=P(7)*TH(1,3)*TH(1,3)+P(8)*TH(2,3)*TH(2,3)+P(9)*TH(3,3)*TH(3,3); 


Q1  =  0.5*(Wl+i*W2); 

Q2  =  0.25*(W4-W7)+0.5*i*W5; 

Q3  =  0.25*(W4+W7); 

Q4  =  0.5*(W6+i*W8); 

%  Ham  is  the  matrix  representing  the  spin-Hamiltonian: 


Ham  =  zeros(4); 

Ham(l.l)  =  0.5*W3  +  0.25*W9  -  0.5*gbn*HH; 
Ham(2,2)  =  0.5*W3  -  0.25*W9  +  0.5*gbn*HH; 
Ham(3,3)  =  -0.5*W3  -  0.25*W9  -  0.5*gbn*HH; 
Ham(4,4)  =  -0.5*W3  +  0.25*W9  +  0.5*gbn*HH; 
Ham(2,l)  =  0.5*Q4; 

Ham(3.1)  =  Ql  +0.5*Q4; 

Ham(3,2)  =  Q3; 

Ham(4,l)  =  Q2; 

Ham(4,2)  =  Ql  -  0.5*Q4; 

Ham(4,3)  =  -0.5*Q4; 

Ham(l,2)  =  conj(Ham(2,l)); 

Ham(l,3)  =  conj(Ham(3,l)); 

Ham(l,4)  =  conj(Ham(4,l)); 

Ham(2,3)  =  conj(Ham(3,2)); 

Ham(2,4)  =  conj(Ham(4,2)); 

Ham(3,4)  =  conj(Ham(4,3)); 

EE  =  sort(real(eig(Ham))); 

if  ic==l 

freq(nn)  =  abs(EE(4)-EE(l)); 


elseif  ic-=2 

freq(nn)  =  abs(EE(3)-EE(2)); 
end 
end 

summ=0; 

for  ii=l  :datapoints 

sumrn  =  summ  +  (FRQ(ii)-freq(ii))A2; 
end 


%  EPR_linepositions_gandA_P 

% 

%  Calculation  of  EPR  linepositions  for  perturbed  Ag2+  hole  center  in  Li2B407 
%  (This  program  is  titled  EPR_linepositions_gandA_P.m) 

%  This  program  calculates  the  EPR  line  positions  as  a  function  of  angle 
%  for  the  x-ray-induced  perturbed  Ag2+  hole  center  in  LTB.  It  is  uses 
%  a  subroutine  entitled  AgEnergyLevels_P.  Both  the  g  matrix  and  A  matrix 
%  are  included  (i.e.,  a  4x4  spin-Hamiltonian  matrix  is  diagonalized). 

clear  all 

%  Constants: 
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h  =  6.626069; 

B  =  9.274009/h; 
CTR  =  pi/180; 
FREQQ  =  9400; 
InitialHH  =  3300: 


%  Planck's  constant 

%  Bohr  magneton  divided  by  Planck's  constant 
%  Conversion  constant,  degrees  to  radians 
%  Measured  microwave  frequency  (in  MHz) 
%  Initial  guess  for  line  positions  (in  gauss) 


gbn  =  -0.00018637;  %  averaged  value  of  gn*bn  for  107Ag  and  109Ag 


%  Spin-Hamiltonian  parameters 

%  Six  for  the  g  matrix  (three  principal  values  and  three  angles). 
%  Six  for  the  A  matrix  (three  principal  values  and  three  angles). 


P(l)  =  2.0319; 

P(2)  =  2.1212; 

P(3)  =  2.3604; 

P(4)  =  35.5*CTR; 
P(5)  =  9.7*CTR; 
P(6)  =  -9.6*CTR; 
P(7)  =  28.18; 

P(8)  =  35.81; 

P(9)  =  65.12; 

P(  10)  =  35.2*CTR; 
P(1 1)  =  25.9*CTR; 
P(  1 2)  =  -38.7*CTR; 


%  G  is  the  3x3  rotation  matrix  which  takes  the  principal 
%  axes  of  the  g  matrix  into  the  crystal  coordinate  system. 

%  Taken  from  "Classical  Mechanics,  2nd  ed."  by  Goldstein,  pp.  146-147. 

%  H  is  the  3x3  rotation  matrix  which  takes  the  principal 
%  axes  of  the  A  matrix  into  the  crystal  coordinate  system. 

%  R1  is  the  3x3  rotation  matrix  which  takes  sites  1  through  8 
%  into  the  crystal  coordinate  system. 

%  R  is  the  3x3  rotation  matrix  which  takes  the  crystal  coordinate 
%  system  into  the  magnetic  field  coordinate  system. 

G(  1,1)  =  cos(P(6))*cos(P(5))  -  cos(P(4))*sin(P(5 ))*sin(P(6)); 

G(l,2)  =  cos(P(6)  )*sin(P(5))  +  cos(P(4))*cos(P(5))*sin(P(6)); 

G(l,3)  =  sin(P(6))*sin(P(4)); 

G(2,l)  =  -sin(P(6))*cos(P(5))  -  cos(P(4))*sin(P(5))*cos(P(6)); 

G(2,2)  =  -sin(P(6))*sin(P(5))  +  cos(P(4))*cos(P(5))*cos(P(6)); 

G(2,3)  =  cos(P(6))*sin(P(4)); 

G(3,l)  =  sin(P(4))*sin(P(5)); 

G(3,2)  =  -sin(P(4 )  )*cos(P(5)); 

G(3,3)  =  cos(P(4)); 

H(l,l)  =  cos(P(  1 2))*cos(P(  1 1 ))  -  cos(P(  1 0))*sin(P(  1 1  ))*sin(P(  12)); 
H(l,2)  =  cos(P(12))*sin(P(l  1))  +  cos(P(10))*cos(P(l  l))*sin(P(12)); 
H(l,3)  =  sin(P(12))*sin(P(10)); 

H(2,l)  =  -sin(P(12))*cos(P(l  1))  -  cos(P(10))*sin(P(ll))*cos(P(12)); 
H(2,2)  =  -sin(P(12))*sin(P(l  1 ))  +  cos(P(10))*cos(P(l  l))*cos(P(12)); 
H(2,3)  =  cos(P(  1 2  ))*sin(P(  10)); 

H(3,l )  =  sin(P(10))*sin(P(l  1)); 

H(3,2)  =  -sin(P(10))*cos(P(l  1 )); 

H(3,3)  =  cos(P(10)); 


%  Select  a  plane  of  rotation  before  running  the  program. 

%  Plane  =  1  corresponds  to  rotation  from  a  to  a. 

%  Plane  =  2  corresponds  to  rotation  from  c  to  a. 

%  Plane  =  3  corresponds  to  rotation  from  c  to  [1 10]. 

%  Enter  the  number  below  for  the  plane  to  be  used. 

forPlane=l  %  Rotation  from  a  to  a. 

for  n=l  :91  %  Rotation  increment  is  one  degree. 

Alpha  =  90*CTR; 
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Beta  =  (n-l)*CTR; 
N(n)  =  n-1; 


Rf  1,1)  =  cos(Alpha)*cos(Beta); 
R(l,2)  =  -sin(Beta); 

R(l,3)  =  sin(Alpha)*cos(Beta); 
R(2,l)  =  cos(Alpha)*sin(Beta); 
R(2,2)  =  cos(Beta); 

R(2,3)  =  sin(Alpha)*sin(Beta); 
R(3,l)  =  -sin(Alpha); 

R(3,2)  =  0; 

R(3,3)  =  cos(Alpha); 

for  k=l:8 


if  k==  1 


Rl(l,l)  =  1;R1(1,2)  =  0;R1(1,3)  =  0; 
Rl(2,l)  =  0;R1(2,2)  =  1;R1(2,3)  =  0; 
R1(3,1)  =  0;R1(3,2)  =  0;R1(3,3)  =  1; 


elseif  k==2 

Rl(l.l)  =  0;R1(1.2)  =  1  ;R  1(1,3)  =  0; 
Rl(2,l)  =  -1;R  1(2,2)  =  0;R1(2,3)  =  0; 
Rl(3,l)  =  0;R1(3,2)  =  0;R1(3,3)  =  1; 
elseif  k==3 

Rl(l.l)  =  - 1  ;R1  (1 ,2)  =  0;R1(1,3)  =  0; 
Rl(2,l)  =  0;R1(2,2)  =  -1;R1(2,3)  =  0; 
Rl(3,l)  =  0;R1(3.2)  =  0;R1(3,3)  =  1; 
elseif  k==4 

Rl(l,l)  =  0;R1(1,2)  =  - 1  ;R1  ( 1 ,3)  =  0; 
Rl(2,l)  =  1;R1(2,2)  =  0;R1(2,3)  =  0; 
R1(3,1)  =  0;R1(3,2)  =  0;R1(3,3)  =  1; 
elseif  k==5 

Rl(l,l)  =  1;R1(1,2)  =  0;R1(1,3)  =  0; 
Rl(2,l)  =  0;R1(2,2)  =  -1;R1(2,3)  =  0; 
Rl(3,l)  =  0;R1(3,2)  =  0;R1(3,3)  =  -1; 

elseif  k==6 

Rl(l,l)  =  0;R1(1,2)  =  1;R1(1,3)  =  0; 
Rl(2,l)  =  1;R1(2,2)  =  0;R1(2,3)  =  0; 
Rl(3,l)  =  0;R1(3.2)  =  0;R1(3,3)  =  -1; 
elseif  k==7 

Rl(l,l)  =  -1;R1(1,2)  =  0;R1  ( 1 ,3)  =  0; 
Rl(2,l)  =  0;R1(2,2)  =  1;R1(2,3)  =  0; 
Rl(3,l)  =  0;R1(3.2)  =  0;R1(3,3)  =  -1: 


elseif  k==8 

Rl(l,l)  =  0;R1(1,2)  =  - 1  ;R1  (1 ,3)  =  0 
Rl(2,l)  =  -1;R1(2,2)  =  0;R1(2,3)  =  0 
Rl(3,l)  =  0;R1(3,2)  =  0;R1(3,3)  =  -1 

end 

RT  =  R1  *  R: 

TG  =  G  *  RT; 

TH  =  H  *  RT; 


for  ic= 1 : 2  %  Specifies  the  particular  EPR  transition, 

if  ic==l 

HH  =  InitialHH; 

EE  =  AgEnergyLevels_P(TG,TH,HH,P,B,gbn); 

FREQ  =  EE(4)-EE(1); 
while  abs(FREQQ-FREQ)>  1 ; 

HH  =  HH  *  (FREQQ/FREQ); 

EE  =  AgEnergyLevels_P(TG,TH,HH,P,B,gbn); 
FREQ  =  EE(4)-EE(1 ); 
end 

MagneticField(n,k)  =  HH; 
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elseif  ic==2 

HH  =  InitialHH; 

EE  =  AgEnergyLevels_P(TG,TH,HH,P,B,gbn); 
FREQ  =  EE(3)-EE(2); 

while  abs(FREQQ-FREQ)>  1 ; 

HH  =  HH  *  (FREQQ/FREQ); 

EE  =  AgEnergyLevels_P(TG,TH,HH,P,B,gbn>; 
FREQ  =  EE(3)-EE(2); 
end 

MagneticField(n,k+8)  =  HH; 


end 

end 

end 

end 

for  Plane=2  %  Rotation  from  c  to  a. 

for  n=l  :91  %  Rotation  increment  is  one  degree. 

Alpha  =  (n-l)*CTR; 

Beta  =  0*CTR; 

N(n)  =  n-1; 


Rf  1,1)  =  cos(Alpha)*cos(Beta); 
R(l,2)  =  -sin(Beta); 

R(l,3)  =  sin(Alpha)*cos(Beta); 
R(2,l)  =  cos(Alpha;*sin(Beta); 
R(2,2)  =  cos(Beta); 

R(2,3)  =  sin(Alpha)*sin(Beta); 
R(3,l)  -  -sin(Alpha); 

R(3,2)  =  0; 

R(3,3)  =  cos(Alpha); 

for  k=l:8 


if  k==l 

R1(1,1)=1;R1(1,2)  =  0;R1(1,3)  =  0; 
Rl(2,l)  =  0;R1(2.2)  =  1;R1(2,3)  =  0; 
Rl(3,l)  =  0;R1(3,2)  =  0;R1(3,3)  =  1; 
elseif  k==2 

Rl(l,l)  =  0:R1(1,2)  =  1  ;R1(  1,3)  =  0; 
Rl(2,l)  =  -1  ;R1(2.2)  =  0;R1(2,3)  =  0; 
Rl(3,l)  =  0;R1(3,2)  =  0;R1(3,3)  =  1; 
elseif  k==3 

Rif  1,1)  =  -1;R1(1,2)  =  0;R1(1,3)  =  0; 
Rl(2,l)  =  0;R1(2.2)  =  -1;R1(2,3)  =  0; 
Rl(3,l)  =  0;R1(3,2)  =  0;R1(3,3)  =  1; 
elseif  k==4 

Rl(l,l)  =  0:R1(1,2)  =  - 1  ;R1  ( 1 ,3)  =  0; 
Rl(2,l)  =  1:R1(2,2)  =  0;R1(2,3)  =  0; 
Rl(3,l)  =  0;R1(3.2)  =  0;R1  (3,3)  =  1; 
elseif  k==5 

Rl(l,l)  =  1  ;R1  ( 1 .2)  =  0;R1(1,3)  =  0; 
Rl(2,l)  =  0;R1(2,2)  =  -1;R1(2,3)  =  0; 
Rl(3,l)  =  0;R1(3,2)  =  0;R1(3,3)  =  -1; 

elseif  k==6 

Rif  1,1)  =  0;R1(1.2)  =  1;R1(1,3)  =  0; 
Rl(2,l)  =  1;R1(2,2)  =  0;R1(2,3)  =  0; 
Rl(3,l)  =  0;R1(3,2)  =  0;R1(3,3)  =  -1; 
elseif  k==7 

Rl(l,l)  =  -1;R1(1,2)  =  0;R1  ( 1 ,3)  =  0; 
Rl(2,l)  =  0;R1(2.2)  =  1;R1(2,3)  =  0; 
Rl(3,l)  =  0;R1(3.2)  =  0;R1(3,3)  =  -1; 
elseif  k==8 

Rl(l,l)  =  0:R1(1,2)  =  - 1  ;R1  ( 1 ,3)  =  0; 
Rl(2,l)  =  -1  ;R1(2,2)  =  0;R1(2,3)  =  0; 
Rl(3,l)  =  0;R1(3,2)  =  0;R1(3,3)  =  -1; 
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end 


RT  =  R1  *  R: 
TG  =  G  *  RT; 
TH  =  H  *  RT; 


for  ic=  1 ; 2  %  Specifies  the  particular  EPR  transition, 

if  ic==l 

HH  =  InitialHH; 

EE  =  AgEnergyLevels_P(TG,TH,HH,P,B,gbn); 

FREQ  =  EE(4)-EE(1); 

while  abs(FREQQ-FREQ)>  1 ; 

HH  =  HH  *  (FREQQ/FREQ); 

EE  =  AgEnergyLevels_P(TG.TH,HH,P,B.gbn); 
FREQ  =  EE(4)-EE(1); 
end 

MagneticField2(n,k)  =  HH; 

elseif  ic==2 

HH  =  InitialHH; 

EE  =  AgEnergyLevels_P(TG,TH,HH,P,B,gbn); 

FREQ  =  EE(3)-EE(2); 

while  abs(FREQQ-FREQ)>  1 ; 

HH  =  HH  *  (FREQQ/FREQ); 

EE  =  AgEnergyLevels_P(TG,TH,HH,P,B,gbn); 
FREQ  =  EE(3)-EE(2); 
end 

MagneticField2(n,k+8)  =  HH; 


end 

end 

end 

end 

end 


for  Plane=l 


h(  1  )=3042.965  ;Beta(  1  )=0; 
h(2)=3058.819;Beta(2)=0; 
h(3)=3294.133;Beta(3)=0; 
h(4)=3304.680;Beta(4)=0; 


h(5)=3042.308;Beta(5)=15; 

h(6)=3261.255;Beta(6)=15; 

h(7)=3272.969;Beta(7)=15; 

h(8)=3288.831;Beta(8)=15; 

h(9)=3299.390;Beta(9)=15; 

h(10)=3 1 1 0.997  ;Beta(  1 0)=25 ; 
h(  1 1)=3224.281  ;Beta(l  1)=25; 
h(12)=3237.190;Beta(12)=25; 
h(13)=3266.529;Beta(13)=25; 
h(14)=3277.682;Beta(14)=25; 

h(  1 5)=309 1 .055  ;Beta(  1 5  )=35 ; 
h(16)=3105.149;Beta(16)=35; 
h(17)=3 136.245  ;Beta(17)=35; 
h(  18)=3 150.346 ;Beta(18)=35 ; 
h(19)=3193.166;Beta(19)=35; 
h(20)=3231.323;Beta(20)=35; 
h(2 1  )=3243 .043  ;Beta(2 1  )=35 ; 

h(22)=3132.1 14;Beta(22)=45; 
h(23)=3 146.786;Beta(23)=45 ; 
h(24)=3185.534;Beta(24)=45; 
h(25)=3 197.270;Beta(25  )=45; 
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h(26)=3042.965;Beta(26)=90; 

h(27)=3058.819;Beta(27)=90; 

h(28)=3294.133;Beta(28)=90; 

h(29)=3304.680;Beta(29)=90; 

h(30)=3042.308;Beta(30)=75; 
h(3 1  )=326 1 .255  ;Beta(3 1  )=75 ; 
h(32)=3272.969;Beta(32)=75 ; 
h(33)=3288.83 1  ;Beta(33)=75; 
h(34)=3299.390;Beta(34)=75; 

h(35)=3 1 10.997  ;Beta(35)=65; 
h(36)=3224.28 1  ;Beta(36)=65 ; 
h(37)=3237. 190;Beta(37  )=65 ; 
h(38)=3266.529;Beta(38)=65 ; 
h(39)=3277.682;Beta(39)=65 ; 

h(40)=309 1 .055  ;Beta(40)=55 ; 
h(4 1 )=3 105 . 1 49;Beta(4 1  )=55 ; 
h(42)=3 1 36.245  ;Beta(42)=55 ; 
h(43)=3 150.346 ;Beta(43  )=55 ; 
h(44)=3193.166;Beta(44)=55; 
h(45)=3231.323;Beta(45)=55; 
h(46)=3243 .043  ;Beta(46)=55 ; 

plot(N+90,MagneticField,' l,Beta+90,h,'oVcolor',,b') 
ylabelf'Magnetic  Field  (gauss)') 
xlabel('Angle  (degrees)’) 
end 

for  Plane=2 

h(  1  )=293 1 .543  ;Alpha(  1  )=0; 
h(2)=295 1 .301  ;Alpha(2)=0; 

h(3)=2889.426;Alpha(3)=10; 
h(4)=29 1 1 ,328;Alpha(4)=10; 
h(5)=293 1 .28 1  ;Alpha(5)=10; 
h(6)=2981.737;Alpha(6)=10; 

h(7)=2858.8 1 8;Alpha(7)=20; 
h(8)=2945 .435  ;Alpha(8)=20; 
h(9)=3008.362;Alpha(9)=20; 
h(10)=3057.329;  Alpha(  10)=20; 

h(  1 1  )=3094. 1 74;Alpha(  1 1  )=30; 
h(  12)=3 1 1 0.0 1 6;Alpha(  1 2)=30; 

h(13)=3047.808;Alpha(13  )=40; 
h(  14)=3 120.576;Alpha(  14)=40; 
h(  1 5)=3 1 34.660;  Alpha(  1 5  )=40; 
h(16)=3149.529;Alpha(16)=40; 

h(17)=3056.118;Alpha(17)=45; 
h(  1 8)=3073 .326;  Alpha(  1 8  )=45 ; 
h(  19)=3 1 35 .903 ;  Alpha(  1 9)=45 ; 
h(20)=3150.763;Alpha(20)=45; 
h(2 1  )=3 150.763  ;Alpha(2 1  )=45 ; 
h(22)=3 163.281  ;Alpha(22)=45; 

h(23)=3089.195;Alpha(23)=50; 
h(24)=3105.625;Alpha(24)=50; 
h(25)=3 159.608;Alpha(25)=50; 
h(26)=3172.122;Alpha(26)=50; 
h(27)=3172.122;Alpha(27)=50; 
h(28)=3 186.2 1 4;Alpha(28)=50; 

h(29)=3 172.171  ;Alpha(29)=60; 
h(30)=323 1 .667  ;Alpha(30)=60; 
h(3 1  )=3244.939;Alpha(3 1  )=60; 
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h(32)=2932.574;Alpha(32)=70; 
h(33)=2952.130;Alpha(33)=70; 
h(34)=3139.145;Alpha(34)=70; 
h(35)=32 1 5 .042  ;Alpha(35)=70; 
h(36)=3276.065;Alpha(36)=70; 
h(37)=3287.795;Alpha(37)=70; 

h(38)=2988.303;Alpha(38)=80; 
h(39)=3005.893;Alpha(39)=80; 
h(40)=3093 .348  ;Alpha(40)=80; 
h(4 1  )=3 108.603  ;Alpha(4 1  )=80; 
h(42)=3267.043;Alpha(42)=80; 
h(43)=3279.369;Alpha(43  )=80; 
h(44)=3298.719;Alpha(44)=80; 
h(45)=3309.865;Alpha(45)=80; 

h(46)=3042.965;Alpha(46)=90; 
h(47)=3058.819;Alpha(47  )=90; 
h(48)=3294. 1 33;  Alpha(48  )=90; 
h(49)=3304.680;Alpha(49)=90; 


X(l)=30; 

Y(l)=2841; 

X(2)=30; 

Y(2)=2867; 

X(3)=40; 

Y(3)=2838; 

X(4)=40; 

Y(4)=2863; 

X(5)=45; 

Y(5)=2841; 

X(6)=45; 

Y(6)=2867; 

X(7)=50; 

Y(7)=2856; 

X(8)=50; 

Y(8)=2882; 

X(9)=60; 

Y(9)=2903; 

X(10)=60; 

Y(10)=2924: 


hold  on 

plot(N,MagneticField2,'  '.Alpha, h, 'o', 'color', 'b') 
hold  on 

plot(X,Y,'rsVMarkerSize',5,'color','k') 

ylabel('Magnetic  Field  (gauss)') 

xlabelCAngle  (degrees)’) 

end 

end 
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